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ABSTRACT 


A large variety of physical and chemical systems are characterized by the 
repetition of identical units. The periodic structure of these systems allow their 
handling by similar mathematical methods. It is the purpose of this thesis to 
present certain techniques for the evaluation of analytic functions of matrices 
associated with such systems. The method consists of Poisson-type transforma- 
tions of finite sums, which lead to rapidly converging expressions for the problems 
considered. Examples related to lattice dynamics and molecular orbital theory 
are discussed. Certain restricted sums of perturbation theory are calculated 
exactly by a method of independent interest. An extensive list of matrices arising 
from diverse periodic systems, together with their eigenvalues and eigenvectors, 
is presented. This extends and generalizes results in the literature. 




iii 


TABLE OF CONTENTS 


Chapter Page 

INTRODUCTION 1 

I. DYNAMICS OF A CRYSTAL LATTICE 5 

1.1 One-Dimensional Lattices 13 

1.2 Two-Dimensional Lattices 25 

1.3 Three-Dimensional Lattices 41 

II. OTHER APPLICATIONS 53 

IH. METHOD OF CALCULATION AND SPECIFIC RESULTS 61 

3.1 One-Dimensional Sums 64 

3.2 Two-Dimensional Sums 92 

3.3 Three-Dimensional Sums 107 

IV. LATTICES WITH DEFECTS AND PERTURBATION THEORY .... 117 

APPENDIX A. EQUATIONS OF MOTION FOR RECTANGULAR 

LATTICES 133 

APPENDIX B. EVALUATION OF AN INTEGRAL 141 

APPENDIX C. FOURIER SERIES FOR A S -FUNCTION 149 

APPENDIX D. PARTIAL FRACTION DECOMPOSITION AND 

SUMMATION OF FINITE SUMS 157 

APPENDIX E. MATRICES ASSOCIATED WITH PERIODIC 

SYSTEMS 169 

ACKNOWLEDGMENT 245 

SELECTED BIBLIOGRAPHY 247 


v 



INTRODUCTION 


There are a large variety of phenomena in physics and theoretical chemistry 
which are characterized by the repetition of identical units. This type of model 
arises in the study of lattice structure as in lattice dynamics, exciton theory, 
polymer chains and other similar situations. Many of the equations which appear 
in the analysis of such periodic systems have similar features, and allow sim- 
plifications due to the periodicities involved. 

The first one to call attention to the class of mathematical problems which 
arise from periodic phenomena was Rutherford [1,2], who calculated the eigen- 
vectors and eigenvalues of particular continuant and circulant matrices useful 
mainly in the solution of one dimensional problems. Some of this work was con- 
tinued in a mathematical context by Egervary [3]. Of course, particular problems 
involving periodic structures have been solved from Bernouilli onwards. A good 
deal of this work is summarized in the books by Brillouin [4] and by Parodi and 
Brillouin [5], Recently there has been some interest in this genre of problems, 
as it arises from the analysis of molecular systems, [6,7]. 

It is the purpose of this thesis to present certain techniques for calculating 
analytic functions of circulant and continuant matrices, appropriate to the solu- 
tion of one- and higher-dimensional problems. In doing so we shall also extend 
some of the results of Rutherford, particularly in regards to matrices character- 
istic of multi-atomic problems and lattices with defects. 
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In Chapter I we introduce the basic equations describing the vibrational 
motion of crystal lattices in the harmonic approximation and the matrices as- 
sociated with particular atomic configurations are derived. These are a good 
example of the type of problem arising in studies of the physics of periodic units. 
The main model employed is that of a rectangular lattice with a variety of boundary 
conditions. The ground work needed for the elaboration of these systems can be 
useful for the treatment of models with other symmetries. The statistical me- 
chanics of lattices is shown to lead naturally to the concept of function of a 
matrix and the need for an explicit evaluation of its elements. 

Chapter II describes problems selected from several fields (theoretical 
chemistry, molecular physics, etc.) which give rise to types of matrices similar 
to those discussed in Chapter I. 

Chapter III presents a technique for calculating analytic functions of the 
particular matrices introduced in the previous chapters. The task at hand is 
shown to be the evaluation of certain finite sums. This is accomplished by a 
summation method somewhat similar to the Poisson sum formula. The formalism 
thus developed is then applied to several examples in one-, two- and three- 
dimensions which lead to specific analytic results. The method provides a scheme 
for numerical work in cases which are analytically intractable. 

The systems dealt with in foregoing chapters possessed unperturbed perio- 
dicities. When perturbations are introduced (such as isotopic impurities, holes, 
etc.) the treatment has to be modified. Montroll and coworkers have shown [8,9] 
how to evaluate sums of analytic functions of the unknown eigenfrequencies of 
such systems by using matric Green's functions of the unperturbed systems. 

While the same approach could be adopted here by employing the results of 
Chapter HI, one can instead find approximate eigenfrequencies and then evaluate 
the appropriate sums of these. This procedure is illustrated in Chapter IV 
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where closed form expressions are found for the restricted sums of ordinary 
perturbation theory , by using a method of some intrinsic interest. 

In order to increase the coherence of presentation some of the lengthier 
derivations have been relegated to appendices. In addition, a special appendix is 
devoted to the listing of matrices which appear most frequently in applications. 
This appendix contains generalizations and extensions of Rutherford’s work 
[1,2] and each matrix listed in it is accompanied, whenever possible, by its 
eigenvalues, eigenvectors and the closed form of its characteristic polynomial. 
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CHAPTER I 


DYNAMICS OF A CRYSTAL LATTICE 

A solid is by definition a group of atoms arranged in a regular array (lattice), 
which perform small oscillations about their equilibrium positions. It is usually 
assumed that the forces acting on each particle in the lattice can be derived from 
a potential. 

In most cases the exact form of this potential is not known and arbitrary 
constants have to be introduced in the theory. On the other hand, stability con- 
siderations show that the potential energy of the lattice can be expanded in a 
Taylor series about the equilibrium positions, at least for temperatures below 
those at which a state transition takes place. If only the first non-zero term is 
retained in this expansion, the potential energy becomes a quadratic form in 
the displacements of the atoms. This procedure is known as the harmonic approxi- 
mation, the nomenclature referring to the fact that here the solid is regarded as 
a set of coupled harmonic oscillators. The harmonic approximation describes 
adequately many of the vibrational properties of a solid although it cannot explain 
characteristics such as thermal expansion, temperature variation of the elastic 
constants, etc. For the treatment of these phenomena higher order terms in the 
series expansion have to be included in the potential energy. In this work all the 
systems considered will be treated in the harmonic approximation. 



The general theory of lattice vibrations is discussed in detail in several 
books [10, 11] ; here only an outline of the theory will be given and specialization 
will be made to those models for which the methods of calculation described in 
Chapter III apply. 

To write down the Hamiltonian for the lattice we assume for definiteness that 
there are N unit cells each containing n particles. We define r°(^) and r 
to be the position vector at equilibrium and the actual position vector, respectively, 
of the K-th particle in the l-th unit cell, both related to a suitable origin of co- 
ordinates. The vectors r°^j represent the lattice sites and can be written 
also in the form 


r °(>t) = r °^ + r °^ ^ 

with 

r 0 ^) = -Cj a x + ■i 2 a 2 + a 3 (2) 

where the a. are three noncoplanar vectors called the primitive translation 
vectors of the crystal and the l. are three integers, positive, negative or zero. 
The position vector of the x-th particle in the unit cell with respect to some 
origin there, is denoted above by r° (x) . All crystals the unit cells of which 
contain only one particle are called Bravais crystals — other crystals are called 
nonprimitive. 


In the harmonic approximation the potential energy of the lattice is given by 
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The zero subscript above denotes evaluation at the appropriate equilibrium 
position. By definition 




(i) • 


= 0 


If we define the displacements from equilibrium by 


■0 


( 4 ) 


(5) 


and the (a, /3) component of the tensor | by 


<D. 







a, /3 = 1, 2, 3 


( 6 ) 


then the potential energy can be written as follows: 


-••i? CHIN::) 


(7) 


'CHa 

Vk'p 


This is a quadratic surface in 3Nn dimensions and one can investigate its in- 
variance properties under a variety of transformations. This will not be done 
here as it is excellently presented in [ 11 ] . We want to remark that while 
important, these invariance properties pertain mostly to infinite lattices and 
lattices with boundaries which have been removed by use of the Born-von Karman 
cyclic constraints. 


The coefficient of $> a/8 is the force exerted in the a-direction on the par- 
ticle at when the particle at r°^'j is displaced a unit distance in the 


fi -direction, and from eq. (6) we see that 
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•-(13 --Cl) 


( 8 ) 


If the potential energy consists only of two-body interactions, with each pair 
of particles interacting via a potential function <£ KK » (r) which depends only on 
the magnitude of their separation, the atomic force constants will take the form: 



The number of particle species is at most equal to the number of particles 
in the unit cell. Therefore M K , x = 1, . . . ,n , will denote the mass of the x -th 
particle in the unit cell, and we can write the equations of motion 


»Mt)- £* 4 ::)^:) 


( 10 ) 


•£'x'/ 3 


There are two different approaches usually adopted for solving these equations: 

1. The Plane-Wave Method 

If in eqs.(10) one proceeds to the limit of an elastic continuum, one obtains 
a wave equation for the displacement vector of an arbitrary point in the medium. 
An expansion of the displacement in plane waves of the type exp {2 7rik'r-io>t} 
leads in most cases to a solution of the problem, the amplitudes being determined 
from the boundary conditions. 

This treatment of the wave equation suggests a similar procedure for the 
discrete crystal lattice. If we make the substitution 

) = — - — v a (x) exp [-ia>t + 2iri. k* r 0 ^)] (11) 

I 
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where r° (t) is as in eq. (2), k is a three-dimensional vector to be subsequently 
determined from boundary conditions, and v a (x) is independent of we find 

= 22 ( 2 ) < i2 > 

k'i 3 

where 

D ^(xt') = ^(^,) exp {2771k. [r°(^)-r°(^]} (13) 

^ M x ra x' 

If one can at this point make the quantities D defined in eq. (13) independent 
of i, then a considerable reduction of the problem is achieved. To this end specific 
assumptions about the second order force constants have to be made: either the 
lattice is infinite and its periodicity insures invariance under rigid body trans- 
lations produced by a lattice vector r ° (-t) , or it is finite with periodic boundary 
conditions. Both cases lead to the relations 


3-vC-a 


(14) 


i.e., the atomic force constants depend only on the difference l - V . 

This being the case the superscript l can be suppressed in eqs. (12) and (13). 
Thus eq. (12) can be rewritten 


a> 2 v a (x) 


£ D - (»•) 

X /3 


Vg(x') 


(15) 


with 


(xx') = £ ** (xx>) exp [ - 2wi k • r0<4)1 (16) 

O' 
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Therefore there are only 3n equations instead of 3Nn or an infinite number of 
them. The set of equations (15) will have a solution if the determinant of the co- 
efficients will vanish 


D 


a/3 



S a/3 8 XK' 


= 0 


(17) 


It is easily seen from the definition (16) that the matrix of the coefficients D 
is Hermitian. The 3n solutions cJ? (k) are then real and the stability of the lattice 
requires that these be positive. 

To a certain extent the reduction effected above is illusory since the vectors 
k which satisfy the boundary conditions for a finite lattice have still to be found. 
This brings us back to the original problem. For an infinite lattice the 3r functions 
cJ (k) can be regarded as the branches of a multivalued function oj 2 (k). 

The plane-wave method is indispensable in all those cases for which the 
normal mode frequencies cannot be explicitly found. Moreover it is the starting 
point of the quantum-mechanical treatment of solids, scattering of waves and 
particles, etc. 

2. The Normal Mode Method 

In this method the following substitution 


u 


a 




e -i*n 


(18) 


is used in eq. (10) to yield 


- L vC"') 





(19) 
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The condition that the set of equations (19) have a nontrivial solution is 


| 4 > 


a/3 


ill) 


KK a/3 


= 0 


( 20 ) 


The roots of this equation are the normal mode frequencies of the crystal. The 
3Nn x 3Nn symmetric matrix in eq. (20) is known as the dynamical matrix of the 
system, usually denoted by A (co 2 ) or simply A. The boundary conditions are 
automatically introduced through the specific form of the dynamical matrix. 
Generally, if the roots of eq. (20) cannot be found analytically, the numerical 
solution is more complicated that that of the reduced eq. (17) and this approach 
is not useful. But if one considers short-range interactions (i.e., interactions 
between distant neighbors can be neglected) and certain simple boundary condi- 
tions, there exist several models which possess explicit analytic solutions. 


It is this latter aspect of the normal mode method which will be utilized 
exclusively in the sequel. We wish to emphasize that the application of the plane 
wave method to these soluble models, though possible, is much less convenient 
than that of the present method. 


Before we proceed to the actual models, a further question has to be answered: 
how to connect the atomic force constants with the geometrical and physical struc- 
ture of a given lattice. It readily appears that knowledge of the invariance proper- 
ties of $ is not sufficient for this purpose. If we restrict the discussion to two- 
body interactions, a more physical approach is through the forces acting between 
pairs of particles. This approach has been used by Bora [ 10 ] and subsequently 
amplified by de Launay [ 12 ] , whose treatment we follow. 

Let 3 1 A be the force exerted on the particle located at (-6 ,k) by the 
particle located at {t ‘ , x') . Then by definition 
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( 21 ) 




where $ is the 3x3 matrix with elements given by eq. (9) , and we see that $ 
is also a force density tensor possessing only 6 independent components. 

The procedure of de Launay is essentially to write for <$ a linear combination 
of dyadic s 


<D 

I 



= a € £ + a ' § § + a " i n 


where 


and 


a, a* , a" = a, a', a" 



c 


5, n = e. n 



One of the three unit vectors, say e , is chosen as follows 



( 22 ) 


(23) 


while the remaining two are taken to form with e a right-handed orthogonal 
triad. There are not enough conditions to make the choice of these two unique 
and therefore a certain degree of arbitrariness remains. Despite this fact, it 
is not difficult to choose these vectors for actual models, as will be described 
below. 
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The coefficients a and a' , a" are known as the central force and non-central 
force constants, respectively. These depend on the type of the two-body inter- 
action and also on the separation distance — the constants will vary for different 
neighbors even when the same type of interaction is involved. 

In the following we present the equations of motion and the dynamical matrices 
for several models of finite lattices in one, two, and three dimensions. 


1.1 One-Dimensional Lattices 

In this section we shall discuss finite linear chains and classify these according 
to their symmetries and the boundary conditions imposed. 


The model assumed is that of N particles M x , M , . . . , 11^ arranged on a 
line and constrained to vibrate longitudinally, as in Fig. 1. The equilibrium 
positions x° obey the conditions 


x !+i ~ x ° = a j ; i = 1 N - 1 


(24) 


M, 




% 


Figure 1 


where the ads denote the spacing constants of the linear lattice. If none of the 
a. ’s coincide there will be N(N - l)/2 independent force constants. This case 
is of no interest here since no periodicity is involved. We consider below the 
case where the two-body interactions are independent of the particle masses and 
assume that the particles, at equilibrium, are equi distantly spaced. Then the 
number of independent force constants reduces to N - 1 . According to our model 
the tensor 4> of eq. (22) reduces in this case to one constant a (•£, V) since e e 
becomes 1. Moreover the properties of the two-body interaction assumed lead to 
force constants a with the property 
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a a, l') = a(K -4'b 3 V j = 1 N - 1 


(25) 


In order to write down the dynamical matrix the boundary conditions have to 
be explicitly introduced. This is done below. 

1.1.1 Linear Chains with Free Ends 

The forces acting on the two end particles and M 2 , according to eqs. (21) 
and (25) are respectively: 




N-l 


J=1 


> 


N-l 

F N “ a N-j ^ u j ” U N^ 

j=l 


J 


Then the dynamical matrix is 


' A 1 - q ? , 


- a i> • 


’ ” a N-l 


V- 2 ) = 


- a, 


A 2 “ M 2 U>2 ’ N ^ 


V 


VrVi^- - a i 

\ 

. N - a, , A 1 - M kT a? 

1 1 N 


(26) 


(27) 


in which 



j-i 



r=l 


j = 2, 


N-l 


k >- H a ' 

r=l 


~\ 

N - 1 


> ( 28 ) 
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The eigenvalues and eigenvectors of the matrix (27) are in general unknown. 
But if the chain is monatomic M. = M, j = 1, . . . , N, one eigenvalue is zero and 
its associated eigenvector is 



\'l 


(29) 


This of course reflects the invariance of the lattice against rigid translations and 
as a result in matrix theory it is a particular case of the somewhat more general 
theorem: If the sums of the elements in each row (column) of a matrix coincide, 
the matrix has one eigenvalue equal to this sum and the associated eigenvector 
is that of eq. (29). 


A monatomic chain for which all interactions can be neglected except those 
between nearest neighbors yields a matrix A of the form 



with 

a = 2a.j - Mw 2 ; b = - (31) 


This matrix appears in eq. (63) of Appendix E and is discussed there. Use 
of eq. (31) above gives then the frequencies 


15 





4a, 

a ? - — - sin 2 — ; r = 0 , 1 , • 

r M 2N 


. , N - 1 


(32) 


The eigenvectors are presented in the same Appendix. 

A monatomic chain with two distinct spacing constants a 1 and a 2 regularly 
alternating is characterized, for nearest neighbor interactions, by two different 
force constants /3, y and leads to the matrix 


= 


/ 


a + c b 
b 


a c 

\ \ 

„ \ \ 


. \ 


\ \ ' 

0 \ a N b 


\ 


b a + b 




(33) 


in which 


a = (3 + y - M o?\ b = — /3; c = - y 

and b , c alternate regularly along the minor diagonals in eq. (33). 


(34) 


This type of matrix is treated in the section following eq. (79') of Appendix 
E . It is shown there that if N = 2n then the eigenvalues and eigenvectors can 
be found explicitly. Using eq. (34) we obtain the frequencies 




77k 


= 


+ y ± ] 

1 10 2 + 2f3y cos + y 2 

1 

f n 


M 


r = 1 , • 


n 


- 1 




= 0 


2 M 

w 2 n - M 


( 35 ) 
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The frequencies to 2 and co 2 n coincide, for /3, y a x , with the two frequencies 
of eq. (32) obtained by the substitution of r = 0 and r = N/2, respectively. We 
note also that the frequency is independent of the number of particles. As 
such it is called a "surface" frequency, a term borrowed from the theory of an 
elastic continuum with free boundaries. At present no analytic results are 
available for the odd case N = 2n + 1. 


The diatomic chain with nearest neighbor interactions and masses alter- 
nating regularly along the chain has the dynamical matrix 


in which 


/ 


u + b b 


V" 2 ' = 


\ 0 

b v N 
\ 

\ 

\ 




\ 

\ 

u \ 

\ \ X b 

0 \ \ 

\ \ 

b * 


u = 2a 1 - Mj oo 2 ; v = 2a x - M 2 co 2 ; b = - 


(36) 


(37) 


and the last term on the main diagonal is either u + b (for N odd) or v + b (for 
N even). As the odd case does not possess at present analytic solutions, we 
consider only N = 2n. Then the eigenvalues and eigenvectors are given in eqs. 
(99), (101) and (102), (103) of Appendix E and we can write the frequencies 


w k = a i 


M x + M 2 ± j/m x -M 2 ) 2 + 4M x M 2 cos 2 ^ 


M i M 2 


k — !»•••> n — 1 


=» 




a t (Mj + M 2 ) 


M x M 2 




> 


( 38 ) 
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Again co 2 is the "surface" frequency and as in the previous case there are 

° 2n 

two frequency branches or bands corresponding to the ± signs in eq. (38). 

For odd N , N = 2n + 1, the only information available is the dispersion 
relation 

u sin(2n + 2) 0 + 2 /uv sin (2n + I)6 l + vsin2n6 l _ q 

sin 2 6 

in which /uv = 2b cos 6. This is of course the characteristic equation of the 
matrix in eq. (36) and can be put in the form 


{/u sin (n + 1) 6 + /y sin (n<9)} {/u cos (n + 1) 6 + /v cos n 6 } = 0 (40) 


which exhibits the existence of two frequency branches also in this case. 

1.1.2 Linear Chain with Fixed Ends 

Here again we assume mass-independent interactions and equidistant spacing 
of the equilibrium sites (unless specified otherwise). The two end particles are 
supposed to interact with rigid walls via the same force constants employed through- 
out the chain. Though more realistic, the inclusion of different force constants 
for the end particles will render the problem tractable only be perturbation 
methods. We restrict therefore the discussion to the former situation. 


The dynamical matrix for the polyatomic case is 


- Mj co 2 , - a x 


= 


- CL 


A 2 -M 2 oK 


■Vi 


X 


\ 


\ 


\ 


\ 


\ 


\ 


\ 


\ 


\ 


\ 


\ 


(41) 


\ 


•" Vi’ 


\ 

• ~ a v 


\ 


X N - ai 
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in which 


22 a > + a i ; 


r = 1, N 


£ “i + 22 V 


r = 2, . . . , N - 1 


As for a chain with free ends, the polyatomic case cannot be treated analyti- 
cally. The cases for which closed analytic expressions of the eigenfrequencies 
are available will be listed below. 


The monatomic chain with nearest neighbor interactions only, has the dynami- 
cal matrix 


& N (cJ) = 


a b 
\ \ 

\\\ 0 

\ \ \ 

\ \ \ 

\ \ b 
0 \ \ 

b a 


= 2a 1 - Moj 2 ; b = - a 


Using eq. (44) and the result (10) of Appendix E, we obtain the frequencies 


cci 2 = — - sin 2 — — , r = 1 , • • • , N 

r M 2(N + 1) 


The eigenvectors are exhibited in eq. (12) of the Appendix. 
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The inclusion of interactions among more distant neighbors in the dynamical 
matrix makes the problem of calculating the eigenvalues analytically extremely 
difficult, and all efforts in this direction have been so far unsuccessful. 


The diatomic chain, with particles M x and M 2 alternating regularly along 
the chain, can be regarded for N = 2n as a lattice with a basis the unit cell of 
which contains two particles of different species. Because of the simplicity of 
the situation, one can disregard this feature and proceed directly. 


For equidistant spacing and mass-independent interactions, the dynamical 
matrix is 





\ \ \ 
\ \ 


\ 


b 



(45) 


in which 


u = 2a x - Mj a ? ; v = 2a.j - M 2 a? 


(46) 


and the last element on the main diagonal is either u (for N odd) or v (for N even). 
The eqs. (36) and (42) of Appendix E show, on using eq. (46) above, that the 
frequencies for these two cases are given by: 


N = 2n 


k Mj M 2 


Mj + M 2 ± 


JwT- 


M 2 ) 2 + 4Mj M 2 cos 2 


7rk 

2n + 1 


k = 1, . . . ,n 


( 47 ) 
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N = 2n + 1 


Cl)} - 


k MjM 2 


Mj + M 2 ± yeM^ -M 2 ) 2 +4M J 


M cos 2 -■ ^ 

2 2n + 2 


] 


k - 1 , . . . , n 


> ( 48 ) 


5 a. 


a;. 


2n + l 


M, 


The dynamical matrix for a diatomic chain with two different spacing constants 
similarly alternating is 





(49) 


in which 


u = + 7 - M l cJ 2 ; v = / 6+y-M 2 oj 2 ; b = - fi; c - - y (50) 


and (3, y are the two different force constants between nearest neighbors. Eqs. 
(56) and (59) of Appendix E show the frequencies to be 


co} = — L— 

k 2M X M 2 L 


(y8 + y) (M x + M 2 ) ± j/oS + y) 2 (M x + M 2 ) 2 - 16 M 2 sin 2 ** 


2n + 2 J 


k — 1 y ■ • • y n 


> 


Cl), 


2n+l 


- @_Ll 

M. 


( 51 ) 
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Here again we observe the separation of the frequencies into two distinct 
branches. If this type of diatomic chain consists of an even number of particles. 

N = 2n , exact results are no longer available. 

1.1.3 Linear Chain with Periodic Boundary Conditions 

The chain is assumed here to form a closed loop. This means that there are 
no end particles and the interactions proceed along the shorter separation dis- 
tances on the loop. Use of the forces defined in eq. (21) with mass independence 
and the assumption of equidistant spacing, yield for the polyatomic case the dynami- 
cal matrix 


f di ag (A- MjM 2 , A-M 2 cJ, . . . , A-M^) -(a^ . . .a n . . -aj) cyc 


N = 2n 




(52) 


di ag (A-MjW 2 , A - M 2 w 2 A-M^ 2 ) - (a } a 2 . . .a n a n a n _j . • •a 1 ) cyc 


^ N = 2n + 1 

in which diag (**...*) denotes a diagonal matrix with the elements specified in 
parantheses, and (**. . . *) cyc denotes a matrix the first row of which consists 
of the specified elements *, the other rows being obtained from the first by 
(counter-clockwise) cyclic permutations. Also 


f 2 


n-1 

E 


a i + v 


N = 2n 


A = "S 


2 


n 



j=l 


cl . ; 

j 


N = 2n + 1 


(53) 
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We note that the number of independent force constants is here reduced to 


either N/2 or (N - l)/2. 


In contrast to previous cases, the frequencies and eigenvectors of the monatomic 
chain are known explicitly for the general case and can be given from eqs. (153) 
and (154) of Appendix E as follows: 


r 


n - 1 


1 r - 

M La J 


sin" 


77k j 


a 


2 n M 


+ _2_[l_(-l) k ] 


i = i 


k = 0, 1 , • • • , 2n - 1 ; N = 2n 




(54) 


n 

4 V 1 

— / a. sin' 

M Z__i 1 


nk j 
2n + 1 


j = l 


k = 0,l, •••, 2n; N = 2n + 1 




It is seen that one frequency (k =0) vanishes just as for a chain with free ends. 


The frequencies and eigenvectors of a diatomic chain with an even number of 
particles can also be found. Using the eqs. (135) and (140) of Appendix E, we 
find the frequencies 


00t 


(Mi +M 2 ) 2 
SM^ 


k k + n 


(^-^ + n> 2 +4 


Mj - M 2 
M x + M 2 


k k + n 


k = 0, * • • , n - 1 


(55) 


The corresponding result in the literature is given only for nearest neighbors interactions and is 
not as transparent as the above. 
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in which 



2 

+ m 2 



n - 1 



. 2 7T r 1 r , 

a. sin' + a LI 

> 2 n n 



* 


r = 0, 2 n - 1 (56) 

It is of interest to note that when the odd-indexed a's vanish, the radical in 
eq. (55) disappears and a decoupling of the motion takes place, i.e., the resulting 
frequencies are those of two separate loops with particles M x and M 2 respectively. 

Here also the diatomic chain with an odd number of particles cannot be solved 
analytically. 

Other situations can be also treated exactly. For instance if, everything being 
as before, we introduce two distinct spacing constants, as in Fig. 2, the dynamical 
matrix for a chain with only nearest neighbor interactions, will be 


A 


2n 



(57) 


in which 


u=/3 + y- M 1 o) 2 ; v=/3+y-M 2 ai 2 ; b = -/3; c = -y 


These are just the frequencies for a monatomic chain with 


particle mass(Mj +M 2 )/2. 


( 58 ) 
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Figure 2 Figure 3 


On using (58) above and eq. (161) of Appendix E we find the frequencies 



1.2 Two-Dimensional Lattices 

The plane lattices to be considered in this section will be rectangular. The 
model assumed is shown in Fig. 3: the unit cell is based on the vectors a x = 
a x i; a 2 = a 2 j , where i, j are the unit vectors in the x- and y- directions, 
respectively. We assume the lattice to consist of N 2 horizontal rows equidistantly 
separated, with N x particles equidistantly spaced in each row. 

The discussion will be restricted here to at most eight neighbors (for parti- 
cles not on the boundary). For instance, the particle located at the site a x + a 2 
is assumed to interact via two-body forces only with the particles at the sites 
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0 (this is labeled 1 in Fig. 3), a x , 2 a x , a 2 , 2 a 2 , 2 a x + a 2 , a x + 2 a 2 and 2(a x + a 2 ) 
For these sites to represent indeed the first, second and third neighbors, we have 
to assume 


a i- a 2 > \ F 


+ a 


2 ’ 


The particles will be constrained to vibrate only in the plane of the lattice. 
The force density tensor $ will then be that of eq. (22), with a" suppressed. 
The eight unit vectors t needed in eq. (22), and their counterparts §, are given 
here by 


~\ 


£ 1 = " e 2 = 1 


E 3 = " £ 4 


§1 = -§ 2 = j 

*3 = -*4=- i 


" £ 6 = 


E 7 = " £ 8 = 


+ a 2 


fi\ 


= 


" a i + “2 


> 


(60) 


2 +a2 




2 + a 2 


" a i + a 2 


F\ 


2 , o 2 

+ a 2 




8 1 + a 2 


/ 3 1 +a 2 

Though immaterial in the two-dimensional case, we have used here the convention 
that angles between unit vectors are measured counter-clockwise from t to § . 

There are three central and three noncentral force constants associated with 


the unit vectors of (60), in the following manner 

£ 2 } ~ a i ; < £ 3’ £ 4 } ~ a 2 : {e 5’ V V £ 8 } ~ a 3 

and similar relations for 5 with a' . 


(61) 


We consider below different boundary conditions. 


1.2.1 Plane Lattice with Free Boundaries 

The equations of motion and the dynamical matrix are found in Appendix A. 
Here we quote the results. 
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A monatomic lattice has the dynamical matrix 



where G denotes the transpose of G, and A t , A E are defined in eqs. (6) 

and (7) of Appendix A. 

Two eigenfrequencies can be seen from the equations of motion to vanish, 
corresponding of course to a rigid translation of the lattice. The two independent 
eigenvectors for this case are obtained when putting in turn 
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( 64 ) 






for all (^,m ), withU defined in eq. (4) of Appendix A. 

The remaining eigenfrequencies cannot be found analytically in the general 
case. But if next-nearest neighbor interactions are neglected, then the following 
simplifications take place 

D = E = 0 ; A X =A + B; A 2 = A + C ; A 3 =A 4 =A + B+ C (65) 

and 

G = G; F' = F" = F + G (66) 


Note that for this case the motions of the atoms in the x- and y-direction are 
independent. 


Using (65), the dynamical matrix simplifies to 


A ('a; 2 ') 

U 2Nj N 2 k > 



(67) 


As F and G are now commuting matrices, we can use eq. (203) of Appendix E 
and the definitions of the 2x2 matrices A, B , C to write the frequencies: 




k‘>i 

2 _ 4 a 

2 1 yQ 4 P 
sin 2 k & + — — 

sin 2 j cp 





L kj J 

~ nr 

M 

O 

it 

• *, Nj 

- 1 






y 



(68) 


2 _ 4a' 

sin 2 k 8 + — — 

sin 2 j cp 

j =o,- 

• *, n 2 

- 1 


L k i - 

” IT 

M 






J 


in which 8 = 7T / 2 N x ; cp = 77-/2N 2 . 
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It is seen therefore that a rectangular plane lattice, in this model, has 2N X N 2 - 1 
distinct normal frequencies. This property persists even for the case N x = N 2 
or when the unit cell is a square (for which a. = /3 and a' = ft ). On the other 
hand, if N x = N 2 and the unit cell is square, the number of distinct frequencies 
is reduced to N x N 2 , each being doubly degenerate. 


It is not difficult now to write down the dynamical matrix for the diatomic 
lattice. In addition to the two masses one has to take into account the three types 
of interaction possible: 


Mj— Mj ; Mj— ■ - m 2 ; M 2 -~-M 2 

The case Mj*~M 2 corresponds here to nearest neighbors and the other two to 
next nearest neighbors. Since we cannot find the eigenvalues for more distant 
neighbors even in the monatomic case, we shall restrict the discussion to nearest 
neighbors interactions. The dynamical matrix for this case is given then essentially 
by eq. (67), with the following modifications: 1. The 2 x 2- A matrices alternate 
regularly along the main diagonal of F < > : 




A (l) + B B 


B A< 2 > 
\ \ 


F^ 1 ) = 


\ 


\ 


\ 


\ 


\ 0 
\ 

\ 

\ 


0 \ 


B 


\ \ \ / 

\ B A< > + B / 


(69) 


in which 


A<»> 


A (2) 


/ 2 (a + ft 

\ o 

/2 (a +/3')-M 2 u> 2 

\ 0 


0 

2 (a' + ft) - MjO) 2 
0 


2 (a' + ft) - M 2 u> 2 



> 


J 
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and 


p(2) 


’ A (2 ' + B B 


\ 


\ 


B A<‘>\ 0 

\ \ \ 

X \ N 

0 X \ X B 
\ \ 

B V > +B / 


(71) 


The dynamical matrix has the form 


A (oo 2 ) = 


/ F<‘>+G 

C x 0 

\ \ \ 

\ \ \ 

n \ N N 

0 \ 


\ 


\ ^ N<5 

X G F< > + G / 


(72) 


where the matrices F (1) , F (2) alternate regularly along the main diagonal. 

Though certain reductions of the dynamical matrix in eq. (72) are possible, 

analytic expressions for the frequencies cannot be found. 

• • • • 

• • • • 

• • • • 

• • • • 

o o o o 

\ 

• • • • 

Figure 4. •-.Mjjo-Mj 


On the otherhand if we consider the lattice shown in Fig. 4, which consists of 
alternating rows of two types of particles ,M 2 , and assume N 2 = 2n 2 , the dynami- 
cal matrix will have the form 
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I F< 1 >+G G 


A (oj 2 ) = 


A 


0 


G F (2) \ 

\ \ \ 

\ x \ 

\ \ G 

0 x \ 

N G F< 2 >+G/ 


(73) 


in which F ( 1 > and F (2) alternate along the main diagonal of A and 


'A<j>+B B 

\ 

B A<>^ 


F (i ) = 


\ 


x \° 

X \ 


X 

0 \ 


B 


B 


A< j > + B / 


j = 1,2 


(74) 


The matrices A^\ A (2) ,B and G are as defined previously. As F (1) , 
F (1) and G are commuting in pairs, it is shown in Appendix E that the eigen- 
values of A can be found. Using then eq. (209) of that Appendix, we obtain two 
groups of frequencies, each group consisting of four branches: 


Group 1 


o* - 
kj 


= Ti~7r” i( M i + m 2 ) (£' + 2 a sin2 j 

Mi M 2 I 


1 


y (M x -M 2 ) 2 (/?' + 2 a sin 2 j <9) 2 + 4M x M 2 /3' 2 cos 2 kcp J 




= i( M i +m 2> (^ + 2a ' sin2 J 0 ) 

M ^ M 2 ^ 

]/ (Mj - M 2 ) 2 (J3 + 2 a' sin 2 j 0) 2 + 4M x M 2 /3 2 cos 2 kcp j" 


(75) 
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in which 


j-0, --, ,N 1 -l;k = l,***,n -l; 6 = — — ; c p = 


2N, 


77 

2n„ 


(76) 


Group 2 


0)? = 


MjM 2 


JcMj +M 2 )^y + 2 a sin 2 ± j/(M^ -M 2 ) 2 ^-+ 2 a sin 2 j^ 2 + MjMj/3' 2 


oj? = • 


MjM 2 


i(M J +M 2 )^-y-+ 2a' sin 2 j ± |/(Mj -M 2 ) 2 ^-+ 2a'sin 2 j + MjM 2 /3 : 


(77) 


with 0 as in eq. (76) and j = 0, . . . , Nj - 1. 

The frequencies in the last group are those associated with the surface modes. 
Two of these vanish (for j = 0), corresponding to rigid translations of the lattice. 


1.2.2 Plane Lattice with Rigid Boundaries 

Here we assume the marginal particles to interact with the rigid boundaries 
via the same two-body forces operating inside the lattice. Therefore using the 
same notation as before, the only equations of motion that change are those for 
the marginal particles. Proceeding as in§ 1.2 of Appendix A we find the dynamical 
matrix for a monatomic lattice to be of the form: 


F G 

\ \ 
r.\ \ 


\ 


A (a) 2 ) 


2 \ - 


\ 


\ \ 

\W 

\ ' \ 

\ \ 

G F 


/ n,xn, 


(78) 
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in which 


F = 


B 

\ \ 
B x x 


\ 


\ 


\ \ 

\ \ 


\ 


\ \ 

\ x x 

\V B 


; G = 


\ B Na /n i xm i 




C D 


\ \ 

E \ \ 

\ \ \ 

\ \ \ 

\ \ \ 


\ \ x 
\ \ 

\ \ 

E 


D 

C 



X N, 


(79) 


and A, B , . . . , E are as defined in eq. ( 6 ) of Appendix A. 


Again we note that F, G ,G do not commute. This stems from the fact that 
D and E do not commute with one another or with A, B. The coupling between 
the motions in the x- and y-directions is given by D and E , more specifically by 
the matrix element 8 . It is clear therefore, from the expression defining 8 , that 
if ~ < 2.3 , then S ~ 0 and all the 2x2 matrices involved become diagonal, with 

f-y 0 \ 

D = E = I ; G = G (80) 


Assuming the condition 8 ~ 0 to be satisfied, the motions in the x- and y-directions 
become independent and the frequencies will be obtained by equating to zero the 
A kj given in eq. (175) of Appendix E, namely 


A 


A + 2 C co^ 77 ^ 

+ 2 

B + 2 D cos 77 ^ 

n 2 + ! _ 


L n 2 + 1 J 


TTk n 
cos = 0 

Nj +1 




k = 1 , • • • , Nj ; j = 1 , • • • , N 2 
Finally two groups of frequencies are obtained 


(81) 


J 


kj M 


o 2a /„ 77 k 

- 1 1 - cos 


Nj +1 


2 P 

M 


1 - cos 


77 j 


/ 1-, 


7rk 7T i 

, +— ^|l-cos cos 

N 2 + 1 J M ^ Nj + 1 N 2 + 1 
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( 82 ) 



and 


CO 


2 _ 
kj 



1 - COS 



2/3 

M 


1 - cos 



4 y 

"IT 


[1 - cos 


7rk 

nTTi 


. cos 


7T ) } 

N 2+ iy 


(83) 


with k and j as in eq. (81) . 

We remark that the boundary conditions under consideration do not lead to 
degeneracies in the spectrum, even for a lattice with a square unit cell (in which 
case a = /3; a' = ft ) . The degeneracy that exists for a square lattice (N x = N 2 ) 
with square unit cell, is removed by the introduction of the next-nearest neighbor 
interactions, represented here by y . 

For the diatomic lattice with regular alternation of M x and M 2 , the changes 

r 

in the dynamical matrix are likewise obvious and we can write: 


A (<^> 2 ) = 


F ' \ 

« F, \ 


\ 


\ 


\ 


\ 


X \\ 

\ \ \ 

\ \ \ p 

G F 


/ 


(84) 


in which F x and F 2 alternate regularly along the main diagonal and 


B 

B A( 2> >\ ' 0 
\ \ \ 

\ \ \ 

N \ \ 

x \ B 
\ 


\ 


0 


\ 


\ 


\ 


B 


\ 


A< > 


/ 


; F 0 


/A< 2 > 

B 

\ 

\ 

0 

B v 

Ad) \ 

\ 

\ \ 


\ 

\ \ 


\ 

\ 

\ 


\ v \ 

\ 

l 0 

\ 

\ 

\ 

\ 

\B 

\ 


B 

X A 


\ 


<>/ 
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The matrices A ( 1 \ A (2) which alternate regularly in Fj , F 2 are as defined 
in eq. (70), D and E as in eq. (80), G as in eq. (79), andB ,C as in eq. (6) of 
Appendix A. 

The eigenvalues of b(a?) can be found explicitly, for all parities of ,N 2 , 
as shown in the treatment following eq. (177) of Appendix E. Here we quote re- 
sults only for the case N 2 = 2n 2 , and use eq. (195) in Appendix E (with A r re- 
placed by A ( r ) , r = 1, 2) to write down the frequencies. These appear in four 
distinct branches: 


CO 


?k ~ TTT7 V M i 4M 2 ) C a + /3' + 2y (1 - cos j 0 cos kcp] ± 

M ^ 2 ^ 

± |/(M 1 -M 2 ) 2 [a + /? + 2y (1 - cosj0coskcp] 2 +4M 1 M 2 [a.cosj/9 + /3'cosk9] 2 j' 


00 


• k j( M i + M 2 ) +fi +2 y (1 -cos j 6 cos ktp) ] ± 

± y^(M 1 -M 2 ) 2 [a'+/S + 2y (1 - cos j 0COS kcp)] 2 +4M X M 2 [a'cos j 0 H-yScos kcp] 2 I" 


with 


k = 


1 , • • * , n ; j = 1 , • • • , Nj ; 6 - — — ; 

* Nj + 1 


<P = 


77 

2n 2 + 1 


It is of interest to compare the frequency spectrum of this model with the one 
obtained when a rearrangement of the particles has taken place. The rearrangement 
envisaged is the one already considered in Fig. 4 of the previous section, namely 
alternating rows of two types of particles, each row consisting of identical particles. 


The dynamical matrix is (where for simplicity we assume N 2 = 2n 2 ); 
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and A (1) ,A (2) ,B ,G are as in eq. (84) . 

The matrices F x , F 2 and G commute in pairs, and so do all the 2X2 matrices 
involved. This being the case, eq. (183) of Appendix E shows the frequencies to be 
the roots of the matric equations 



j~A (1) + A (2) + 2B cog j q + 2 [C + 2 D cos j 9 ] cos k cp| 
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Finally, the frequencies we obtain, appear in four distinct branches: 


"kj |( m i + m 2 ) t a (l- cos i e ) +/?+2y] ± 


± |/(M X -M 2 ) 2 [a(l - cos j 0) + /S' + 2y] 2 + 4 MjM 2 [/3* + 2y cos j (9] 2 cos 2 kcp 


"kj |( M i +M 2) [“'C 1 " cos j 6 ) + Z 3 + 2 7 J ± 


± ^/(Mj - M 2 ) 2 [a' (1 - cos j(9) +/3 + 2 y ] 2 +4M X M 2 Q3+2')'cosj£] 2 cos 2 kcp 

(92) 

with k , j , 6 and cp as in eq. (91) . 

Similar expressions can be easily found also in case N, is odd. 


1.2.3 Plane Lattice with Periodic Boundary Conditions 

Here particles on opposite boundaries interact with one another and therefore 
we can envisage the lattice particles as being located on the surface of a torus. 
This means that there are no marginal particles. 


One can write down the dynamical matrix taking account of all possible inter- 
actions, since in fact the eigenfrequencies can be explicitly exhibited for this 
general case. For the sake of simplicity and comparison with preceding results, 
we consider below only nearest and next-nearest neighbor interactions. 


It is not difficult to show that the dynamical matrix for the monatomic lattice is 


A (co 2 ) = 


' F G G 

~ \ 

G \ \ \ 0 
\ \ \ 

\ \ \ 

• \ w 


\ 


\ 

-V, \ 

G F 



(93) 
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cv 


where 


F = 


A B B 

B \\ 

\\\° 

\ \ \ 
° \ \ B 

\ \ 

B BA 


M. 


; G - 


\B BA/, 


E ' 

\ N 


X \\° 


' ' \ 
o \ \ ' 


; G - 


NjXNj 


\ \\ D 

\d e c/ BiXNi 


C E D 

■w . 

\ \ \ 

\ v \ 

\ ° \ \ E , 

\e d \ / 


(94) 


The matrices A, . . . , E are as defined in eq. (6) of Appendix A. 

It is seen that A (&j 2 ) is a two-dimensional circulant matrix and that F, G and 
are circulant in their elements. Then we can use eq. (216) of Appendix E to 
write down the following matric equation the roots of which are the frequencies 
of the system: 

0 = A + 2 B cos j 0 + 2 C cos k <p + 2 D cos (k <p + j 8) + 

+ 2E cos (kcp - j 8) = A k j ; k = 0, * * * , N 2 - 1 ; j = 0, * * * , Nj - 1 (95) 

in which 8 - 2 77/Nj ; cp = 2 tt/N 2 . 

The explicit form of the 2x2 matrix A kj is 

( 4 a sin 2 (j 8/2 ) +4/3' sin 2 (k cp/2) + 4 / (1 - cos j 8 cos kcp) - Mo?', 48 sin j 8 sin kcp'' 

4 8 sin j 8 sin kcp ; 4 a' sin 2 ( j 8/2) + 4 0 sin 2 (kcp/ 2) + 4/ (1 - cos j 8 cos kcp) -Mw 2 

(96) 

Finally the frequencies are given by 

a > 2 =— -|(a + a* ) sin 2 ( j 0/2) + (/3 + /3* ) sin 2 (kcp/ 2) + (7 + /' ) (1 - cos }8 coskcp) ± 
kJ M 

± /[(a-a')sin 2 (j(9/2)+ (/S' -( 3 ) sin 2 (k<p/2) + (7 - 7* ) (1 - cos j 8 cosk<p) ] 2 + (2 8 sin j <9sinkcp) 2 j* 
j = 0, • • • , Nj - 1 ; k = 0 , • • • , N 2 -l, 8 = 2 77/Nj ; <p = 2 tt/N 2 (97) 
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We note that here we obtain two branches, a result to be expected from the in- 
clusion of next-nearest neighbor interactions. If S is negligible, then the x- and 
y-motions are independent and the frequency spectrum reduces to one branch 
only. We observe that for j = k = 0 we obtain a vanishing, doubly degenerate 
frequency. 


When we consider the diatomic lattice the only solvable case is when both N x 
andN 2 are even. Then, if N x = 2n l , N 2 = 2n 2 ,the dynamical matrix takes the form: 


F x G n 


AO 2 ) = 


G F N 0 
\ r 2 \ 

\ ' x 

\ \ \ 

° \ F, 'g 


\ 


(98) 


n 2 xn 2 


in which G and G are as for the monatomic case, while 


F, = 


/V 1 ) B B 

\ 0 

B, A<*>\ 

\ ' V 

0 \ v \ 

U \ A* 1 ) B 


B 


\ 


B A (2) i 


; f, = 


/A< 2 > B v 

B A* 1 ’ X 


0 

\ 


B 


\ \ \ 

0 \ V 2 ) X B 

B B A^> 


(99) 


and A (1 0 A< 2) as in eq. (84). 


The matrices F x , F 2 and G do not commute. Therefore the appropriate 
treatment is that preceding eq. (241) of Appendix E where it is shown that 
AO 2 ) is reducible to two types of block matrices on the main diagonal. These 
matrices denoted K (0) and K (1) in eq. (241), become here 4x4 matrices. 
These matrices cannot be further simplified and the resulting equations are of 
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the fourth degree in co 2 . Hence, for this model, we obtain eight frequency branches. 
Analytic expressions for the frequencies can be found only for 4 (Nj + N 2 - 2) 
of them — corresponding to the cases for which the coefficients of the coupling 
parameter S vanish, making all the 2x2 matrices diagonal. Below we exhibit 
the 2 NjN 2 frequencies obtained when S is dropped from the start: 




cot 


. - — — •[ (M + M ) [a + /?* + 2y (1 - cos j 6 cos kcp) ] ± 

J Mj M 2 L 

± J (Mj -M 2 ) 2 [a + /3* + 2y(l - cos j ~6 cos kcp)] 2 + 4 MjM 2 [a cos j ~d + /3'cos kcp)] 2 ^ 

( 100 ) 

j = m"1T {^ m i + I’ 0 '* + ^ + 2y ^ ” cos 1 e cos kcp ^ ^ 1 

± ^CMj-Mj) 2 [a' +/3 + 2 y (1 - cos j ~6 cos kcp)] 2 +4 MjM 2 [a' cos j 0 + cos kcp] 2 | 


( 101 ) 




' Mj M 2 { (M ‘ 


+ M ) [a + /?' + 2y (1 + cos j 6 cos kqp) ] ± 


0? 

k 


± y^(Mj -M ) 2 [a + /3' + 2y (1 +cosj ~6 cos kcp)] 2 + 4M X M 2 [a cos j 0-/3' cos kcp] 2 ^ 

(102) 

■ = ttAt {( m i + M o ^ t a ' + P + 2 7 (! + cos i 19 cos k, p)] ± 

1 M 2 L 

+ - M 2 ) 2 [a' +/3 + 2y (1 + cos j ~Q cos kcp)] 2 + 4 MjM 2 [a' cos j ~Q -/3 cos kqp] 2 ^ 

(103) 

In the last four expressions we have 


j = 0, . . . , nj - 1 ; k = 0, . . . , n 2 - 1 ; 9 = ; <p = 77/n 2 (104) 

The 4(Nj + N 2 - 2) frequencies mentioned previously can be obtained from 
eq. (100) - (102) by letting j and k take in turn the value 0. Two frequencies, 
obtained from eqs. (100) and (10i) for k = j = 0, vanish. 
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1.3 Three-Dimensional Lattices 


In this section we consider a rectangular space lattice. The model assumed 
is shown in Fig. 5. The unit cell is based on the vectors = a 1 i, a 2 = a 2 j ; 
a 3 = a 3 k , where i , j , k are unit vectors in the x-, y- and z-directions, respectively. 
We further assume the lattice to consist of N x particles in each row parallel to 
the x-axis, N 2 and N 3 particles in rows parallel to the y- and z-axes, respectively. 


z 



Figure 5 


No constraints are imposed here on the vibrations of the particles and hence 
the force tensor <t> will be used in the form given by eq. (22). The discussion 
will include 26 immediate neighbors, located at the vertices of 8 unit cells. For 
these to represent first, second and third neighbors, certain not too stringent 
conditions have to be imposed on the relative values of a 3 , a 2 and a 3 . 
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As before, the two-body interactions considered are assumed to be 
mas s -independent. 

The unit vectors t , § , n and the force constants associated with the twenty- 
six neighbors are exhibited in section 2 of Appendix A. 


1.3.1 Lattice with Free Boundaries 

Following the procedure described in Appendix A, the dynamical matrix for 
the monatomic lattice could be written down so as to include all the neighbors 
considered in the present model. Since such generality precludes an exact 
solution for the frequencies, the discussion below will be restricted to the case 
of nearest neighbors interactions. The dynamical matrix for this case is: 


in which 


/ 8tT \ 0 \ 


A (v> 2 ) = 


T S N 
\ \ 

\ \ 

\ \ 


S T 


\ 


T S 


+ N 


2 ‘2 


(105) 


/F + G 


s = 




G F 

\ \ 


. \ 


H 

\ 


\ 


\ N 
\ \ 


0 \ F G 

\ 

G F + G 




\ 


\ 0 

\ 

\ 

\ 

\ 

0 \ 


(106) 


3 '’3 


H * 3**3 
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and 



The matrices B r are as in eqs. (12) of Appendix A and A as in eq. (18) of the 
same, with all C r and D r suppressed. All the 3x3 matrices involved are now 
diagonal. 


The frequencies can be found immediately on using eq. (268) of Appendix E: 


r a. sin 2 j 6 + a" sin 2 kcp + a' 3 sin 2 T\p 

2 = — J a' sin 2 i 9 + a. sin 2 kcp + a" sin 2 ri p 

jkr Ml 1 J 2 T 3 T 

^a" sin 2 j 8 + a' sin 2 kcp + a 3 sin 2 r 0 


in which 


j = 0 N _ 1 ; k = 0,...,N-l; r = 0 N, - 1 | 

> (109) 

0 - 7t/2N x ; cp=7r/2N 2 ; </j = 7t/2N 3 J 

The zero frequency a> 000 is triply degenerate, while the rest are distinct. 
These frequencies will remain distinct even for a cubic unit cell, if at least one 
dimension (N. ) of the lattice is different from the other two. For a cubic lattice 
with a cubic unit cell all the frequencies (108) become triply degenerate in this 
approximation. 


The dynamical matrix for the diatomic lattice can be readily written down, 
but as in the corresponding plane lattice the frequencies cannot be found exactly. 
Here we mention only that exact solutions exist for a lattice composed of alternating 
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planes of two particle- species, each plane containing like particles. In the approxi- 
mation of nearest neighbors, this type of lattice possesses twelve distinct frequency 
branches. 


1.3.2 Lattice with Rigid Boundaries 

The dynamical matrix for the monatomic lattice can be written down when 
interactions up to third neighbors are included, but the frequencies of the system 
cannot be found analytically. Therefore we restrict the discussion to the following 
soluble case: 1. Third order neighbors are neglected; 2. The off-diagonal elements 
of the matrices C r defined in Appendix A are dropped. This last requirement is 
equivalent to assuming / 3 . /3 ' ; j = 1, 2, 3. 

Assuming these conditions to be satisfied, the dynamical matrix becomes 


with 


A(a) 2 ) = 


/S T 




-\\ 

\ \ ^ 

\ x \ 

\ \ ' 

\ \ T 

\ \ 

\ \ 

T S 




N 2 XN 2 


( 110 ) 


S = 


F G 

' ' 0 
G X X 


\ 


\ 


\ \ 

\ ^ 


\ 


\ 


\ 


\ X G 
0 A 

G ' F 


/ 


\ 


'H K 

\ \ 0 
K v \ 

\ \ \ 

\ \ \ 

\ \ \ 

\ \ K 


0 \ \ 


« / 


( 111 ) 
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Here the matrices Bj are as in eq. (12) of Appendix A while according to 
requirement (2) above we obtain 


and 


C, =C 2 = -diag FJ 

C 3 = C 4 = - diag C/3 2 /3" /3 2 ) 

c s =c 6 = -di ag (/r 3 e 3 fi 3 ) 




> (H3) 


J 


A = diag {2 [a^ + a " 2 + + 2(^ + /S 2 + /^)] - Ma; 2 ; 2 + a 2 + a" + +/?' + /8 3 )] 

2 [a" + a' + a 3 + 2(y8J + yS 2 + /8 3 )] ■ - } < 114 ) 

On using eq. (244) of Appendix E the frequencies can be written down 
immediately: 
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r 




cij sin 2 j 6 + a" sin 2 kcp + a 3 sin 2 rt// + /3j (1 - cos 2 j 6 cos 2kcp) 

+ /? 2 (1 - cos2ji9 cos2 ri/0 + / 3 “ (i - cos 2kcp cos 2ri/d 
a' sin 2 j # + a 2 sin 2 kcp + a" sin 2 r0 + /3 1 (l-cos 2j£cos2k<p) 

+ p‘ 2 (1 -cos2j<9 cos2r0) + /3 3 (1 -cos2kcpcos r\p) 


(115) 


a" sin 2 j 8 + a' sin 2 kcp + a 3 sin 2 rv/z + ^'C 1 -cos 2jd? cos 2kcp) 

+ /3 2 (l-cos2j£cos 2r< /») +/3 3 (1 - cos2kcp cos 2r</0 


in which 


j =1,..., N x ; k =1,...,N 2 ; r = l N 3 ;0 = 


77 


; cp = 


77 


; 0 


77 


2(N t +l) 2(N 2 + 1) 2(N 3 + 1) 


(116) 


If instead of neglecting the third neighbors interactions entirely we assume 
the force constants to be nearly equal, 


7 ~ y ~ y" 


(117) 


then all the matrices D defined in Appendix A become equal to -7 1 . The fre- 
quencies for this case can again be found and will be as given in eq. (115), to 
each being added the term 

2y(l - eos2j0cos2kcpcos2ri//). (116) 


Whether this term is added or not, we see that this model posses three distinct 
frequency branches, a property that persists for a lattice with cubic unit cell. 
On the other hand if such a lattice is cubic (i.e., = N 2 = N 3 ) the frequencies 

become triply degenerate. 


We consider now a diatomic lattice and shall assume the arrangement to be 
such that no two adjacent particles are alike. Then the dynamical matrix reads: 
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in which 



The matrices Sj , F. , A 4 ( i = 1, 2) alternate regularly along the main 
diagonals of the appropriate matrices. G and T are as in eq. (Ill), while A i is 
as given in eq. (114) with M replaced by M., i = 1, 2. We have chosen , N 2 , 
N 3 , to be all even only to simplify the presentation, but in fact all possible choices 
lead, in this model, to exact solutions. 
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The frequencies of the system can be obtained by equating to zero the A* 

j K r 

in eq. (265) of Appendix E and solving for w 2 . Here we make use from the 
start of eq. (117) to include the third neighbors interactions. For brevity, we 
present the frequencies a, 2 as elements of a 3 x 3 diagonal matrix £2* kr » 

j- («i ♦ P lk , * )/ (Mj -M 2 ) 2 P? kr +4M 1 M a (^ fcr | < 121 > 

in which 

P jkr =B, + B 2 + B 3 + 2Cj (1 - cos j 6 cos kcp)+ 2C 3 (1 - cos ]6 cos r 0) 

+ 2C S (1 - cos k cp cos r i//) - y 1(1 - cos j 6> co s kcp co s r i/<) (122) 

Q jkr = Bj cos j 8 + B 2 cos k cp + B 3 cos r i p . (123) 

All of the matrices appearing above are as in the monatomic case. On the 
other hand, we have here 


j = 1, • • • , 2n t ; k = 1, . . . , n ; r = 1, . . . , 2n 3 ; Q - — ; q, = — IL_ ; 0 = — — - • 

Nj + 1 N 2 + l N 3+ l 

Eq. (122) shows the existence in this model, of six distinct frequency branches. 


1.3.3 Lattice With Periodic Boundaries 

The dynamical matrix for the monatomic lattice to be displayed here contains 
all the interactions assumed in the model. 


A(w 2 ) = 


/ S T T 

( t\\o 

l '°'\\ T 

\ T T S 


\ 

/ 


N. 


(123) 
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in which 


S = 


F G C 

\ \ 

G\ \ 0 

\\\ 

o\\'g 

0 Vr. 


\g 'g X f /„ 


T = 


H K I 

\ \ 

L \ \ 0 

\ \ \ 

\ \ \ 

\ \ \ 

0 \ \ K 

\ \ 

K L 


.L 


F = 


/b 2 c, 


H 


\ \ 

Vx> 

\'\ 

0 \ \ c . 

\ \ 

Cj C 2 B 2 


'a Bj Bj 
B i X \ 0 

\\\ 

l ° '' \ *' / 

\ B i B i A / 

A 


G = 


/ B 3 C 3 C 4 

C 4 \ \ 0 

\ \ \ 

\ \ \ 

° \ \ X 

\ \ 


C 4 B 3 / N, 




D A 


; K 


C 5 D, 

“A \ 0 

\\ \ 

X X \ X °i 
0 \ ' 1 


; L 


/ C 6 B 4 D 2 \ 

f D 2 \\ 0 ' 
\ \ N 
\ \ \ 

N N X I>4 

0 \ \ 4 


D 4 0 2 X C 


a, y, ^ x c 5 /», 

All of the matrices A , Bj , . . . , D 4 are as defined in § 2 of Appendix A. 


y (124) 




The frequencies are obtained by equating to zero the determinant of the 
matrix in eq. (274) of Appendix E, thus leading to a third degree equation for co 2 . 
Explicit expressions can be derived if we assume the validity of the simplified 
conditions exhibited in eqs. (113) and (117), for which all the 3x3 matrices 
become diagonal. 


For this case the frequencies w 2 are the elements of the 3x3 diagonal 
matrix £2 jkr given by 
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Q jkr = --j^- j B i ( 1 “ cos j <9) + B 2 (1 - cos k(jp) + B 3 (1 - cos r 0) + 

+ 2Cj (1 - cos j 0 cos kcp) + 2C 3 (1 - cos ji9cosri//) + 2C 5 (1 - cos kcpcos r< p) 
-4yl (1-cos j 8 cos k cp cos r \p) } (125) 

in which 

j =0,- • Nj-1; k = 0,- N 2 -l; r = 0,- • N 3 -l; 6 = 2 ^ ; cp = 2 tt/N 2 ; <A= 2 tt/N 3 

(126) 

and the matrices Bj are as in eq. (12) of Appendix A, while C as in eq. (113) 
above. 


Three distinct frequency branches appear also here. 


The only soluble case for a diatomic lattice is when N x , N 2 , N 3 are all even. 
This is in contrast to the diatomic lattice with rigid boundaries. The changes 
in the dynamical matrix of eq. (123) are as follows: instead of the S given in 
eq. (124) there are Sj and S 2 alternating on the main diagonal of A(w 2 ) , each 
one with alternating Fj and F 2 on their main diagonals, 


F, G G 
1 \ 

»F 1 \ o 

\ \ \ 

\ \ \ 
o \ \ \ 

\ \ G 

\ 'Xj \ 

G G X F„ 


F, G 


\ 


; S„ 


G 

0 


V 


G F, \ 

\ X ' 

\ X \ 

\ \ G 

~ \ 

G F, 


(127) 


F i = 


/\ B i B i\ 

B 1 A 2 \ 0 

\ \ \ 

o \ \ \ 

\ \ B > 

, B . B, \/ 


F 2 = 


a 2 b, b, 

\ 

8. A, \ 

X ' \ 0 


0 \ 


\ 


\ B 1 

A B, X A ,/ 


(128) 


in which again A i is as given in eq. (114) with M replaced by M i5 i = 1, 2. 
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Once more we assume the conditions given in eqs. (113) and (117), and make 
use of the results (306) and (307) in Appendix E to solve for the frequencies. It 
is seen that there are 24 distinct branches and we exhibit below six typical fre- 
quencies as the elements of a 3 x 3 diagonal matrix Q* kr , 


. 


PJi, } (129) 


in which P jkr and Q jkr are as defined in eqs. (122), (123), and 


j = 0 , * • • , rij - 1 ; k = 0 , * * * , n 2 -l; r = 0 , • •*, n 3 -l; 8 = n/n l ; cp = 77/n 2 ; \p = 77/n 3 

(130) 

Three additional matrices of the same type as in eq. (129) but with different 
signs accompanying the cosines, complete the frequency spectrum of A(^ 2 ). We 
note that the functional form of the frequencies for periodic boundary conditions 
coincides with that for rigid boundaries. 


We conclude this chapter with two remarks: First, the treatment of harmonic 
lattices described in the opening section is fully adequate for attacking lattices 
with symmetries different from the rectangular model considered here. Second, 
although the linear chain and the plane lattice particles have been assumed to 
vibrate longitudinally and in the plane respectively, it is readily seen that these 
limitations can be dropped if certain simple reductions are carried out in the 
three-dimensional lattice. 
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CHAPTER II 


OTHER APPLICATIONS 

The physics of periodic units is by no means restricted to lattice-dynamical 
problems. In this chapter we shall present also situations which are unrelated 
to calculations with normal mode frequencies, yet require handling by similar 
mathematical techniques. The functions of matrices that arise and the finite 
sums associated with these, are not necessarily connected with the statistical 
mechanics of the systems. The functions appear naturally in most of the situa- 
tions considered. The examples below, selected from diverse fields, will illu- 
strate these ideas. 

The theory of electrical lines is entirely analogous to lattice dynamics theory, 
and one can indeed set up a one-to-one correspondence between the quantities of 
interest in both theories. An example of this is given byBrillouin [4] in the 
treatment of an electrical line analogous to a diatomic linear chain with two 
distinct spacing constants. 

A different example from classical physics was discussed by G. N. Watson 
[27]. He became interested in the evaluation of sums of the form: 


N - 1 


(r) = y ' [cosec (ki9/2) ] ~ r ; 6 - 2 77 /N ; r = 1 , 3 , • • • 


( 1 ) 


k = 1 
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Apparently, these sums occur in a classical treatment of a Bohr type atom: 
electrons, equidistantly spaced, are assumed to move in the same circular orbit, 
with the positive nucleus at the center. The small oscillations of the system are 
found under the assumption of Coulomb interactions and the sums enter the re- 
sulting frequencies. The analogous gravitational problem is the motion of 
satellites rotating in a circle about a planet This is a simplified form of the 
problem of Saturn’s rings considered by Maxwell. 

The sum above, laboriously evaluated by Watson, can be immediately brought 
to the form of a trace of a known matrix, on using the identity 

2 sin 2 (k (9/2) = 1 - cos k (9 (2) 

Introducing a parameter e , which is ultimately made to vanish, we can write 


-r/2 

( 3 ) 


This sum is treated in eq. (90) of Chapter HI. 


The statistical mechanics of a finite plane lattice fully packed with rigid 
dimers (these are pairs of particles connected by bonds) provides an example 


of a situation where normal mode vibrations do not appear. The problem which 
has been treated independently by Fisher [28] and Kasteleyn [29], is to evaluate 


the configurational grand partition function of the system 


z mn = 2^ g mn ( x -y) < 4 > 

m, n 

1 

m+n=-MN 

2 


N ~ 1 

S N (r) = lim V 
e-0 f rri 


k = 1 


11 1/3 

e + cos k 0 

2 2 
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in which g mn (x , y) represents the ninnber of ways of placing m horizontal (or 
x-dimers) and n vertical (or y-dimers) on a lattice with square unit cell and 
MN sites. Thermodynamically x ,y are the activities of the x -,y -dimers re- 
spectively. Physically this problem is a simplified version of a model, including 
also monomers, considered in the thermodynamics of adsorbed films and mixed 
solutions [30, 31 ]. Fig. 1 shows a simple situation. 

X 


Figure 1 


Using topological methods to enumerate the possible configurations, Fisher and 
Kasteleyn succeeded in showing that Z 2 N is equal to the determinant of a two- 
dimensional continuant matrix (or a two-dimensional circulant matrix if the 
dimers are placed on the surface of a torus), as follows: 


Z 2 N = 2 mn 1 1 {x 2 cos 2 mO + y 2 cos 2 ncp) 

m = l h = l 

6 - 77/ (M + 1); 9 = 77/ (N + 1) 
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A simple transformation brings this to the form 


Z 2 

a mn 


1 im 

e -*0 


1 

7 M 


n n {e + a 2 + x 2 cos 2m 6 + y 2 cos 2ncp} 

m = 1 n “ 1 


( 6 ) 


in which 

2 2 2 
a z = + y 

If we denote by Z (e) the expression following the limit sign in eq. (6), we can 
write 


N 

y i (?) 

e + a 2 + x 2 cos 2 m 6 + y 2 cos 2 n cp 

This is a sum which can be readily approximated by using the techniques of 
Chapter III. 

In the two-dimensional Ising model of a ferromagnet, the traces of certain 
matrices represent the thermodynamic functions associated with the system. The 
fundamental paper by Onsager [32] exhibits several such sums, the simplest 
being similar to that in eq. (1), 

N 

S N (1) =1 y cosec (k - 1/2) d; 6 = tt/N ( 8 > 

k =1 


Be 


l 

T M 


logZ(e) 


•L 

m - 1 


in which N is the number of parallel chains (or rows). This sum enters the 
specific heat expression for the lattice of spins. The sum can be treated in a 
fashion analogous to that described for eq. (1). 
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A more complicated type of sum arises in this model when correlations 
between spins located at different sites are considered. Bruria Kaufman and 
Onsager [ 33] show the appropriate ensemble averages to yield sums of the form 



cos [2 7rak/N + S' fe ] ; a = 0, 1 , 2, • • • 


( 9 ) 


in which S' is given, for a quadratic lattice, by 


S' k = cot 


'f 


cosh 2H [1 - sinh 2H cos (27rk/N)] 


sin (2 wk/N) 


( 10 ) 


where H=J/kT andj is the interaction energy between nearest neighbors. 

Similar sums appear also in the spherical model of a ferromagnet proposed 
by Berlin and Kac [34] . 

The theory of random walks on multidimensional lattices has lately attracted 
attention because of their mathematical equivalence with certain physical situa- 
tions, notably the motion of defects in crystals and the theory of spin-wave inter- 
actions [35]. Extensive use is made there of Green's functions, which appear as 
multidimensional sums of the type shown in eqs. (163) and (207) of Chapter HI. 

In a study on the electronic states of a one- dimensional crystal under an 
applied electric field, P. Feuer [36 ] is led to sums of the type 


S N (n) 


N 

2 V ' cos n k d_ 

N / i yi + 2 S cos k 0 
k =o 


6 = 2 tt/N 


( 11 ) 
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These sums represent the coefficients needed to construct the localized Wannier 
functions from Bloch orbitals. This is a slightly more general sum than the one 
in eq. (1) (for r = 1). 

The reverse problem of constructing Bloch orbitals from atomic orbitals, 
i.e., from localized functions centered around the lattice sites, leads to the im- 
portant concept of overlap matrix A [37]. If 4>^= 4>(t - r (i) ) is the normalized 
atomic orbital connected with the nucleus at the lattice point r (£)= x a x + 

: 2 a 2 + ‘ 2 a 3 , then A is defined as follows 

L zi' = JVe < V dr (12) 

The Bom-von Karman cyclic conditions imply here that A is a circulant matrix. 

If instead of the original set {<p^} , one requires a new set {cp^} of orthonormalized 
atomic orbitals, the transformation is 

cp = 0 A~ 1/2 ( 13 > 

Lowdin and al. [6] have evaluated A' 1/2 for a linear chain by using a Chebishev 
expansion method. This approach of linear combinations of atomic orbitals (LCAO) 
is known also as the tight-binding method. 

Molecular theory abounds in systems characterized by repeated units. First, 
the problem of calculating the normal mode frequencies of molecules possessing 
periodicities: e.g., linear and zig-zag chains of atoms [5] lead to matrices 
(continuant and circulant) already considered. More interesting are the calcula- 
tions of the electronic structure of polyenes and aromatic molecules. Lennard- 
Jones and al. [38] have treated these extensively by use of the molecular orbital 
(MO) method, and were lead to continuant and circulant type matrices, the eigen- 
values of which gave the tt -electron energies. Moreover, these authors encountered 
sums of the form 
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S N = 2 


(14) 


N 



+ 2 a ft cos k 6 + /3 2 ; 6 = 2t7/(2N + 1) 


representing the total orbital energy. 

Work along this line has been carried on by several authors [39] , while more 
recently Salem [40], in a study of bond alternation in long polyenes obtains sums 
of the type 


ft 

S 2N + i ( r ) " 2N + 1 


cos r k 9 


Vl + 2t 2 cosk(9 + t 4 


; e = 2rr/ (2N + 1) 


(15) 


In conclusion we remark that most of the authors mentioned above evaluate 
the sums they obtain by passing to the limit N - ® and calculating the resulting 
integral. Lowdin [6] and Gilbert [7] are the exception - they evaluate their 
sums exactly. Salem uses a technique of changing the order of summations, but 
then expands in powers of t 2 and stops with the first term. A generalization of 
this technique is the basis of the mode of calculation presented in Chapter Eli. 

It is to be remarked further that while the conversion of the sums into 
integrals is not critical for large N, for moderate or small number of particles 
this is no longer so. It is therefore important to be able to assess the corrections 
stemming from finite N. Our method has precisely this advantage. 
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CHAPTER m 


METHOD OF CALCULATION AND SPECIFIC RESULTS 


In the preceding chapters we have presented several physical and chemical 
problems that lead to the classes of matrices exhibited there. In this chapter we 
discuss functions of these matrices, in particular analytic functions. Such functions 
appear naturally in the statistical mechanics of the systems previously discussed 
and in perturbation treatments of chemical systems. One matrix function that 
already appears in classical physics is the inverse of ACw 2 ): when external forces 
act on the system, knowledge of A -1 (oj 2 ) is necessary for the complete solution of 
the dynamical problem. The inverse of A (co 2 ) (or the Green's function of A) is 
essential also in the method developed by Montroll et al. [8,9] for the calculation 
of characteristic frequencies of lattices with defects. On the other hand all of the 


thermodynamic functions of the systems previously discussed are essentially 
traces of analytic functions of the dynamical matrix A (w 2 ): e.g., the partition 
function from which all thermodynamic functions can be deduced, is defined by: 


W / 

i- + log f 1 - e kT 

kT V 


or in trace form 



where is the electronic ground state energy of the system. 


61 



The vibrational zero-point energy is just 

E = — — 3r (a 1/2 (a*)}- (3) 

0 2 kT V 

We note that to find the elements of the function of a matrix it is necessary to 
calculate certain finite sums involving the same function of its eigenvalues. Let 
us assume that F(z) is a function defined on the spectrum of A and analytic in 
some neighborhood of the origin, and let k k , k = 1, .... N ,the eigenvalues of the 
matrix A be all contained in the circle of convergence of F (z). Then if T is the 
matrix diagonalizing A, we can write 


where 


Fw = L 

n = 0 


r ( ">(Q) 4 „ 

n ! 


-L 

n = 0 


F (n) (0) T A n T* 1 = TF(A) T 1 
n! 


A kj =A k S k . 


Hence 


N 

t F (A)] kr = Tkj T3rl - 

j = l 


(4) 


(5) 


( 6 ) 


Sums of this type are somewhat more general than the trace-type sums 
associated with the thermodynamic functions, 

N 

2^ F(\.) = 3r {F(A)> ■ ( 7 ) 

j=l 
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It is of interest at this point to present alternative expressions for fini te 
sums of the type discussed. First the contour integral representation based on 
the Cauchy residue theorem 

N ~ 

zL F(Xj)= ih. J F(z) d7 logj9 N( z > dz < 8 > 

j = l J C 

in which 

£> N (z) = det(A N - zl) (9) 

if the k. are eigenvalues of a known matrix A N , or J9 n (z) a polynomial with k 
for roots if no such matrix is known. The contour C has all k. in its interior. 

The representation (8) is useful more for asymptotic approximations than 
for exact calculation of the sums. 

A different representation arises when the sum is converted into an integral. 
Let the L's be real numbers, then one can write for sufficiently well- 
behaved functions F (k) 

F(A) = I ¥(k) S(k - \.) dX. (10) 

j - 00 


where S is the Dirac S -function. 

Assume now that interchange of summation and integration is permissible. 
Then we obtain 




F(k) 



S(k - \.) V dk 


( 11 ) 


We define now 


QW - - 

N 


N 



j=l 




( 12 ) 
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and hence 


N ^ 

^F(\ j )=N | F(M Q(A) dk 

: i - 00 


( 13 ) 


If we denote byU (A) the number of k- } 's which are less than or equal to k , 
it can be shown that 


N 

ri(\) =£ h(^- V 

j=i 

where H (x) is the Heaviside unit function, and 

Q(M =— — 

N dk 


(14) 


(15) 


The function Q (A.) is called the eigenvalue distribution function. It is not difficult 
to show also that 

Q(IM) = 2|\| Q(A 2 ). (16) 


Later on in this Chapter we shall evaluate Q ( A 2 ) for certain special forms 
of A . . The above formalism will be useful also for the conversion into integrals 
of sums other than the trace-type ones. In the following we shall treat both types 
of sums by a method which is different from both the complex and real integral 
representations. Although, for the sake of continuity, explicit reference is 
given to the matrices from which the eigenvalues arise, it will be apparent that 
the treatment is independent of this knowledge. 


3.1 One-Dimensional Sums 

In this paragraph general results regarding circulant and continuant 
matrices will be presented. We have seen that these matrices are associated 
with problems where periodic and rigid boundary condition, respectively, prevail. 
It is not difficult to write down formulas for the matrix connected with free 
boundaries, and also for more general situations — the only requirement being 
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that the eigenvalues have the same functional form. We shall not do this here 
as the method is sufficiently well represented for the cases mentioned. 

(i) Periodic Boundary Conditions 

Instead of treating the symmetric case alone we shall consider here the 
more general case of an asymmetric circulant matrix, defined as in eq. (116) of 
Appendix E 

A = Of) S 1 ' ' ' S N-l>cyc. (17) 


with eigenvalues 


N-l 


v - 1 

- J s. e N ; k = 0, . 
i=o 


. . , N - 1 


(18) 


As shown in eq. (124) of Appendix E , the diagonalizing matrix U is such as to 
reduce the (m, n)-element of F (A), to the form: 


N L — i N 

k-0 


(m - n) k 


(19) 


Therefore we consider in the remainder of this paragraph sums of the type 

N-l 

Sn F( ^ k) eXP 

k=0 

where r is an integer. Without loss of generality we can restrict r to the range 
0 < r < N - 1, as all other cases reduce to the form (20) by periodicity. When 
r = 0, reduces to the simpler trace-type sum of eq. (7). 

As it stands, the sum (20) cannot be computed unless N is small or unless 
F (X,^ has special features that enable the sum to be evaluated in closed form. 

U sing a technique which is the generalization of a method quoted in [ 6 ] , we shall 
convert this sum into a form which is more convenient for accurate calculation or 
for approximate evaluation. 


27ri rk 
N 


( 20 ) 
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Our technique is based on the fact that the A. k ! s have the same functional 
dependence, namely, there exists a function k (9), which generates the eigen- 
frequencies under discussion 

^ k = M 0 k)- (21) 


Here we define 


\(d) 


N-l 

L 

j=0 


s exp(i j O'). 


(22) 


Then clearly 


k. = k 

k 


2 rrk 
. N 


(23) 


The form oik (6) suggests a Fourier series expansion for the function 
F {k (8) } . Therefore we shall assume that such an expansion exists and write 


CD 

F{\(0)> = ^ A j ex P(i j 0 ) 
j =o 


(24) 


with 


r 2n 

v-L f 

J 277 J 


FfA-C#)} exp( -i j 0) dd . 


(25) 


In eq. (24) only positive j ’s are required since k ( 8 ) contains only positive 
powers of exp(i^) and f(z) only positive powers of z. 

Using the expansion (24), we can write 


F (\) = A i ex P 



( 26 ) 
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where the A* are expressible in terms of the A’s: 


A* 

J 



Aj +^N ‘ 


( 27 ) 


The result in eq. (26) is a consequence of the periodicity of the trigonometric 
function exp(i j ( 9 ). We note that the rearrangement of the series is legitimate 
since the analytic character of F(z) guarantees uniform convergence within the 
circle of convergence. 

Insertion of (17) into (11) leads to 


k=0 



j=0 


exp 


/ 2 77 ik(j 
\ 




(28) 


This result shows that a finite summation of this type acts like a filter on 
the Fourier coefficients of FfA.^)}, "sifting" out an infinite number of them. 

In this connection a historical note is of interest: the British mathematician 
Thomas Simpson discovered in 1758 [13] that if a function f(x) has the Taylor 
series expansion 

CD 

f (x) f n x " 

n=0 

and if 


then 


N-l 

V r>0),-L£ f 

k=0 


f £ 2ttU. \ 2ni rk 

e N x e" N 


uj 

^ = x r 22 ( xN > n . 

n=0 
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This is a consequence of the fact that the sum of the roots of unity vanishes. 
Later on, De Morgan [13] extended this result to the roots of x N = -1. 

2?rk 

It is obvious that for \ = e l ~ z , a Fourier expansion of F (z) will not be 
needed in our treatment. 

Returning now to the result (28) we see that the evaluation of the finite sum 
is reduced to the evaluation of the Fourier coefficients and the infinite summation 
of eq. (27). While this may seem to complicate matters, it must be remembered 
that Fourier coefficients decrease with increasing index, by virtue of the Riemann- 
Lebesgue lemma, and the convergence of the series can be quite rapid. Fre- 
quently the approximation 

S* - A r (29) 


will be sufficient for computational purposes since when F (z) is analytic, the 
Fourier coefficients will fall off exponentially with the index. Moreover the 
approximation of eq. (29) is precisely the conversion of the sum (20) into an 
integral (i.e., passage to the limit N - oo ) . This will be a particularly effective 
approximation when N is large. For those problems in which N is not large 
enough to warrant the approximation by an integral, the formula (28) allows one 
in principle to estimate the correction terms successively. 

Procedures similar to those given above suffice to discuss the important 
special case when a is a symmetric matrix. We can write in this case 

f ( S 0 S 1 • ' ' S M S M-1 • ' • S l>cycl. N = 2M 


( S 0 S 1 ' ' ' S M S M S M-1 ' ' S 1 ^cyc 1 . N = 2M + 1 . 


The eigenvalues can be written then 


K = s o + 2 



/ 277k j \ 

s j cos ) + 26 s m cos 


/ 277Mk \ 

\~N~7 


(31) 
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where 


- for N = 2M 
2 

1 for N = 2M + 1 


(32) 


The function k {8) defined in eq. (21) is replaced in the present case by 


M-l 

k(6) = s 0 + 2 L s. cos j 6 + 2 e s M cos M 8 . 
j=i 


(33) 


Then k(d) being an even function of 1 9, so is also F{\(0) } and we can assume the 
Fourier expansion: 


F{\(6?)} = 


= — A + \ A cos n 

2 0 / . n 
n= 1 


e 


(34) 


where 


J JLTf 

F{\(0)} cos n 6. 


(35) 


As for the asymmetric circulant matrices, a rearrangement of the series 
for F (A^) leads to 


N-l 

f<\) = 2Z A i 


j=0 


* cos Izii i 
N 


(36) 


with 


00 

A 't« 


(37) 




and 


A* 


00 

Zv 


j>0 ^ , "j+^N 

-C=o 


(38) 
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( 39 ) 


Then the sum (20) for the of eq. (31) becomes 



If the exponential in eq. (20) is replaced by the cosine term cos 27rrk/N, 
eq. (39) will be replaced by 

K ♦ *».,}• < 4 °> 

In both cases r = 0 leads to the formula for the trace-type sum, 

S„ = K ■ < 4 « 


The results of the preceding analysis can be easily generalized to the calcu- 
lation of sums of the form 


N-l 


k=0 


F(\.) exp 


[" i7r(2k + 1) r 

L N . 


(42) 


with r as in eq. (20) and 



j=o 


exp 


i 


7T(2k + 1) j 

N 


(43) 


are the eigenvalues of the matrix A + , called a skew-circulant in eq. (164) of 
Appendix E. 

The argument leading from eq. (22) to eq. (27) can be repeated with a slight 
modification to yield in this case the result: 

S* = A* ( 44 ) 


in which 


*: = 2Z ( ' 1)i a '*‘ n ' <45) 

j=0 
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Similar results can be obtained for other values of r in eq. (42). If on the other 
hand the first row of A + can be written as 


( s 0 Sj . . S M , S M _j ... Sj) 


A + : 


( s o s i ' ' ' S m’ " S M ' • • * s iW- N = 2M + 1 


skew- cycl . 
skew- cycl . 


N = 2M 


(46) 


then the eigenvalues are 


M-l 

j=i 


\ o > 7T(2k + 1) j - 77 

r s o + 2 / S j cos — — + 2 e s M cos 


N 


(2 k + 1) M 
N 


where 


r 


Z- < 


— N = 2M 
2 

1 N = 2M + 1 . 


The Fourier expansion then contains only cosine terms 


(47) 


(48) 


A = 

n 


= — f F{\(<9)} cos n& dd 

77 J 

J o 


(49) 


and the sum 


N-l 

k-0 


F(A.. ) cos 


7(2k + l)r 

N 


(50) 


with r as in eq. (20) , can be represented by the expression: 

[*♦<-.] 


(51) 
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in which 


00 

A* ( _1 ) J A j N 

j=i 

> 

00 

A !>„ -J2 ( - 1)iA ^i» 

j =0 J 


( 52 ) 


An analysis of the preceding exposition shows clearly that the following general 
result can be stated: 

From a given eigenvalue-generating function k (6) and a given analytic func- 
tion F(z)— which is defined on and its domain of convergence contains the set of 
values {\(<9), 0 <6 < 2 n}— one can construct a whole class of finite sums by 
choosing N values 6 , such that e 10k be the roots of an algebraic equation of 
the N th degree. Since the set of Fourier coefficients for this class of sums is 
unique, the values of these sums will depend solely upon the filtering properties 
of the <9 k 's. 

We proceed now to apply the formulas developed above to cases of general 
interest. The discussion will be restricted to the asymmetric and symmetric 
circulants. 

(i) The Inverse A ~ 1 

The function F (z) = z' 1 , strictly speaking, is not analytic in the neighbor 
hood of the origin, yet the formulas developed above can still be applied if the 
integrals associated with the Fourier coefficients are taken as principal value 
integrals, whenever needed. 
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In the asymmetric case the integrals to be evaluated are 


A = 

1 2tt 


.277 


i jl9 


d& 


'o s o + s i e + s 2 e + • • • + s N . 1 e 


i (N- 1 


2- f 

2 77 1 J z j+l 


dz 


(s 0 + s x z + . - . + z N1 ) 


(53) 


where c is the unit circle. The contribution to the integral from the pole at the 
origin is the coefficient of z j in the Taylor expansion of 


[P(z)]' 1 = [s 0 + SjZ + . . . + s N _ x z^ 1 ]' 1 A( 

In the form of a determinant this contribution is 


p) 7 i 


( 54 ) 


j=o 


A< p > 

j 


(-1) J 

( s n) j+1 


S 1 s 0 
\ \ 

S 2 \ 

\ ' 


\ o 


\ 


\ 


: % 


(p) = 


( 55 ) 


\ \ \' 

N ' \ ' \ 

\ x N N S 
. \ \ \ 0 

Sj-~ “ “ “ S 2 S 1 

Unless a computer is used for the numerical evaluation of the determinant, 
this expression for A ( . p > is not very practical. An equivalent form can be found 
for A ( j p) provided the roots z k of P (z) are known. Suppose that P (z) has no 
multiple roots (the extension to the more general case is trivial). Then we can 
write 


N- 1 


p O) = s o | | (1 + \ z) 


(56) 


k=l 


with cr = - l/z. . In terms of the a k , the value of A. i s calculated to be 
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(57) 


Aj = (-l) j 


N-l 

L 

k= 1 



Since the only dependence on j is through the term c^ + 1 we may find an explicit 
expression for the A* (p) : 


N-l r+l 

* <P) = (-!)' T -1 ^ 

kTi p»L±\ 1 “(-°k> N 


(58) 


where we have summed the geometric series involved under the assumption that 


kj < i- 

In particular, A (p) * will coincide with A* whenever P (z) has no roots 
inside or on the unit circle. There exist several sufficiency conditions for this 
to happen. One condition given by Landau [14] , is 


s i > S 2 


> S 


N- 1 


> 0 ; s! 


s u /s o 


(59) 


A second simple condition can be obtained by using Rouche's theorem*. For 
example if 


*\\ + I s' | + ... + |s'J <1 


(60) 


then the roots of P(z) will lie outside of the unit circle provided 

(1+sj + s' + ...s' e )>(S e + 1 + ...+S N1 ). 


(61) 


In particular this is true if 


N-l 



i=i 


(62) 


* Rouche’s theorem states that if f(z) and g(z) are analytic within and on a closed curve C, and if 
jf ( z ) | > | g ( z ) | on C, then f (z) and g (z) + f (z) have the same number of zeros inside the region 

bounded by C. 
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An exact result can be written down immediately for the case s k =0, k > 3 if 
it is assumed that only the pole at the origin contributes. Then the easiest way 
to evaluate the A. 's is to use their determinantal expression of eq. (54). Using 
the value of this determinant given for the appropriate form of eq. (15) in Ap- 
pendix E, we find 



Therefore the sum 



on evaluating A* of eq. (28) with the A. ’s of eq. (63), becomes 



This expression has a finite, well-defined limit when N - oo . 

The elements of the first row of A -1 are obtained from eq. (65) by letting 
r take the values 0, 1, • • • , N - 1. These elements characterize A -1 completely 
since a function of a circulant matrix is also circulant. 

It is clear that if all the s k ’s vanish with the exception of (any) two of them, 
the sum again is easily evaluated in a closed form. In some instances it is more 
convenient to use a full Fourier expansion. For example if the ’s are of the 
form: 
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2tt(N- 1) k 
K = S 0 + S 1 « N + S N-1 e N 


. 277k 

1 - 


. 277k . 277k 

= S 0 + *1* N + S N-1 e N 


Then the appropriate k(6) is 

k(9) = s 0 + s x e ie + s N _j e~ ld 




k = 0 N — 1 > 


( 66 ) 


(67) 


and the Fourier coefficients associated with the sum 


1 V 

Sn “ N / . 

k = 0 


N-l 

e N 




are 


A = 


1 C 2n e~ ing d8 1 f dz 

277 Jo s 0 + Sj e iG + S N _J e~ id 2jtx J c z n (s 1 z 2 + s Q z + s N _j) 


n > 0 


j_ r 2n 
n ' 2 - 1 


A. 

n > 0 


e ind dd 


o s Q + Sj e id + s N _ 1 e 


— f 

s ~i 0 2771 J 


z n dz 


s, z x + s n Z + S-. 


N-l 


( 68 ) 


(69) 


The values of these coefficients, which can be calculated exactly, depend on the 
location of the zeros of the denominators with respect to the unit circle. Since 
in Appendix D we exhibit an alternative way of computing such simple sums 
exactly, we shall neither calculate here the integrals (68) and (69), nor perform 
subsequently the summation for S N . 

Similar expressions will be obtained now for the symmetric case . Let us 
start with the expression (35) for 


r 2n 

x =± f 

" 77 Jn 


cos nd dd 


M-l 


s 0 + 2 


^ ' s. cos j 8 + 2e cos Mi 9 
j = i 


-I f 27T cos n 
_77 J 0 V co 


e de 

cos 8) 


(70) 
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The polynomial P M (cos 9 ) can be factorized such that 

M 

P M (cos 9) = s Q a. cos 0) ( 71 ) 

j = i 


where again a. are the reciprocals of the roots of P M (z). 
discussion to the previously investigated case of 

M 


T hi < w 


If we restrict the 


(72) 


which means |cr. | < 1 for all j , then the integrals A n can be evaluated by 
partial fraction decomposition: 



Hence the appropriate sum of eq. (40) can be summed explicitly 


L 


-rf (1 -cr 2 ) 

J-l ' ] 


2 ) 1/2 p , / 1 \ ! 
j 


- 1 \ !-{[(! - cr }) 1 ' 2 -1 ]/ a .}“ 


(1 - a}) 1 ' 2 - iy /(l-q2)l/2_i 


M-r 


cr. 

J 


(74) 


A number of exact results are available also in this case. For example if MS) 
is of the form 


*The value of this integral is exhibited in Appendix B 
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(75) 


k(9) = S Q + 2Sj cos ti 

then again using the integral ( 4 ) of Appendix B, we can write 


7 / 

•'n 


r 


a2tt 

cos nd dd 


1 

s 0 + 2s, cos 6 



1 / /s o - 4s i " s oY 

^ s o - 4s ? - s o 


\ 2s, ) ' 

2s 

-4s 2 ' 1 > 

l 


< 1 


< 


(- D r 




/ §2 _ 4 S 2 
/ Sq tij 


//s 2 - 4 s 2 + s 0 V 

7s 0 - 4s 2 + s 0 

\ 2s i / 

2s i 


< 1 


Then the first case in (76) leads to 


S =— 
N 2 


/ _ 4 — s 

27 


/ s o - 4s i - s 0 y- r 

2s, 


/? — 4s 2 


1 - 


2 s, 


(76) 


(77) 


The second case leads to an identical form since the value of a finite sum is 
independent of the mode of summation as long as no terms are neglected. 
Another case which admits of an exact treatment is the following: 

. _ 2nk „ 4rrk (78) 

\=s 0 + 2s 1 cos — + 2 s 2 cos—. ' 

Instead of evaluating the Fourier coefficients it is simpler in this case to 
relate the sum associated with (78) to the sums associated with the eigenvalues 


A.* given by 


A^ = - 2 /ij x ± +2 /7~ 2 cos^ 


(79) 


where 




/s 2 - 4 


S o S 2 + 8 S 2 


4s„ 


(80) 
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and = \ k . It is easy to see then that 


%.=-r-r is -si- 

X - X 


(81) 


(2) Other Functions of A 

We have seen that the calculation of the inverse of the general circulant 
matrix already presents great difficulties. It is therefore not surprising that 
the evaluation of more involved functions complicates the computation consider- 
ably. Hence such cases will be best treated by numerical methods (assuming of 
course the elements s . to be known). But even for numerical treatments the 
general discussion above is of value, since it exhibits clearly what are the 
dominant terms and therefore provides a practical scheme of computation. 

For certain special cases some results can nevertheless be given in terms 
of known functions. These special cases are treated below. 

(i) A v , v_£ 0, 1, 2, ••• 

The elements of A v give rise to the sums 


N-l 

*-z: 

k=0 



+ s . 


277-ik 47Ti k \ v 

e N + s 2 e N ) 


- jZLi 
e n 


(82) 


where we have assumed 

k(d) = s Q + Sj e ie + s 2 e i20 . 


(83) 


The Fourier expansion of [A. (i 9)] v can be effected as follows: 


(s 0 + Sj e iS + s 2 e i2S ) v = Sp [1 - 2xz + z 2 ] v 

where 


(84) 


X 


2 




(85) 
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If it is true that |s 2 /s 0 | < 1 and |s 1 / s_ I < 1. the following expansion is known 
[15] to exist 


(1 - 2XZ + Z 2 ) V = ^ C n V ( X ) Z " 

h = 0 


( 86 ) 


where (x) are Gegenbauer polynomials with the explicit representation 

W 


c; (x> =<-D" 2] („?.)( V) (2x) " 

m = 0 


- 2m 


(87) 


Therefore the Fourier coefficients of eq. (84) are 

,2tt 


n 2 


hi 


(s Q + s 1 e 1 ^ + s 2 e l2(9 ) v e * in61 d 6 


S 2 


vn/2 


c v 


2)/ s 0 s 2 


( 88 ) 


Inserting these values in previous formulas we obtain the sum S,^ of eq. (82). 
When A is symmetric the pertinent sums are of the form 

N- 1 

l ~ 2 77 k Y 2 77 r k 

(89) 


cos ^ir) 

k = 0 


CO S 


N 


The evaluation of this sum proceeds as follows: we write 


(s 0 + 2 Sl cos ey = ( S 0 y 


<2su 


(l-2/3cos0 + /3 2 ) 

(i +/e 2 r 


(90) 


where fi is that root of the equation Sj x 2 + S|) x + s ( =0 which stays finite for 

Sj 0, 


P 2 s, + 1/(25,) 1 


(91) 
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To find the Fourier expansion we require \fi\ <1. Assuming this condition 
to be satisfied, we can use tables [16] to write down the result 


r 2n 

K =± f 

" - Jo 


(s Q + 2s J cos dy cos n0 dd 


o \ 2^ n r(-v + n) , 

2 ) “7", a ~ ~ F( ' V ' -v+n, n + 1: p 2 ) 


\i + P 2 j r (-v) r (n +i) 

where F denotes the ordinary hypergeometric function. Then 


(92) 


l s ° 1 

1 1 

\1 +fi 2 ) 

r (- v ) 4ro 1 


/3 r+nN V (-v + r + nN) 


F(- v\ - v + r +nN; r + nN + 1; fi 2 ) 


r * 0 


/S (n + l)N-r r( _ v + (n + 1)N _ r) 

+ — — F( -V, -v + (n + 1)N - r ; (n + 1)N - r + 1; / 3 2 ) 


Sn = 

r = 0 


s 0 \* 


,1 +/3 


F((n + 1) N - r + 1) 

F( -v, - v , 1; /3 2 ) + 


(93) 


V 1 /S nN F( -v + nN) M .. .. o 2 . 

/ — - F(-v ; -v + nN; nN + 1 ; p- 4 ) 

Z__i F(nN + 1) 


F( - v) /—i F(nN + 1) 

n- 1 

For v = ± 1/2, these results are somewhat analogous to those obtained by 
Lowdin, Pauncz, and de Heer [ 6 ] . 

(ii) Exp (rA) 

The sums involved for the asymmetric case are of the form: 


(94) 


N-l 

exp|T 

k — 0 


s o + s i e 


2rriU 

N 


2'77'i k r 
N 


(95) 


The Fourier coefficients are found by using the well-known expansion: 

ie y a" 

- ae = / — e in 6 (96) 

L 1 n! 

n=0 
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Hence 


1 T s ^ 

= 2 ^ e J e~ inS exp (rSj e ld ) d<9 
0 


C^-s ) n 

1 e o 


( 97 ) 


Therefore 


S N = e 


Ts o V 1 (rs i ) 


r+nN 


Z I (r + nN)! 

n = 0 


(98) 


Similar results can be easily derived from the above for sin (tA), cos (tA) and 
c rA (c a constant). For the symmetric case the pertinent sums are of the form: 


N“1 

E e 

k = 0 


xp |r^s 0 + 2Sj cos^jj 


CO s 


27zrk 

"N 


(99) 


The Fourier coefficients are 


r 2 7T 

1 f Ts 

A n =— I e o exp{2SjTcos0}cosn#d# 

77 o 

l c 27T 

= e 0 _ exp {i (- 2s ir) cos 6 } cos nd di 9 (100) 

7 T J 1 


Using the expansion [16] 


exp (i z cos 6 ) = J 0 (z) + 2 i n J n (z)cosn£ 

n= 1 


(101) 


we immediately obtain 


A n = 2e S ° I„ (2 S j t); n = 0, 1, . . . 


( 102 ) 
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where I (x) is the modified Bessel function of the first kind. Then 

n ' 


N 

r40 


uu 

= eT ’°Z] (2s i*’ I (n .ON-- < 2s . T >}(103) 

n = 0 


S N = e ™°Uo (2Si r) +2 T* (-l) nN I nN (2 Sl r) 

l 


(104) 


Again the functions cos (r A), sin (tA) and c tA can be evaluated by using the 
results above. 


4. Calculation of the Eigenvalue Distribution Function 

Here we want to evaluate the function (]( w 2 ) defined in eq. (12) for two cases. 
For a monatomic linear chain the eigenvalues are 

o?. - a + 2b cos — 1 ; j = 0, . . . , N - 1 (105) 

J N 


We need therefore the Fourier expansion of the function b (a? - a + 2 b cos 9). 
Following the procedure of Lighthill [17] for generalized functions, it is easy 
to show, by making use of the transformation properties of the S -function, that 

r 


l 

2-77 |b| | sin <p| 



cos m(p cos 


md 


} 


S(co 2 - a - 2b cos 6) = ] 


a — 2 |b | < a? < a + 2 |b| 


(106) 


0 elsewhere 


In the above we have assumed for definiteness that a > 2 |b| and b < 0, the 
case most frequently met in practice.* Then tp is given as the principal branch 
of the function 


*0ther cases can be treated with similar ease. 


A 

V 


83 



(107) 


(p - arc cos 


co 2 - a 


2b 


i.e., 0 < cp < 77 . 


Writing 


co 2 = a - 2 |b| ; = a + 2 |b| 


(108) 


the above can be rewritten as follows: 


j (CO 2 - Cl?) (cl)^ - CO 2 ) 


\ V" 1 

1+2 J cosmcpcosm # 


77 


S(O) 2 - a- 2bcos 0) 


m=l 


O ) 2 < OJ 2 < CO 2 

s L 


(109) 


elsewhere 


Performing the finite summation over j , we finally find: 


Q(^) =< 


00 

— 4 — J 1 + 2 y 1 cos mNcp L ; a> 2 < a? < a> 2 

TT /(w 2 _ W 2 ) (aj£ _ w 2) "1 m=l J 


(110) 


0 


els ewhe r e 


On comparing this result (for w. = 0) with that of Montroll et al. [11], we see 
that the infinite sum represents the correction terms* which depend explicitly on 
N, the number of particles in the chain. 

It is of interest to note that Q(oj 2 ) can be also written in the form 


r 


Q(^) = < 




CO 2 - OJ 2 ) (co 2 - co 2 ) k = 


elsewhere 


00 

S (Ncp - 


27rk) 


(111) 




*When used as the kernel of an integral. 
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We note that the first term in the eigenvalue distribution function depends 
solely on the form of the eigenvalue generating function, while the correction 
terms depend also on the particular set of 6 k ' s associated with given boundary 
conditions. Since the first term is the more important one we do not explicitly 
exhibit correction terms for other boundary conditions. 

We quote below the corresponding results for a diatomic chain the eigen- 
frequency generating function of which can be written as 


\ + (6) = a) 2 - a ± V0 2 7 y 2 cos 2 6 


( 112 ) 


Then the Fourier expansion of the function 


F(0) - Hk + (9)} + S{\_(£)} 


( 113 ) 


is given by 

r 


a? - a! 


c ~ 

— — ^ 1 + 2 ^ ' cos2mcp 

I t m=i 


F(^) 


77 |(a; 2 -a) 2 — /3 2 j 1/2 I/3 2 +y 2 -(co 2 -a) 2 I 

for a. - I/3 2 + y 2 1 1 ' 2 < a? < a + j/3 2 + y 2 ] 1 ^ 2 
0 elsewhere 

with cp given as the principal branch of 


cos 2m0 * 


( 114 ) 


cp = a rc cos 



( a ? - a) 2 - /3 2 


Then the frequency distribution function Q (or 2 ) 

2N-1 

^ (a)2) = h L [ 8 {X + ( o' . ) } + S{\_(0. )}] 

1=0 


( 115 ) 


( 116 ) 
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is given by 


r 


\a > 2 - a! 


1+2 y cos 2mNcp 

m=l 


=< 


77 |(oj 2 -a) 2 _ /3 2 | 1/2 |/3 2 + y 2 - (c«J 2 + a) 2 1 1/2 

for a - |/3 2 + y 2 | 1/2 < oj 2 < a + |/3 2 + y 2 | 1/2 


0 elsewhere (117) 

It is apparent that these results can be extended to more general situations, 
the sole requirement of this method being explicit knowledge of the real roots of 
the eigenvalue generating function (if such a one exists) . 

We conclude this section by noting that if the function F {A. ( 9 )} can be ex- 
pressed as a sum of products of simpler functions whose Fourier coefficients 
are known, then repeated use of Parseval formula will lead to the desired Fourier 


expansion, and subsequently to the required sums. 

(ii) Rigid Boundary Conditions 

In this section we treat sums associated with the matrix 


A = 



(118) 


shown in eq. (2) of Appendix E. 

Before proceeding to the actual calculations we extend somewhat the scope 
of this discussion by including matrices (not necessarily symmetric) of the type 


A' = 



0 


b 


i 



\ \ 


\ ' 

\ \ 


\ 

\ 


\ \ 


,\ \ 
N- 



(119) 
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where it is assumed that for all i = 1, . . . f N - 1 

b. c. = const = b 2 . (120) 

In Appendix E the discussion following eq. (15) shows that a' and a are connected 
by a similarity transformation. Because of this property we will work henceforth 
with A exclusively. It can also be verified that the elements of an arbitrary 
analytic function F(A') are related to those of A as follows: 


r 


c 

C . - . 

c , 

n 

n + 1 

m- 1 

b 

b , . . . 

b , 

n 

n+l 

m- 1 


1/2 


x t p ( A )Ln m > n + 1 


= < [F(A)1 


m = n 


( 121 ) 


b b , . . . b 

m m+1 n-1 

Cr^ p 

_ m m+1 ‘ ‘ ’ n -1 


1/2 


X [F(A>] n > m + 1 


Since the eigenvalues of A are 


k = a + 2b cos — ; k = 1, . . . , N (122) 

N + 1 

their generating function will be 

k(6) = a + 2b cos 0. (123) 

The elements of an analytic function F(A) can be written as 


N 

Ntl ra 


7r(m - n)k 7r(m + n)k 1 

Ceos - cos / (124) 


N + 1 


N + 1 
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where use was made of eq. ( 6 ) and the form of T shown in eq. (12) of Appendix E. 
Assuming as before that F{/V( 6 >)} can be represented by a Fourier series 


^ ^ 

F{\(6>)} = — + L A n ® ^ ^ 

n= 1 


(125) 


where 


= i r 

^ J 


F (X(<9) } cos n0 d<9 


(126) 


it is seen that the F(^ k ) result when Q is set equal to 77 k /(N + 1). Substituting 
the value of F(\ k ) from eq. (125) into eq. (124), we obtain 


[F(A)1 


— L [ 

2(N + 1) 4^ I 


77(m - n)k 77(m + n)kl 

COS - cos > + 


N + 1 


N + 1 


♦ _!_ r a, » 7>. 

2 < N + 1> fei 1 l N + 1 


^-1 f 77(m - n - j)k 7T(m - n + j)k 

-cos — 


N + 1 


7r(m + n - j )k 

-cos 

N + 1 


” ( m + n + 3 ) k \ (127) 

N + 1 

J 


To evaluate the finite sums appearing in this equation we note that they are 
all special cases of the prototype sum 
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77 S k 


L 

k = l 


COS - 


N + l 


= NS 


s t 


± 2 r(N + 1) 


+ I [(-D"-l][l-S,. 12 , (N4l) ] 


2 + s - i(2r + 1) (N + i) 


( * ^s, ± 2 r (N + 1) ^s, ±(2 r + 1)(N + l)^s,±2r 


(128) 


where r = 0,1,2 The use of this formula in eq. (127) results in 


[F(A) ] mn - — [B m _ n + ®-(m- n ) ®rn+n ®-(m+n)-l 
m > n 


(129) 


where 


B j = £ \ . 2,<»„> < 13 °> 

r = 0 


and a coefficient A ^ is zero if -i is negative. The restriction m > n is dictated 
by the symmetry of A (which in turn implies the symmetry of F(A) ). 

In the following paragraphs we shall give some specific applications of the 
formula of eq. (129). 

1. Calculation of A' 1 

We have already found the Fourier coefficients needed here in eq. ( 4 ) of 
Appendix B. When their expressions are inserted into the definition of the B's, 
the resulting series are geometric and the matrix elements can be written in 
closed form as 
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m >n 


(131) 


(A' 1 ) 


mn 


V 2n - 1 V 2 (N + 1) _ y2m 


Uv m + n i_y2(N+D 


where 


/a 2 - 4b 2 - a 
2b 


U = /a 2 - 4b 2 


(132) 


2. Calculation of a v , v ¥ 0, 1, 2, . . . 

If v as obtained in eq. (132) satisfies 

t V | < 1 (133) 


then the Fourier coefficients are given by eq. (92), on identifying V with /3, s Q 
with a and s ; with b : 


A 


n 


/ a \ v 2T T\-v + n) 
\1 + V 2 J f(- v) r (n + 1) 


F(- v, - v + n, n + 1; V 2 ) 


2( -by V n-V F( , - +n) F(-74 -V + n, n + 1 ; V 2 ) ; n > 0 . 

r( - j/) r (n + i) 


(134) 


Then 


= (- — - J Y" 1 F(-v, -^+m-n + 2r(N + l), m-n + 2r(N + l) + 1; V 2 ) x 


mn I v 


yj r(-v) i Z — i 

r = 0 




P( -v +m -n + 2r(N + 1) . ... .. 

X v ' ym-n + 2r(N + l) + 

r(m-n + 2r(N + l) +1) 
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00 


+ J F( - v , -v + 2r(N + l)-m+n, 2r (N+l)-m+n + l; V 2 )13 m + n + 2r (N + l)) v . m+n+2r ( N+1 ) 


r = l 


P(n-m + 2r(N+l) + 1) 


- ^ F(-v, -i'+m+n+2r(N+l),m+n+2r(N + l) + 1; V 2 ) — ^ +m+n + ^ r (^ + 1) v m+n+2r ( N+1 > 

Z — ' r(m+n+2r(N + l) +1) 


r=0 


CD 

-E r( 

r=l 


T( -v + 2r (N + l)-m -n) | 

-v, _v + 2r(N + l)-m-n, 2r(N + l)-m-n + l; V 2 ) _ v 2r < N+1 ) m -n ^ 

F(2r(N + l)-m-n + l) 


(135) 


3. Calculation of exp (rA) 

The Fourier coefficients are given by 


A n = 2e Ta I n (2br) ; n = 0, 1, 2, ... 


(136) 


Then 


[e rA ] = e Ta 

J mn 


m > n 


UU 

y ^m-n+2r(N+l ) (^br) 


r=0 


UU vaj 

+ y~' i I 2r(N+l)-(m-n) ( 2bT ) ~ I m+n+2 r (N+l ) ( 2br ) 

r - 1 r-0 


UU 

~ y *2r(N+l) - (m + n) ( 2br ) 

r=X 


(137) 


We can write down also the following result for the sum 


N / 7Tk\ 

s = -L_ V e T V a+ 2 bC ° S ^)cos^i 
N N+l Z—j N+ 1 

k = l 


(138) 
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We can restrict m (because of periodicity) to the range 0£m< 2 ( N + 1) • Then 


I °° 

s; = e Ta Jl 0 (2br) t 22 (2b’-) 


U 1 


00 

1.1 


(N+l)+m 


(2 br) 


e 

N + l 


I ft (2br ) S + 

I 0 v s m, even 


L 

*-ta- 


< 2bT > 


“ [2] I 

Z, X - + 2^ ( 2b ^) + 7] Va* (2br) j> 

^ = 1 - - 4*1 J 


m * 0 


(139) 


This result simplifies for m = 0 to: 


Vo = eTa ^o ( 2bT ) + 2 ^ 2 4( N + 1 ) ( 2bT Y - 


_e cosh (2br) 

N + 1 


(140) 


3.2 Two-Dimensional Sums 

In this section we present results relating to sums arising only from periodic 
and rigid boundary conditions. Other cases, whenever tractable, can be treated 
in a similar fashion. 
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It is clear that as we proceed to higher dimensions, the computational diffi- 
culties multiply. Therefore fewer specific results will be available. Nevertheless 
our technique will still be advantageous for numerical calculations. 

The matrices A touched upon here will be two-dimensional matrices, i.e., 

Mx M matrices the elements of which are N x N matrices. To specify the elements 
of a (regarded as an MN x MN matrix) one needs four indices, or two vector 
indices i = ( i , i 2 ) » j-(j x » j 2 )- The first component of the vector index repre- 
sents the row (column) of blocks in which the element is located and the second 
the row (column) location inside the block. 

If F ( z) is an analytic function of z over the entire eigenvalue spectrum of A , 
then the generalization of eq. (6) is as follows: 


< 141 > 

i 

where h- are the generalized eigenvalues of a (namely, M diagonal NxN -matrices) 
and T is the (generalized) matrix which brings a to diagonal form. 

i) Periodic Boundary Conditions 

In this subsection we shall consider sums arising from eigenvalues connected 
with asymmetric and symmetric two-dimensional circulant matrices. 

The general asymmetric case leads to eigenvalues A. of the form: 


K 


mn 


M-l N-l 




m = 0, . . . , m-l; n = 0, . . . , n-l 
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Then, using the form of T given in eqs. (212), (216) of Appendix E 


[F(A')] pq 


1 

MN 


M-l N-l 




exp 



m(Pj - Qj) 
M 



Pj. qj = 0 M-l, p 2 , q 2 = 0 N - 1 • (143) 

Eq. (143) shows that we are concerned here with sums of the type: 

2ni rm 2ni 

F(\ r £) e” M e” N (144) 

where m , n are integers which can be restricted to the ranges 0 ^ m < M - 1 ; 0 < 
n < N - 1 because of periodicity. For m = n = 0, S reduces to a simpler, 

MN 

trace-type sum. We define now the frequency generating function A. (0, 9 ) by 


M-l N-l 



He, 9) 


M-l N-l 

L L 

j=0 k=0 


C e i(j0+k9) 
jk e 


(145) 


Clearly 


X , = X. 1^1, (146) 

V M N / 

The analytic function F{\(0, 9 )} is assumed to have a double Fourier series 
expansion 


FO(£, qp)} = 


00 00 

A e i(p<9+q<P) 
pq 

p = 0 q=0 



(147) 
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where 


A pq = _L_ I f F{\(0, cp)} e-^+^dtfdtp. (148) 

4tt 2 Jo Jo 

Again, only positive p, q are considered, since \(d, q>) and hence F{\(6 1 , cp )} 
contain only positive powers of e 10 , e i(p . In cases where A ( 6, <p ) contains nega- 
tive powers of e l£> and e 1 ® , the full expansion must be used. 

Utilizing the expansion (147) we find 


OU CD 

f c\<> = £ £ 


P=0 q=0 


A e 

pq 


27irp 



M- 1 N-l 

■LL 


p=0 q=0 


A* 
p q 


.277_Tp 



(149) 


where 


A* = 
A pq 


CO CD 

LL 

j=0 k = 0 


P+jM.q+kN. 


(150) 


Again, the analyticity of F (z) insures the legitimacy of rearrangement of the series. 
Insertion of (149) into (150) leads to 

Ah = A *..„ < 151 ) 
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The discussion preceding and following eq. (29) concerning the properties of 
the Fourier coefficients, is pertinent also here and need not be repeated. 

The general symmetric case leads to eigenvalues \ r ,£ of the form; 


+ s ok c ° s ^ + 2 , > s on cos 2 ^ + 2 £ S J° COS ^ 

k= 1 j=l 


M- 1 N- 1 

+ 4 y ' S. k COS “ — iL_ COS 

j=l k= 1 


M-l 


2 ILLI cos T S. N — 277j r — 2vfl * 

M' N' Z_i JN 


i=l 


!. k , COS COS 

M' N' 


+ 2es M0 cos 


2ttM 


M' 


N- 1 
k = l 


277k -i 277Mr 


S Mk cos cos 


N 1 M 1 


where 


. 277Mr 27 tN-^ 

+ 4e rj s MN cos cos 


(152) 


M 1 N l 


r 


e = 


r 


M' = 2M 


5 V = ^ 


1 M' = 2M + 1 


N' = 2N 


(153) 


1 N' = 2N + 1 


We obtain the frequency distribution function k(6, cp ) by substituting 8 - 27rr/M' 
cp = 277^/N' in eq. (152). Then 
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CD 00 

F{^(/9, cp} = i A^ + 

A p0 COS P 6 + j 2^ A 0 q COS ^ 

p-l q=l 


'LL 

P=1 Q- 1 


A pq cos p9 cos qcp 


(154) 


with 


277 -277 


A = _ 

Pq -rr2 



F{k(9, cp)}cos p&cos d£dcp . 


o •'o 


(155) 


Then 


M'-l N'-l 


F <vt>= £ 


* 277pr 2710^ 

i COS — cos _ 

pq 


M f 


N # 


p=o q=o 


(156) 


where 


00 GO 




~ Aoo + A °- kN ' + 2^Z A j M '-° + 2Z A iM'> kN ' 

k=l i=l j=l k = l 


W 

A 0 q = — 2~! A 0’ q + kN 
_ . « k=0 


q > 0 


00 00 

' V ” 1 


• + 


/ , 7 , A iM # , q+kN' 


j=l k=0 


y (157) 


00 00 00 

+ 2Z 2--j A P + i M, ' kN ' 


p>0 i=° 


j=0 k=l 


00 00 


K -'LL ^ 

p,q>0 


Pq / , / , P + jM', q+kN 

j =0 k = 0 
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Then a sum of the type 


i “A 1 sw 2™ 

= L w‘“' 

r=0 


27rrm 2iAjx\ 

e n' 


(158) 


in which m, n are as in eq. (144), can be rewritten with the aid of eq. (156): 


S M , - = - A* 

™ N a m , n 


(159) 


If the exponentials in the sum (158) are replaced by the respective cosine terms, 
the modified sum is given by 


c _ * 

m'n' " 4 


1 1 .* 


'-m, N'-ri + A M'-rn,n +A m,N'-n + Am , n} 


For the trace sum with m = n = 0, we obtain: 


(160) 


^M'N' A 00 ' 


(161) 


We proceed now to calculations with specific functions F(z). 


1. The Inverse A' 1 

In the p vvmmetric case the integrals to be evaluated are 


n 27T -a 

e' ip 5 e‘ l 

0 M-l N-l 


A Pq = 


iqC P dd d cp 


j=0 k=0 

1 r r d^dz 

(2Tri) 2 J c » Jc 


kCp 


M-l N-l 




“ zP “ H L s “ 


(162) 


z J ^ 


j=0 k=0 
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where C' and C are unit circles in the £- and z -planes respectively. It is 

obvious that even though some formal results could be written down for the 

general case (e.g., a determinantal form for the contribution to A from the 

pq 

poles at the origin), these would have no practical value. On the other hand if 
only a few of the s jk do not vanish, certain exact and asymptotic results can 
be still derived. We shall not pursue this possibility here as illustrations will 
be given for the more interesting symmetric case. 

In the symmetric case we treat only the case 

\(6, cp) = a + 2b cos 6 + 2c cos <p + 4d cos 6 cos cp (163) 

where the integrals to be evaluated are given by 

A =J_ r C cos p e cos qcp dd dcp (164) 

7 t 2 J 0 J 0 a + 2b cos 6 + 2c cos cp + 4d cos d cos cp 

The following conditions are imposed on the coefficients: 

|b| > 1 c | > 2 |d| ; b, c, d < 0 (165) 

It is clear that one integration can be immediately performed to yield an elementary 
function, and we shall do so in the sequel. Here we find it more convenient not 
to use contour integrals, the analysis being simpler in the real domain. Since 
for arbitrary p,q an exact result is not available, we give first the value of A Q0 
where such a result exists. In Appendix B we have evaluated A Q0 for different 
ranges of the parameter a . Here we exhibit only one such result, eq. (22) there: 
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1 


1 


(166) 


A - 4 

00 


K 


: 2 _ 4d 2 ^(u - 1) (v + 1) 


' / 2(u - v) 

Y(u - 1) (v + r 


where 


a + 2 |b| _ a - 2 jb| 

2 ( | c | - 2 |d| ) ’ V “ 2(|c| + 2|d|) 


(167) 


and K (k) is the complete elliptic integral of the first kind the modulus k of which 
is given by 


k 2 _ 2(u - v) 

" (u - 1) (v + 1) 


(168) 


This generalizes results in the literature [11], including recent ones of Mathews 
and al. [18] . 

Appendix B also treats the asymptotic approximations of the general integral 
A pq . For the case d =0, exact results can be given for all A pq by proceeding as 
follows: 


A 

pq 


n 77 

cos p 8 cos q 9 dd dcp 
a + 2b cos 6 + 2c cos cp 


(169) 


Under the same restrictions as in eq. (165) with the additional condition |a| ^ 
2(|b| + | c | + 2 | d | ) , the following identity can be used 


1 

a + 2b cos 8 + 2c cos cp 



£~ax e “2bxcos# 


e -2cxcos Cp dx 


(170) 
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Since 


e -*cos</« _ J 0 ( iz ) + 2 


i k J k (iz) cos k\p 


(171) 


k=i 


inserting (170) and (171) in eq. (169) we obtain: 


A pq = 4 i S po 8 q0 | e_aX V 2Mx) V 2cix ^ dx 


r 
■-‘•r 


+ 2(1 - S q0 ) S p0 i* | e~ a * J 0 (2bix) J q (2cix) dx 


f 00 

+ 2(l-S p0 )S q0 iP e-“J p (2bix)J 0 (2cix)dx + 4(l-S po )(l + S q0 )iP + <J x (172) 

•'O 

r< 


x | e~ ax J p (2bix) J q (2cix) dx| 


According to Erdelyi [19] , the Laplace transforms in eq. (172) can be evaluated 
explicitly 


f 


e -ax J p (2bix) J q (2cix) dx = 


= ( P+Q \ iP + «J — C<1 3. (p + q + 1; P + 1/2; q + 1/2, 2p + 1; 

P / (a + 2 |b| + 2 |c| ) p+q+1 


2q + 1 ; 


4 |b| 


4 c 


a + 2 |bj + 2 | c | a + 2|b| + 2 |c| 


(173) 


where 3 2 is a hypergeometric function of two variables [15] . It can be shown 
that the expression for A 00 reduces to a complete elliptic integral of the first 
kind, as expected from the case d 4 0. 
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2. Exp (tA) 


For the asymmetric case the pertinent sums are of the form: 


MN 


M-l ^ N ~1 ^ f 2ttt 2 W'C 2 wr 2 ff ' t ') 

-L exp |ts 0 + t Sl e‘ M + ts 2 e* N + t s 3 e‘ M e N J 

r =0 ^= 0 


27Tirm 2vi^n 


x e M e N 


(174) 


and no constraints need be placed on the s.’s. 


The Fourier expansion is straightforward: 


e ts ° tp +< i s p y - 

1 2 4* j - CP — j 

J - o 


Vq = ^ 


q ^ P 


) ■ (q - j ) ! V s 12 , 


(175) 


ts„ 


+ 4 1 

tP + q S P S« ) 

LtJ j • (P - j ) 

J-0 


(q - j)! \ S i S 2 j 


q > p 


These expressions simplify for s = 0, to 


A - e 

pq 


„ t p+q s p s q 

ts Q 1 s ! s 2 


p! q! 


(176) 


and the summations involved in the evaluation of S UM can be performed for p and 

MN 

q separately. 


The symmetric case can be treated in terms of simple functions only for 


\( 6 , <p) = a + 2b cos 0 + 2c cos cp 
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t 


Then 


e t\(0,cp) = e ta e 2bt cos e e 2ctcosC iP 

and we make use of the expansion (171), such that 


(178) 


A pq = 4e ta I p (2bt) I q (2ct); p,q > 0 (179) 

As in the one-dimensional case, one can extend the above results to functions likp 
^"(tA), cosh^ t A) and c A (with c a constant), without difficulty. 

4. The Frequency Distribution Function Q (u?) 

The sum to be evaluated here is 


11-1 N-l 

Q(- 2 )=^ Y (180) 

r =0 4=0 

Hence we need the Fourier expansion of the generalized function S {a? - k(6, 9) } . 
Since for the asymmetric case the eigenvalues k A can be complex, the appropriate 
eigenvalue distribution function Q(a>) would have to be defined over the complex 
plane. This is not done here and instead we treat only the real, symmetric case. 

The frequency generating function k ( 8 , 9) here is taken to be 

k(6, 9) = a + 2b cos 0 + 2c cos 9 + 4d cos 9 cos 9 (181) 

with the restrictions of eq. (165) . 

In Appendix C the Fourier expansion of the generalized function S {a? - k(9, 9)} 
is effected. Here we present only the first (dominant) term: 
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* 


§{a-6j 2 + 2b cos (9 + 2c cos cp + 4d cos <9 cos cp}^ 


~F= K (/ 

/vT^v W 


(1 - v) (1 + uV 
2 (u - v) 


a — 2 |b| — 2 | c | — 4 |d| < oj 2 < a - 2 |b| + 2 |c| + 4|d| 


:J- O 

4„ 2 Vc 2 - 




K 


'O 

- 


2(u - v) 


4d 2 


v) (1 + u) \y (1 - v) (1 + u) ; 
a - 2 |b| + 2 |c| + 4 |d| < w 2 < a + 2 ]b| - 2 |c| + 4|d| 



K 


( /( 1 - v) (1 + u)\ 

\y 2(u - v) y 


|^a + 2 |b| — 2 | c | + 4 |d| < ai 2 < a + 2 |b| + 2 |c| - 4 |d| 


where K(k) is the complete elliptic integral of the first kind with modulus k(k < 1), 
and u , v are given by 


a + 2 |b| - a? a - 2 |b| - a 2 

2 ( | c | - 2 |d| ) ’ 2 C I c | + 2 |d| ) 


(183) 


The expression for Q(a 2 ) is identical in the first term with that of 8{a 2 - k( 0, <p) } , 
as can be seen from eq. (12). This result in eq. (182) is more general than that 
quoted by Montroll et al. [ 11] in two respects: 1. The addition of the mixing 

term 4d cos 6 cos <p; 2. The constraints imposed on the parameters a, b, c andd 
are weaker than those ordinarily used. Moreover this evaluation exhibits the 
advantages of our method, since alternative integral representations for the 8-function 
lead, when mixing terms are included, to integrals which are not recognizable from 
existing tables. 
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ii. Rigid Boundary Conditions 


As in the one-dimensional case, the formulas of this section differ from those 
of the previous one only because the sets C, }, i$ v > are modified. Again, it 
is clear that the changes will occur solely in the correction terms. 


The sums treated here arise from eigenvalues connected with two-dimensional 
continuant matrices. The eigenvalues considered are of the type 

\ ; •■■■: a T 21 1 cos — —1— . 2c cos — 4d cos — -1_ cos — ( 1 84) 

r - M + 1 N + 1 M . 1 N - 1 

and the matrix elements of F( a), on using the explicit form given in eq. (176) of 
Appendix E for the diagonalizing matrix T. will be 


. F('A')’J 


pq 


(M f 1) (N . 1 , 


V -1 W 


F ( • 


? in 


M • 1 


sin 


rc U 
M - 1 


~ : P 2 

nTT 


^2 
N- 1 


Pj- q, - 1 p 2 . q 2 1. . ... n 

Utilizing the relation 


(185) 


s in v sin v JL ; os ( x » v) - cos (x T 
9 " 


we see that our concern here is with a sum of the type 


MN 


(M r D (N ~ i 


v - 


j : = i 


urm -tn 

) cos — c o s 


M r 1 Nil 


( 186 ) 


where, on account of periodicity the integers m.n can be restricted to the ranges 


0 < m < 2M + 1 ; 0 < n < 2N + 1. For 


0 S v reduces to a trace-type sum. 


n 



Now we assume for F{A. ( 9, 9 ) } a Fourier expansion similar to that in 
eq. (154), with appropriate A on 's, and define quantities U. , V, ; 


U. = 

j 


1 

M + 1 



cos 


TIT j 

M + 1 ’ 


1 

N + 1 



cos 


77k % 

N + 1 


(187) 


The sums appearing in (187) have been evaluated in terms of Kronecker S ’s in eq. 
(128). Note also that U_. = U. and V_ k = V k . Then insertion of the Fourier expan- 
sion of F( A.^) and use of the U. , V k leads to the following form for S MN 


s = — y y a e (u +u + )(v_+v.) 

MN 4 / / pq pq v p“m p+m q n q+n 


p- 0 q-0 


(188) 


where 


r 


- p - q = 0 


e pq = < j P = 0, q > 0; p > 0, q = 0 


1 p > 0, q > 0 


(189) 


If we further define 


00 00 

b = y 1 y .. _ ~ + . + 

r.s £ , £_ pq pq P+r q +s 

p=0 q-0 


A e U x V 


(190) 


then S,„, can be written as follows 

MN 


S un =— + B + B + B } 

MN 4 m , —n ~m , n m , “n m,n 


(191) 
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Since the eigenvalue-generating function K(6 ,cp) considered here is identical 
with that of the symmetric case in the previous section, the Fourier coefficients 
for the various functions F(z) coincide also. Hence there is no need to repeat 
the results, and no particular application will be made to the sums S^. 

3.3 d-Dimensional Sums 

This section is devoted to sums arising from problems in three or higher 
dimensions. The sums will be over analytic functions of eigenvalues of d -dimensional 
matrices. These matrices are straight-forward generalizations of the matrices 
dealt with previously. It must be borne in mind that if one disregards the d- 
dimensional partitioning of a generalized matrix, one is left with an ordinary 
matrix the elements of which can be specified by the usual two indices. On the 
otherhand if one has to keep track of the particular submatrices which contain 
the element in question, the two simple indices p,q are to be replaced by two 
d-dimensional vectors p ,q the components of which take values from d appropriate 
ranges. Examples of these matrices, apart from the (lattice) dynamical matrices 
mentioned, are the matric representations of quantum-mechanical Hamiltonians, 
where the elements are specified by several quantum numbers. This shows also 
that the components p. of the vector indices are sometimes taken from a set 
of N. numbers not necessarily integers. 

For the remainder of this section only the three-dimensional case will be 
considered. The elements of an analytic function of A are given here by 


[F(A)] 


p.q 



k 


F <v Tpk v* 


(192) 
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where A k are generalized eigenvalues of a, T the (generalized) matrix which 
brings A to diagonal form, and p = (p x , p 2 , p 3 ); q = (q x , q 2 , q 3 ), with , q } 
, . . . , r< j) } . In the following the r (i) will be integers. 

i. Periodic Boundary Conditions 

Here we shall treat analytically only the simplest cases, since even for 
these closed form results are generally not available. 


The general asymmetric case leads to eigenvalues of the form 




N 1 N 2 N 3 

w = Z Z Z s - 

j=0 k=0 l =0 


2tt j p 2rrkq 2'n'^r 

1 HnT" 1 “n7“ 1_ n7" 

^ e 1 e 1 e J 


P = o, • • • , Nj - 1, q = 0 N 2 - 1; r = 0 N 3 - 1 


( 193 ) 


J 


Then using the form of T given in eq. ( ) of Appendix , we obtain 


Nj-l N 2 -l N 3 ~l 
X -1 


ivk Z Z Z 


F (k ) e 

pqr 


27rp(u 1 “v 1 ) 2rrq(u 2 *v 2 ) 27Tr ( u 3~ V 3 ) 

‘ ~ 1 n7 


N 1 N 2 
1 e 1 x e 


p-0 q-0 r = 0 


( 194 ) 


Therefore we are concerned with sums of the type 


Nj-i n 2 -i n 3 -i 


27Tipn 1 277'iqn 2 27rirn 3 


5n i n 2 N 3 = NjN N ZIE F( ^ pqr) e 


N 0 N., 

e e J 


( 195 ) 


p=0 q-0 r=0 


where n x , n 2 , n 3 are integers which can be restricted, because of periodicity, to 
the ranges: 
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o < nj < Nj - 1 ; 0 < n 2 < N 2 - 1; 0 < n 3 < N 3 - 1 . 

Proceeding as before, we define an eigenvalue generating function <p ,\p) 
and expand F{\( 6, cp,0)} in a triple Fourier series. The final result is 


Here 


N 1 N 2 N 3 


A* 


n l n 2 n 3 


(195) 


00 CD CD 

Ap l p 2 p 3 ~ L Ap ! + j N r P 2 +kN 2 . P 3 + ^N 3 

j =0 k=0 ^= 0 


(197) 


and 


A jk ^ - — - - JJ J F{A.(<9, tp, 1 //)} e -i ^ j <9+k SP dddydxp (198) 

We have assumed k (0, <p,0) to contain only positive powers of e i6 , e iC P , e 1 ^ and 
therefore the integers j , k , •£ >0. 

The general symmetric case leads to eigenvalues of the form: 


k 

pqr 


N r 1 


N 2 -l 


NS" 1 


L LL 

j-0 k=0 l =0 


2ttp'i 277 qk 
k £ cos -- cos 


N 


x 


N. 


cos 


277 T't 

K 


(199) 


where the s. k ^ ’s satisfy certain identities, similar to those imposed in two 
dimensions. Then 


F{\(0, cp, 0)} = 


CD 00 00 

m: 

p l =0 p 2 =0 P 3 =0 


A 


p l p 2 p 3 


£p l p 2 p 3 


cos Pj 6 cos p 2 9 cos p 3 \p 


( 200 ) 
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where 


8 Pi + P 2 + P 3 = 0 


4 Pi + P 2 + P 3 = 1 


j Pi + P 2 + P 3 = 2 


1 Pi + P 2 + P 3 ~ 3 


r 2 tt r 

: — — jj J F{A.(< 9 , cp, \jj)} cos p x & cos p 2 cp cos p 3 0 d 0 dcpd 0 ( 202 ) 


Hence a sum of the type 


Nj-l N 2 -l N 3 -l 


N 1 N 2 N 3 N x N 2 N 3 




"* p n i q n 2 rn 3 1 

_ N 1 N 2 N 3 J 


p=0 q-0 r-0 


with n , n 2 , n 3 as in eq. (195), can be rewritten as follows: 


N 1 N 2 N 3 8 n l n 2 n 3 


and the A* 's are given by 


CO QO CO 


Apqr L LL £ jk^ ^p+jNj . q + kN 2 , r+'CNj 


j=0 k = 0 l =0 
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When the exponentials in the sum (203) are replaced by the respective cosine 
terms the modified sum can be written as 


+1 

s - 1 y a* 

N.N.N, ~ ~Q~ /. 1-“ „ l-£ l-y 

112 8 a "l + — N l^n 2 + — N 2-^3 + — N 3 


(206) 


1 . The Inverse A -1 

Here we shall treat only the symmetric case for the simplest frequency 
generating function k(8, cp , \p ): 

k(9, cp, \p) = a + 2b cos 6 + 2c cos cp + 2d cos \p (207) 

The corresponding Fourier coefficients are 


A 

pqr 




cos kd cos qqp cos ri p d <9dcpdv// 
a+2bcos^ + 2ccostp + 2d cost// 


(208) 


The only exact results available are for A Q00 , and this only for particular 
relative values of a , b , c and d. For a = - 6b and b = c = d <0, the integral has 
been evaluated by Watson [20] , with the result 


ooo 


3 ff 2 7 d 6 dtp di/> 

arj 3 JJ J 3 - cos 6 - cos cp - cos i p 


96 


7T 2 a 


[18 + 12/2 - 10/3 - 7 /6 x 


x K 2 [(2 -/3) (/ 3 -t/2)] 


( 209 ) 


where K (k) is the complete elliptic integral of the first kind. 
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Montroll [21] has generalized this result to the case d,b = c < 0 and a = 
2 (2 |b| + |d| ), such that 


with 


000 


b sn 


b 77' 


dd dcp d\p 

(2 + a 2 ) - cos 8 - cos - a 2 cos 0 


32 


i/|bdj 


[SyTl - v^Tl] K(kj) K(k 2 ) 


77 


( 210 ) 




y = 


4 + 3a 2 


a 


kj = y [v/y - 1 - /y - 3] [/y + 1 - v/y - 1 ] 


k 2 = y [v/y - 1 + v/y - 3] [v/y + 1 - v/y - 1 ] 


■y 






( 211 ) 


andK(k) the first complete elliptic integral as before. 

There exists a formal expression for A pq r in terms of a hypergeometric 
function of three variables. The result arises as follows: 


A 

pqr 


1 



e "aT e “T( 2b cos (9 + 2c cosCp + 2d cos </>) 


x cos p/9 cos qcp cos n p d 6 dcp d0 


( 212 ) 
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where we assume He (a) > 0 . Since for integral n 



we obtain 


r 2n 

I g— /? T co s x cos nx c ( x _ 277 i n J q (/3Ti) 


CD 

A pqr = 8i p+q+r J e" aT J p (2bri) J q (2cri) J r (2dri) dt 


(213) 


(214) 


The integral in eq. (214) is a Laplace transform. Its value as given by [19] leads 
to 


8(-l) p+q+r W 3 c q d r (p+q + r)! 
Pqr (a - 2b - 2c - 2d) p+q + r+1 P ! r! 


(215) 


Jp + ^ + r + 1 , P + , q + „ . r + _ > 2p + 1 , 2 q + 1 , 2 r + 1 , — 


a -2b -2c -2d 


4c 4d 

a - 2b - 2c - 2d ’ a _ 2b - 2c - 2d 


The function 3 3 , known as the Lauricella function, is a three-dimensional power 
series. Its properties have not been thoroughly investigated and its alternative 
representations, if any, are not known. 


Tables for the general integral A pqr are available for certain ranges of the 

parameters which appear in its denominator [22] . The same reference gives 

also asymptotic forms of A for large values of the indices. 

pqr 

Returning to A Q00 it can be shown that two integrations can be performed 
exactly in terms of complete elliptic integrals on using the results of Appendix B 
for the two-dimensional case, and this can be done even when k(d, 9 , 1 p) contains 
mixed cosine terms. But so far, attempts to carry out the third integration 
analytically have not been successful. 
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2. Exp (tA) 


For the asymmetric case we treat the sums 


N r l N 2 -l N 3 -l 


\ n 2 n 3 N,N 0 N, 2_j L—> 2-J 

123 j =0 k = 0 1=0 

277 injj 27 rin 2 k 277 in 3 ''fc 

N. N " n 7 ~ 

x e 1 e 2 e 3 


exp 


t 2nj _ 2wk . 2tt ''C [ 

i_ N7 i_ N7 i N7 r 

ts 0 + tSj e + ts 2 e + ts 3 e 3 J 


(216) 


Since there is no mixing term, the Fourier coefficients can be immediately written 
down 


ts 0 rr 2lT r 

A =- II r e " i <P £,+qC P + r ' / '> exp (ts x e i0 + ts 2 e iC P + ts 3 e 1 ^} d0 dcp d 0 

Pqr (2rr) 3 JJ 0 J 


e tS ° tP +q+r sP s q s' 
p! q! r! 


(217) 


The summations involved in the evaluation of the sum 
separately for p,q and r. 


j NjN 2 N 3 


can be performed 


The symmetric case leads to the Fourier coefficients 


A 


pqr 



exp {2bt cos 6 + 2ct cos cp + 2dt cos 0} 


x cos p 9 cos qcp cos rt p dd dtp d\p 


= 8e ta I p (2bt) I q (2ct) I r (2dt) 
where I v (z) is the modified Bessel function of the first kind. 


(218) 
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3. Frequency Distribution Function 


The remarks concerning the evaluation of the Fourier coefficients related 
to A" 1 are pertinent also here. The triple integrals that are to be evaluated 
in order to find the Fourier expansion of the generalized function. 

S(a? - a - 2b cos 6 - 2c cos <p - 2d cos \ p) 
present the same difficulties. 

Therefore if Q( a? ) is needed only for the conversion of finite triple sums into 
integrals, this approach will not be practical for numerical work because of the 
free parameter a? . On the other hand, our method is still useful. 
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CHAPTER IV 


LATTICES WITH DEFECTS AND PERTURBATION THEORY 

Defects and imperfections in crystal lattices have been the subject of numerous 
and relatively recent investigations [11] . The most widely used method for the 
treatment of such lattices is Montroll's method of Green's functions, developed 
and elaborated in collaboration with Potts, Maradudin and Weiss [23] . 

An outline of this method is presented below in connection with the evaluation 
of the frequencies of an arbitrary imperfect lattice. Later on the discussion will 
be specialized to the case of one isolated mass defect, so as to compare the re- 
sults with those obtained from perturbation theory. 

Let us assume that the dynamical matrix associated with an imperfect lattice 
is 

A = A q -A' (1) 

where A 0 is the dynamical matrix connected with the perfect lattice, while A' con- 
tains the deviations of the imperfect from the perfect lattice. The characteristic 
equations can be written as usual 

(A - A.I) u = 0 or (A 0 -\I - A') u = 0 (2) 

in which u is an N-dimensional column vector. 

Eq. (2) can be rewritten as follows: 

u = GA' u 
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G is the Green's matrix of A 0 , 


G = (A 0 - KIT 1 ( 4 > 

defined for all k not on the spectrum of A 0 . Eq. (3) can be useful in solving for 
uby the method of successive approximations. More important is the fact that if 
only a few of the A' -elements are non-vanishing, then one can find the character- 
istic equation explicitly in terms of the Green's matrix elements fromeq. (3). 

It is essential therefore to possess explicit expressions for the elements g r . of 
G • In eq. (22) of Appendix D these elements are given by 


IN 

L 

k=l 


T T" 1 

rk kj 

\<°> - A. 

k 


(5) 


in which and T pq denote, respectively, the eigenvalues of A Q and the elements 

of its diagonalizing matrix T . For generalized matrices the scalar indices will 
be replaced by vectors of appropriate dimensionality. 

For the case of one-dimensional systems with simple a' perturbations, 
the characteristic equation can be written down directly without using Green's 
functions. For higher dimensional lattices, the method seems to be indispensable 
for finding the characteristic equation of the perturbed lattice, even though the 
exact forms of the Green's functions are not known, i.e., the sums in eq. (5) 
cannot be evaluated in closed form. 

It is clear that once the characteristic equation is known, the mode of its 
solution is independent of the Green's function method. The same applies to the 
evaluation of the eigenvectors u of eq. (2). 
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In the rest of this chapter we apply second order perturbation theory to a 
monatomic linear chain with one isotopic impurity, symmetrically situated, and 
rigid boundary conditions. 

We assume therefore 



( 6 ) 


with 


a = 2 a - Mw 2 ; 
a* = - e M oj 2 ; 


b = - a 



(7) 


The eigenvalues of A 0 are 

= a + 2b cos — -~r ; k = 1 , . . . , 2n + 1 (8) 

* 2n + 2 

In the remainder of this chapter we shall use the notation Q = n/2n + 1. To 
apply perturbation theory we first bring A 0 to diagonal form by performing the 
similarity transformation T given in eq. (12) of Appendix E. Then the matrix 
to be treated by perturbation is 


with 


H = H 0 + H' 


(H.) . = \<°) S . 

v O'rk k rk 



n + 1 


TT't 

sin sm 

2 


77J 

2 


(9) 


(10) 
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Since the eigenvalues are simple we can use ordinary perturbation theory. 

Then to second order in a' the eigenvalues \ will be given by: 


\ = a + 2b cos ki9 

k n + 1 



+ 





2b (cosk# - cos t6) 


(U) 


It is clear immediately from the form (10) of H'^. and eq. (11) that for k = 
2j; j = 1, . . . , n, we obtain to all orders 


k = k( 0) = a + 2b cos — - 

n + 1 


j = i. 


n 


( 12 ) 


Note that the eigenfrequencies in eq. (12) correspond also to a monatomic linear 
chain of n particles with rigid boundary conditions. The result that n frequencies 
are not perturbed by an isotopic mass defect situated at the center of symmetry 
can be arrived at also by purely matric methods. 

Returning to eq. (11) and putting k = 2 j — 1 ; j = 1, . . . , n + 1 we can 
write 


k. 


2 j -1 


= a+ 2bcos(2j-l)0 — 

n + 1 




L 

r-1 

r/i 


1 

cos(2j -1) 0-cos (2r -1) 6 


The restricted sum in eq. (13) can be rewritten as follows 
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cos(2j -l)#-cos(2r-l)i9 / j cos(2j - 1) 9 - cos r9 

r = l 

r/2j-l 


- T i 

cos (2 j - 1) 0 -cos -ILL- 
n + 1 


The unrestricted sum in eq. (14) can be evaluated in closed form by using 
eq. (6) of Appendix D with 


F (k) = ± [V"" 1 - V" n_1 ] 


U = /\ 2 -l; V = \+U 


Then 


F'(V> = - — F(k) +. n + 1 [V n+1 + V" n_1 ] (1 

U 2 U 2 

Substituting cos (2j -l)0for k in F'(k)/F(k) we obtain the value of the un- 
restricted sum 


£ 1 — 

r =i cos (2 j - 1) 0-cos 


_ cos (2 j -1) 0 
~ sin 2 (2j - 1)0 


To evaluate the restricted sum on the r.h.s. of eq. (14) we have to use 
formula (16) of Appendix D, this time with 
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F(k) = A- [V 2n+2 - V' 2n " 2 ] 


(18) 


where U and V are as defined in eq. (15). For this function we have 


F* (k) = - — F (k) + — + - [V 2n+2 + V' 2 "" 2 ] 

u 2 u 2 


(19) 


F" (k) = 


(2n + 2) 2 -1 2\ 2 

U 2 U 4 J 


F(\) - A. F'(\) - 4(n + --— [V 2n+2 - V‘ 2n " 2 ] 


U 2 


u 4 


Substituting cos (2 j -1)6 for k in F"(\)/F'(M w © finally obtain 


2n+l 

I cos (2j - 1) 6 - cos t6 

r-l 

r/2j‘l 


3 cos (2 j - 1 ) 6 
2 sin 2 (2 j -1)6 


Collecting results we can write for 


( 20 ) 


k , = a + 2b cos (2 j - 1) 6 — + A- 

2 J 1 n + 1 4b 



cos (2j - 1) ^ 
sin 2 (2 j - 1 ) 6 


j = 1, . . . , n + 1 ; 6 - 


77 


2n + 2 


( 21 ) 


It can be shown that the contribution to \ ^ from third order vanishes identically 

and it is probably true that all odd order contributions vanish. 

We proceed now to the evaluation of the eigenvectors to the same order of 
approximation. These are linear combinations of the eigenvectors for the 

unperturbed lattice. In our case the components of are 
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* « 


u ^° r ) = — 1 sin - 7 —^ ; r,^ = l 2n + l 


Vn + 1 


2n + 2 


(22) 


Then if is the eigenvector associated with of eq. (11) and its 
components , we obtain on using eq. (22) and second order perturbation theory: 


r 




2n+i sin— sinr^d? 


< 


n + 


sink^0 - 


2b (n + 1) 


7rk \ 1 

sm — > 

2 


cos k<9 - cos r 0 


r=l 

r/k 


1 ( a' \ 2 . 77 k 

( ] sm 

4b 2 \ n + V 2 


2n+i sin— sinp^<9 

K ' 2 


z 

p=i 

L_p/k 


cos k 0 - cos p<9 


2n+l 

Z 

r*l 

_r/k 


(sin — ) 

A 2 ; 


2 — 1 


cos k£ — cos rd 


-i (A) H) £ 


2n + l sin — sin r'C# 
3 * — . 2 


— J 1 (cos k<9 - cos r 0) 2 

r/k 


2 / 2 2n+1 /sinUI^ 

Y — i 2Z_ 

8b 2 \ n + 1/ V 2 / (cos k£? - cos r<9) 2 


r/k 


f 


(9 = • 


2n +2 


; 4, k = 1 , . . , 


(23) 


All of the sums in eq. (23) can be evaluated exactly when use is made of the 
results in Appendix D. We denote these sums by S , m = 1, 2, . . . , 5 according 
to their order of appearance in eq. (23). It is also seen that Sj and S 2 are 
identical. 


From eq. (23) we immediately find that for k = 2 j , j = 1, . . . , n, u 2 . coincides 
with • Again this result can be shown to be exact by using matric methods. 


2n + 1 
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* 

We consider now the sums S m for k = 2 j -1, j = 1 n+1. Then the 

sum S 3 is identical with the one exhibited on the l.h.s. of eq. (14), and we have 


g _ 1 cos (2j - 1) <9 
^ sin 2 (2j - 1) 6 

The sum S s can be written as follows: 

S s =S s - S s 


(24) 


(25) 


where 


2n + l 
r*l 

r/2j~l 

(26) 

The sum S 3 can be easily evaluated on using formula (20) of Appendix D with 
F(V) given in eq. (18) above: 


(cos (2j -1) 6 - cos r (9') : 




r = l [COs(2j-l)5-COS 


7rr^? 


16 (n + l) 2 sin 2 (2 j - 1) 6 - 17 cos 2 ( 2 j - 1) 6 - 16 
5 ” “ 12 sin 4 (2 j - 1) ~6 


(27) 


Similarly the sum S" can be evaluated by using formula (7) of the same 
Appendix with F(\.) as given in eq. (15). The result is 


(n + l) 2 sin 2 (2 j - 1) 6 - 1 
sin 4 (2 j -1)0 ” 


(28) 


Finally 


S 5 = 


4(n + l) 2 sin 2 (2 j - 1) 9 - 17 cos 2 (2 j - 1) 0 - 4 

12 sin 4 77(2 Lli] 

2n + 2 


(29) 
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The sum Sj can be written as 


s . is; + sp 


where 


All » x 

\ ' cos mrg 

/ ■ cos (2 j - 1) 9 - 


cos r# 


r/2j-l 


m = n + 1 - 'i 


m = n + 1 + 'C 


For 1 < m < 2n + 2 the sum in eq. (31) can be evaluated by using formula (15) 
of Appendix D, with 


£ ^ y 2n^2^m ^ y *"* 2 n 2 


F(M = U [V _2n+2 - V 2n ~ 2 ] 


The result is 


\ 1 cos mr<9 _ 1 

ZL cos (2j - 1) 9 - cos Td ~ 2 sin 2 C2j _ D ~e 

r/2j-l 


{l _ (-1)”+ [1 + (-!)■] cos(2 j -1)( 


+ cos(2j - l)m0cos(2j -1) d -2(2n + 2 -m) sin (2 j - l)mi9 sin (2 j -1) 


If we assume < n + 1 , then both values of m in eq. (31) satisfy the restric- 
tion 1 < m < 2n + 2 for which eq. (33) is valid and we can write for Sj : 
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s. = 


1 2 sin 2 (2 j - 1) 6 


1 _(_l)n + i-^ + [ 1 + (-l)n + i-^] cos (2 j -1)6 


- 2(n + 1) sin ^ sin^j-D^cos^j-l.)^#^; -€■ <n + 1 (34) 


For -i > n + 1 , the sum S" becomes 


s; 


2n+l 

L 

r=l 

r/2j-l 


(-l) r cos pr0 
cos(2j -1) 6 - cos r 6 


(35) 


where 


p = -l - n - 1 . 


This sum can be readily evaluated by using formula (15) of Appendix D with 


f (M = VP + V'p 




(36) 


F (k) = U [V 2n+2 - V' 2n " 2 ] 


The result is 


S i' = 


1 


2 sin 2 (2 j - 1) 6 


{l _ (_ 1)P + [1 + (- 1)p] cos (2 j - 1)6 


- 2p sin (2j - 1) p6 sin (2j - 1) 0 - cos (2 j - 1) p#cos (2 j - 1) @1) 


Finally combining eq. (37) with the result in eq. (33) (for m = -C-n-1) we obtain 
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s. = 


2 sin 2 (2 j -1)0 


|l - (- l/'"' 1 + [!+(- l/ - "' 1 ] cos (2 j -1)0 


+ 2(n + 1) sin X) sin (2j - 1) 6 cos (2 j - 1) Ig). ; 'C > n + 1 (38) 


?j- ; •£ > n 


The last sum to be evaluated is S 4 , which again can be decomposed into two sums 


s 4 =-| fs; + s;] 


(39) 


where 


S' | 2n+l 


L 


cos mrg 


m = n + 1 - -t 


t—rf (cos (2j - 1) 9 - cos r0) 2 
S" i. r/2j-i m = n + 1 + 


(40) 


The sum in eq. (40) can be written down by using formula (19) of the Appendix 


with 


and 


r / \ \ T» 2n^"2 m «T”2n*"2^in 

f ( M = V + V 


F(M = U[V 2n+2 - V _2n " 2 ] 




for 1 < m < 2n + 2 


(41) 


f(M = V m + V m 


for 2n+2<m<4n+4 


F(\) = U[V 2n+2 - V 2n_2 ] 
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The results are: 


cos mr0 


2n+l 

Z-J [cos(2j -1) 0 - cos r0] 2 2 sin 4 (2j - 1)6 

r/2j-l 


2 [(- l) m - 1] cos(2j - 1)6 


- [1 + (- l) m ] [1 + cos 2 (2 j - 1) 6 } + (2n + 2 -m) 2 cos m(2j - 1) 9 sin 2 (2 j - 1) 6 


+ 2 (2n + 2 -m) sinm(2j -1) 0 cos (2 j -1) 6 sin(2j - 1) 6 


-A cos m (2 j - 1) 6 cos 2 (2j -1 ) 9 - COsm(2j ~ - [4 sin 2 (2j - 1) 6 

2 3 sin 2 (2 j - 1) 0 


+ 15 cos 2 (2 j -l)0-4(n + l) 2 sin 2 (2j -1)0] 


; 1 < m < 2n + 2 


(43) 


and 


(- l) r cos pr0 


2n + l 

y - 

[cos (2j -l)0-cosr0] 2 2 sin 4 (2j -1)0 

r/2j-l 


2 [(- l) p - 1] cos (2 j -1) 0 


_ [(_1)P + 1] [1 +COS 2 (2j -1) 0] +2p cos (2 j -1) 0sin (2 j -1) 0sinp(2j -1)0 
-p 2 cosp(2j - l)0sin 2 (2j -l)0+icosp(2j -1) 0cos 2 (2j -1)0 

+ C ° Sp(2j ~ X 1 - [4 sin 2 (2 j - 1) 0 + 15 cos 2 (2j - 1) 0 

3 sin 2 (2 j - 1)0 

- 4(n + l) 2 sin 2 (2 j - 1) 9 ] > (44) 

where ' 


1 <p=m-2n-2< 2n and m > 2n + 2 . 


On using these results with the appropriate values of m from eq. (40), and 


inserting in eq. (39), we obtain 


S. = -1 |2 [(-l) n+1 "^-l] cos(2j -1)0- [(-l) n+1 '^ + l] [1 +cos 2 (2j -1)0] 

2 sin (2 j - 1)0 l 


+ 2^(n + 1) sin -Hi 2 ! - sin 2 (2j -1) 0sin ^(2j - 1) 0 


+ 2(n + 1) sinZI^i — 1} sin (2 j -l)0cos (2j -l)0sin^(2j -1)0 


; 1 <'l< n 


and 


S 4 = 


+ 1 

(45) 


2 sin 4 (2 j -1)0 


|2 [(-l)^-" -1 -1] cos (2j -1) 0- [(-l)^"- 1 +1] [1 +cos 2 (2j -1)0] 


+ 2(n + 1) (2n + 2 - l) sin 


. 7T(2j - 1) . 2 /o; 


sin 2 (2j - 1) 0 sin ^(2j - 1) 0 


+ 2(n + 1) sin L — L} sin (2 j — 1) 0cos (2 j - 1)0 cos ^(2j -l)0l-;n + l<4<2n + l 

J (46) 

Returning to the perturbed eigenvalues of eq. (21), we remark that using the 
techniques developed here one can obtain 


2n + 1 


\ = (2n + 1) a - a' = 3r (A Q - A') 

k = l 


(47) 


where A. are as in eq. (12) and \ 2j _j as in eq. (21). This shows that in this 
order of approximation the sum of the approximate eigenvalues is equal to the sum 
of the exact eigenvalues of A 0 - A' . 

The approximate perturbed frequencies of the linear lattice are obtained from 
the equations 
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= 0; j = 1 , • • • , n + 1 


(48) 


\ 


2 j -1 


with X. 2 j from eq. (21). Ail of these are quadratic equations in a? except for 
the particular case n = 2r , s = 1, 2, . . . , which leads for 2j - 1 = n + 1 to the 
particular frequency 


CO 


2 


2a 

M 



(49) 


In all other cases we solve for co 2 to second order in e , choosing those solutions 
of the quadratic equations resulting in real co' s: 


CO 


2 

2j“l 


4a 

M 


sin 4 


rrj 2j - 1) 
4(n + 1) 


r 


<0 + 


€ 

n + 1 




+ 


, 1 

1 + -COS 


1 + cos 


77(2] -1) 
2n + 2 
tt( 2 j - 1) 
2n + 2 


■y 



j 


(50) 


Eq. (50) shows that the perturbed frequencies will be decreased or increased with 
regard to the unperturbed ones, according to whether e < 0 (a heavier impurity 
mass) or e > 0 (a lighter impurity mass). This is in agreement with general 
theorems of Rayleigh [24] concerning the effect of additional constraints on 
arbitrary vibrating systems. 


It is of interest to note that if we put 4a/M = co 2 (a> L is then the top of the 
frequency band) and take e > 0, one of the frequencies in eq. (50) will emerge 
above co L as n becomes large. This can be seen by considering a>^ n+i : 


2 2 2 
W 2n + 1 =aJ L COS 


V 


4 (n + 1 ) 


i 1 77 

1 cos 

, e 2 2n + 2 

1 + : + 


n + 1 


1 - cos 


77 


n + 1 


(51) 


2n + 2 


J 
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Proceeding to the limit n oo in eq. (51) we find 


" Id h lim 

n — oo 


CO 


2n + l 


= cot 



(52) 


This is the frequency of the "localized” mode, and has been obtained here 
directly from perturbation theory, in contradiction to statements that this is 
impossible [11] . Montroll & al. have shown [8] that the exact localized mode 
frequency can be obtained by proceeding to the same limit from the original 
eigenvalue equation. Their result is 


CO 


2 _ 
LO ~ 


2 - e 2 


(53) 


The discrepancy between eqs. (52) and (53) to second order in e is due to a 
difference in boundary conditions. 
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appendix a 


EQUATIONS OF MOTION FOR TWO- AND THREE-DIMENSIONAL 

RECTANGULAR LATTICES 

1. Two-Dimensional Lattices 
1.1 Free Boundaries 
Let us denote by 



the displacement of the particle located at the site la 1 + m a 2 . Then we can write, 
for particles not situated on the boundary, the equations 


M 


d 2 u 


4.. 


dt 2 


a 0 


0 a' 


W + i. m ~ 2 n, m + 


f P o\ 
o yS / 


W.mn ~ 2 u l. m + 


7 S 

+ i ) { n + i, m+ i - 2 n, m + 

S y' 


7 -S\ 

+ 1 W-1..+1 - 2 u l, m + u ^n, m -i } 

,-S W 


( 2 ) 
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m 


which -£ = 2, . . . , Nj - 1; m = 2 N 2 - 1, and where 


a = a ; a' = a' ; j3 = a, ; /3' 




y - a 3 + — ^ 3 ~ “3 

5* I O'* 


3 I + a 2 


a, a, 

8 =— ! i-1- (a 3 -ap 


( a * “ a 0’ T* - a 3 +— — ^ a 3 “ a 3^ 

2 i o ^ 


a i + a 2 


> 


( 3 ) 


3 1 + a 2 


If we write 


U/ t ,m = ,m e i&)t 

the eqs. (2) can be put in the reduced form 


(4) 






(5) 


where 



In writing down the eqs. (2) use was made of eqs. (60) and (61) for finding the 
appropriate force density tensors A, B E. 

One can write down the equations of motion for particles on the boundary in 
a similar fashion. These will differ from eqs. (2) since a smaller number of 
neighbors is involved. Here we shall only incorporate results in the dynamical 
matrix. 
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It is not difficult to find the dynamical matrix if one adopts the particle 
labeling shown in Fig. 3 of Chapter I. Using this one-dimensional labeling and 
the reduced equations of motion, one obtains the dynamical matrix exhibited in 
eq. (62) of Chapter I, with 

Aj =A+B+D+E A 3 = A + B + C + D + 2E 

A 2 =A + C + D + E A 4 =A + B + C + 2D + E (7) 

and A , B , . . . , E as in eq. (6) . 

1.2 Rigid Boundaries 

With the same notation as above, the equations of motion for interior particles 
are as in eq. (5) . The equations for marginal particles are obtained by using 
eq. (2) and the observation that a rigid wall, by definition, cannot move; i.e., its 
displacement vector u is zero. The changes obtained when the appropriate u's 
are suppressed, have been incorporated in the dynamical matrix for a monatomic 
lattice as shown in eq. (78) of Chapter I. 

1.3 Periodic Boundaries 

Again for interior particles the reduced equations are those shown in eq. (5) . 

The periodicity of the boundary conditions imply that there are no marginal particles. 
For instance, the particle labeled 1 in Fig. 3 of Chapter I has as immediate neigh- 
bors the eight particles labeled there 2, Nj +1, N, + 2, N , 2N , (N 2 - 1) Nj + 1, 

(N 2 - 1) Nj + 2, and . The resulting dynamical matrix for a monatomic 
lattice is shown in eq. (93) of Chapter I. 

2. Three-Dimensional Lattices 

The unit vectors e , § , n and the force constants associated with the 26 
immediate neighbors of the model described in §1.3 are given below in three groups. 
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1. First Neighbors 



e , a 

§, a' 

n 

& 

{k, 


(i, -i}— a; 

o* > j>^ a i' 

{j» 

-j>^ a 2 

(k, -k} — a' 2 

{ i . i }'•“*' a" 

(i. 

-i} *" a 3 

<3. -j>— «' s 

{k, k }-*-*- a." 

Second Neighbors 




e I fi 

§ 

fi' n fir 
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3. Third Neighbors 

The eight third-nearest neighbors are equivalent with respect to the central 
particle on which they act, and therefore only three force constants appear: y , 
y ‘ , y“ associated with t , § , n , respectively. 
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- Sj i + a 2 j - a 3 k a 2 i + a x j - a t a 3 1 + a 2 a 3 j + a * 2 k 


- a ! i - a 2 j + a 3 k 


a i i + a 2 j + a 3 k 


“12 

a 2 1 - 3 1 J 

a i2 

a 2 £ - 3 1 J 

a i 2 


3 1 2 3 


^ a 3 1 + a 2 a 3 J + a i 2 k 


a i2 3 


- a i a 3 i - a 2 a 3 J + a ?2 k 


a !2 3 


Here a 12 is as in eq. (8) and 


a = /af+a^ + aj 


(9) 


Just as in the two-dimensional case, the equations of motion for particles 
not on the boundaries are the same for all boundary conditions involved. To 
write these down, we denote by up the 3-dimensional displacement vector of 

Ay run 

the particle located at the equilibrium position <ta 1 + ma 2 + na 3 and introduce the 
time-independent vectors U , 


up = Up e 

Aymn Aymn 


i<vt 


( 10 ) 


For a monatomic lattice the reduced equations of motion read 

*■ B , «Je-i * B 2 (^.-1 „ + 

* V 11 -!™-! tU -t. n + 1 > nie - 1 , 1 m + 1 n ^+ c 2 oifc. lm+ln+ Ue +Jra - ln ) 

+ C 3 +U^ +1 mn+1 ) +C 4 mn + i +l!t+i mn _j) +C 5 n -i +k k.m+l n + 1^ 

+ ^6 ^m-1 n+1 + ^£m+l n~l ^ + ^1 m~l n“l + ^t+1 m+1 n+1 ^ + ^2 ^£"1 m+1 n"l + ^C+l m~l n + 1^ 


+ D 3%ln,+ln+l +l) 4 t l m -l n -l ) +D 4 (U -{-l m -ln + l +U l+l m +ln-l ) =0 


( 11 ) 
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The 3x3 matrices Bj , . . . , D 4 are the appropriate force density tensors 
preceded by a (-) sign and are given as follows 

- = diag (cij aj a"); - B 2 = diag (a 2 a 2 a.'); - B 3 = di ag (a' a" a 3 ) ( 12 ) 


-c = 

r = 1,2 a 


Pi 


12 


-c 


A 


r = 3,4 a 


(“ D f+1 a , a 2 


/ 


2 

2 

13 


~1 

0 


(— l ) r+1 a x a 2 O’ 

^2 Q 


V" l^ r+1 a i a 3 0 


£ 


‘2 

0 


*12 


0, 



(-l) r ai a 2 0 


+ / 8 J diag (0 0 1 ) 

0 (- l) r+I a , a 3 \ / 

0 0 


/ 


$ 


13 


\(-l) r a 1 a a 


+ # diag (0 1 0 ) 



-C = 


A 


r = 5,6 a 


23 



0 

0 

0 



0 

0 

° 1 



0 

a 2 

a 2 

(-l) r+ 1 a 2 

' 

a 3 

7 

f 

0 

a 3 

(-l) r a 2 a 3 

+ 


0 


a 3 

1 

a 23 

1° 

(-l) r a 2 a 3 

a 2 1 



/ 


-D 


_ y 


r = 1,2 & 


(-D r+ 1 a 1 


a i a 3 


(-l) r+1 ai a 2 


(-l) r+ 1 a 2 a 3 


+ / 3 "diag (10 0 ) 

a i a 3 \ 


(15) 


^ r+1 , 
.2 


(-l) r+ 1 a 2 a 3 


/ 


a 2 a 2 
a i a 3 


7 


2 2 
3 «12 


(-l) r a x a 2 a 

- 3 1 a 3 a i2 


I 


y 


/ a i 

a 2 


(_l) r+1 ai a 2 a2 


~ a i a 3 a i2 


12 


\ 


(-D'a^ 

a i 

0 


a 2 a 2 
a 2 a 3 


(-l) r a 2 a 3 a2 

(“lVa 2 a 3 a * 2 a < 2 / 


(16) 
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/ a i “ a i a 3 \ 


-D 


_ y 


r = 3,4 


(- l) r 3j a. 


(- l ) r+1 a, a 


2 3 


\ - a i a 3 


(- l) r+1 a a 


/ 


7 


1 2 


(- l) r+1 a,a 


12 


2 3 




7 


2 2 
8 a l2 


/ a 2 a 2 (-D'a^a 2 a^af, \ 

(-D'a^a 2 (-l) ra 2 a 3 a ?2 

\ a X 8 3 a i2 ^“^ ra 2 a 3 a !2 


1 2 


(-l) r+i ai a 2 


“l 

0 


(17) 


V 

o 

°/ 


The matrix A in eq. (11) is given by 


- A = 2 


L 


B 




D r + Mw 2 I 


(18) 


The number of distinct force constants in this model is 21. For a cubic unit 
cell this number reduces to 9, since by symmetry = a 2 = a 3 ; aj = a' = a 3 ; 

; A = p 2 = ^ : ^ ^ and ^ ^ 

The dynamical matrices for the several boundary conditions are obtained 
without much effort by using the eqs. of motion (11) and by adopting, instead of 
(<£, m,n), the one-dimensional labeling of particles shown in Fig. 5 of Chapter I. 
The matrices are exhibited in the appropriate sections of the same chapter. 


140 



appendix b 


EVALUATION OF THE INTEGRAL 


A 

pq 



cos p 9 cos qcp d0dcp 

a + 2b cos 9 + 2c cos cp + 4d cos 9 cos cp 


We consider first the case p = q = 0, when an exact result can be given. 
We proceed to integrate first over 9 . 

To perform this integration we define 

A = a + 2c cos <p; B = b + 2d cos cp . (1) 

Then 

a + 2b cos 6 \ 2c cos cp + 4d cos 9 cos cp = A + 2B cos 6 (2) 

and the first integral to be evaluated is 


I = 



69 

A + 2B cos 9 


( 3 ) 


The value of I is given in [16] where different parameters are used. Our 
needs here require that we give all possible cases in a different form for the 
more general integral I n : 


I 



cos n 9 
A + 2B cos 9 


69 = 
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r 


77 


/A i) 2 - 4B 2 


-A. g M- fl . | A | , 2 |B| 
l 2B sgn A / 


= < 


77 


2 /A 2 - 4B 2 


/A 2 - 4B 2 - aT / /A 2 - 4B 2 + A\" 


; |A| < 2 |B| ( 4 ) 




V 2B J \ 2B j 

(In this case the integral exists only as a principal value) 


where 


Then 


sgn A = • 


1 A > 0 


- 1 A < 0 


(5) 


77 


sgn A |A| > 2 |B| 


I = 


/A 2 - 4B 2 


( 6 ) 


othe rwi se 


We assume now that the following conditions are satisfied by b, c, and d: 

b, c, d < 0; |b| > | c | > 2|d| . (7) 

Then it is clear that B < 0, and we have to consider the two possible cases for 
all <p: i) A < 0; ii) A > 0. The case A = 0 can be obtained as a limiting value 
from either i) or ii). If A vanishes for some cp's, the same treatment can be 
given for the subintervals where the cases i) and ii) are pertinent. 


i) A < 0 

Here the conditions |A| > 2 |b| and 0 < cp < 77 imply 


- (a + 2c cos (p) > — 2(b + 2d cos cp) 


or using (7) 


cos cp > 


a + 2 |b| 
2(|c| - 2 |d|) 


( 8 ) 

(9) 
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Let us define now 


a + 2 [b| _ 

2(|c| - 2|d|) 


(10) 


The remaining integral will not vanish only if u < 1 since -1 < cos qp < 1 . This 
condition places a restriction on the possible values of the parameter a . 

Then 


f 

•'o /a 2 _ 4B 2 


The substitution x = cos qp in eq. (11) yields 


( 11 ) 


A 


oo 



dx 

/(1-x 2 ) [(a - 2 |c| x) 2 -4(|b| + 2 |d| x) 2 ] 


( 12 ) 


where the limits of integration will be determined by considering the two pos- 
sibilities: u > -1 or u < -1. 


(1) u < -1 


A 


oo 


- 77 

2 / i c | 2 — 4 |d | 2 



dx 


/( 1 - x 2 ) (v - x) (u - x) 


(13) 


and 


a — 2 |b [ 

2 ( [ c 1 + 2|d|) 


(14) 


The integral in eq. (13) can be transformed into a complete elliptic integral 
of the first kind if the location of v with respect to the interval [-1, 1] is given. 
For definiteness let us assume that v < u . Then a substitution prescribed by 
Erdelyi [15] 


u - 1 - 2u sin 2 d> 
x = — 

1 - u - 2 sin 2 


(15) 
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leads to the result 


77 . 1 K J 2(u ..^ 

,/i c i 2 - 4 jd! 2 v/ ' J _ 1 - y rd -u) (-1 -v) 


where K (k) is the complete elliptic integral of the first kind of modulus k . 


(2) u > -1 

Again assuming v < u , we have the subcases v < -land v> -1. Then 
the integral is given by 


/ |C| - t |U| 


ii) A > 0 


(1 - v) (1 +u) 


/(I - u) 
V (1 - V) 


(1 + V) 

(1 + u) 


; v > - 1 


Here the conditions j A) > 2 |b| and 0 < <p < v imply 


a + 2c cos cp > - 2(b + 2d cos q>) 1 

or using the inequalities (7) and definition (14) we get 

cos cp < v . 

Hence if v < -1 the integral will vanish and two subcases have to be con 
sidered: (1) -1 < v < 1 and (2) v > 1. 


(1) -1 < v < 1 

The substitution x = cos cp in eq. (11) yields in this case 


\o = 


2 / c 2 -4 d 


f- 

1 , /i 


1 - x 2 ) (x - v) (x - u) 





As before we assume u > v . The final result will again depend on whether 


u > 1 or u < 1. Performing the appropriate transformations we obtain 

r 


A oo “ 


77 


/|c| 2 - 4 | d | 2 


< 


- — - - k L/-~ v) u < i 

/(i - v) (i + u) \y (i - v) (i + u) j 

== K (y r 

u - V) \1 


(1 + v) (u - 1) 
2(u - v) 


( 21 ) 


(2) y > 1 

A similar procedure produces the result 



dx 

</(l - x 2 ) (x - v) (x - u) 


n 


K 


2 (u - v) 


v/| c | 2 - 4 |d | 2 + v) (u - 1) V O + v) (u - 1) / 


( 22 ) 


We return now to the general integral A pq , for which an asymptotic approxi- 
mation is required when both p and q are large. We follow here Weiss et al.[41] 
who, in a different context, have evaluated a similar integral by making use of 
certain asymptotic results of Duffin [42] . The argument is based on the remark 
that the main contribution to A pq comes from the neighborhood of the stationary 
point of the denominator, located here at (9= 0, qp = 0). Expansion of the 
denominator around this point leads then to 


pq 


J 0 Jo 


cos p<9 cos qqp dd dqp 


a + 2b + 2c + 4d - (c + 2d) <p 2 - (b + 2d - d<p 2 ) 6 2 


(23) 


This can be written also as 


A 


pq 





cos qcp cos p(9 d<9 dqp 
C + D6* 2 


(24) 
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with 


C = a — 2 ( |b| + t c | 2 | d | ) +(|c| +2|d|) cp" 


D - jbj + 2|d| - |d| cp 2 


The range of integration over 6 can be taken as (0, oo) with an exponentially 
small error, since we assume p to be large: 


dcp cos 


,, f 

J n 


cos p 6 d 8 
C + B0 2 


We consider first the case sgnC = sgnD, for all cp's near cp = 0. This implies 
a > 2 |b| + 4 |d| when eq. (7) is used. The inner integral is then an ordinary 
Fourier cosine transform and from Tables [19] we find: 


f dcp cos _qcp tt /T exp j -p/f 1 

D 2 V C | »D 
•'o L J 


7 T cos qcp 


xp X - P J — V dcp . 


Since q is large we expand this result once more, this time in the neighbor- 
hood of cp = 0, and extend the integration over (0,co): 


A ^ f 
” 2 1 


COS H jL exp {- P [p + crcp 2 ] } dcp 
[ R+Scp 2 


in which 


- 2 ( |b | + |cj +_2|d|) 
|b| + 2 | d | 


a d + be 


2(|b| + 2|d|) 3 / 2 [a - 2( |b| + |cl + 2|dhl 1 / 2 


R = [(a - 2( |b| + |c| +2|d|))(|b| +2|d|)] 1/2 ; 

_ (|c| +2|d|)(|b| + 2 |d| ) - |d | [ a — 2 ( |b | + |c| + 2 | d | )] 


We assume now S > 0. Then the integration in eq. (28) can be carried out 
[19]: 



For q — oo the dominant term is the first one in the square brackets. 

The case S = 0 leads again to a well defined result. Since originally we 
have considered values of <p only in the neighborhood of the origin, we disregard 
the singular case arising from S < 0. 

We still have to consider the case sgnC t sgnD, for all cp near the origin. 
This will occur whenever a - 2 (|b| + |c| + 2 |d| ) <0. Then the inner integral in 
eq. (26) exists only as a principal value and we obtain 


A 

pq 




cos qcp 

✓-CD 


sin p 



dqp . 


(31) 


Expanding the functions /- C/D and /- CD around the origin cp = 0 and as 
before letting the upper limit of the integration go to infinity, A pq can be written 
as 


pq 


f co 1 . 9? sin p[o' -cr'cp 2 ] dcp 
2 J 0 R' - S V 


(32) 


in which 

r 


'V s 


'2( |b| + |c| +2|d|) - a 
p" - -I / — ; cr 


aid! + be 


|b| + 2 | d | 


2(|b| + 2 jd | ) 3/2 [2(|b| + | c | + 2 1 d | ") — a] 1 /2 




R' = [(2(|b| + |c| +2|d|) -a) (|b| +2|d|)] 1/2 ; 




o , (|b|+2|dh(|c|+2|d|) + (2(|b| + |c|+2|db-a)|d| 

O — — - 

R' 


(33) 
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The integral in eq. (32) can be evaluated in the sense of a principal value by 


considering the integral of the function 

e iqr sin p |jo' - cr' z 2 ] 
z 2 - a 2 

where a 2 = R' /S' , along the contour shown in the figure. 
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APPENDIX C 


FOURIER SERIES EXPANSION FOR 


S(a - o? + 2b cos d + 2c cos qp + 4d cos 6 cos q>) 

In the following we adapt the one-dimensional procedure of Lighthill[17] for 
expanding a generalized function in a Fourier series, to the two-dimensional case 
and apply this on the function of the title. 

Setting 


S(a - o? + 2b cos 6 + 2c cos cp + 4d cos d cos cp) = 


■L 

m, n-“°° 


»iro 6 i nCp 


(i) 


the Fourier coefficients will be defined by 


f e " inv v £ ds u {£} 


x § (a - o? + 2b cos 0 + 2c cos qp + 4d cos 6 cos cp) (2) 

where the unitary functions U(x), V(y) have the following properties: 


co 

L 


m =-oo 


U(x + nm) = 1; 


CD 

L 


V(y + nrr) = 1 


( 3 ) 


for all x, y, and 

U(x) =0 for | x | > 1 

U(x) + U(x - 1) = 1 for 0 < x < 1 

V(y) = 0 for |y | > 1 

V(y) + V(y - 1) = 1 for 0 < y < 1 
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Let us denote by 6* the real roots of the equation 


a - o? + 2c cos cp + 2(b + 2d cos (p) cos 6 = 0 


( 5 ) 


considered as an equation in 6 . Then 

\Jj + 2 rv 

d\ - ^ r = 0, ±1, ±2, • • • 

- pi + 2 ttt 

where p is taken as the principal branch of the function 

, co 2 - a - 2c cos <p 

W - cos ■ 

2(b + 2d cos cp) 

On the other hand the transformation properties of the 8-function lead to the 
representation 

8(a - cl? + 2b cos 0 + 2c cos cp + 4d cos 6 cos cp) = 


(6) 


(7) 


OD 

•L 


i 


2 b + 2d cos cp 


'we - e + T ) we - e~Y 

| sin 0* | |sin 6 r \ 


( 8 ) 


Inserting eq. (8) in eq. (2) and performing the integration over 6, we obtain: 

A =— I d'5 e" 1 "* Vf-M 

4tt 2 L W”' 


2 I b + 2d cos cp | 


O) 

E 

r =— co 


"im^ + 


/ \ 

~im0~ 

ft 7 \ 

e 

\2tt) 

1 sin 9~\ 


(9) 


Using now the form (6) of 6 * and taking account of the summation formula 
(3) for U (x) , the sum in (9) becomes simply 


2 cos mp 
|sin >/<| 
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Therefore 


A 

mn 



cos mi p 

b + 2d cos cp| |sin \p\ 


e ~in(p y 



( 10 ) 


The function of <p , the nth Fourier coefficient of which is to be evaluated in 
eq. (10), is now an ordinary function and its Fourier coefficients are defined in 
the usual way (Lighthill[17]). Hence 


J 7T 


dcp 


cos m \p 

b + 2d cos cp | |sin \p | 


"inCP 


=j_ r _ 

2t t 2 J 0 l b + 


cos mi p cos ncp 
2d cos cp | |sin \fj\ 


dcp . 


( 11 ) 


There exists a restriction on the integration over cp, since the assumption that 
the principal branch of (7) is to be taken, namely 0 < \p < tt 9 leads to a condition 
on cp: 

_ l < - 3 - 2c cos g < 1. (12) 

2(b + 2d cos qp) 

To find the limits of integration over q> , we must assume something about the 
parameters a, b, c, d (the range of values of a? will follow automatically). We 
make the following choice 

a > 2 |b| + 2 |c| + 4 |d| ; |b| > |c| > 2 |d| ; b, c, d < 0 (13) 


which represents the situation most frequently met in practice. Then solving 
for cos cp in eq. (12) and using eq. (13) we obtain 

v < cos cp < u (14) 


with 

= a - o? + 2|b| . v _ a - a; 2 - 2 | b 1 
2( |c - 2 |d| ) V 2(|c + 2|d|) 


(15) 
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Since the largest possible range of cos cp is -1 < cos cp < 1 , we have to con- 
sider the positions of u, v relative to the interval [-1, l] . A detailed analysis 
of the possible cases (under the assumptions (13) ) furnishes the following re- 
sults (a change of variable cos cp = x is implicit below): 

r 


i-f 

>77 2 I 


dx 


f „( x >; 


^ 77 •'v i/l - x 2 


m , n 


= j 

277 J-j / X _ x 2 


- r 

2T7 2 J~1 


dx 




f (x); 

m , n 


in which 


a — 2 |b| — 2 | c | - 4 |d| < a? < a — 2 |b| + 2 | c | + 4|d| 
-l<v<l<u 

a — 2 |b| + 2 jc| + 4 |d| <cJ < a + 2 |b| - 2 |c| + 4 |d| 
v<-l<l<u 

a + 2 |b| — 2 | c | + 4 |d| < w 2 < a + 2 |b| + 2 |c| - 4 |d| 

v<-l<u<l 

(16) 


f 

mn 


(x) 


cos n (cos 


|b + 2dx| 


x) cos m (co 


sin |cos 1 


- ; - a - 2cx\ 

2' v b + 2dx) / 

a) 2 - a - 2cx \ 

2(b + 2dx) ) 


cos n (cos -1 x) cosm fcos -1 — — ^ CX ) 

_ \ 2(b + 2dx) / (17) 

v/|c| 2 - 4 | d | 2 i/(x - v) (u - x) 

We note that a change in sign of m or n leaves A mn invariant and hence the 
Fourier series for the S -function will be a double cosine series. The theory of 
elliptic integrals [15] shows that with f mn (x) as in eq. (17) the integrals A mn are 
in principle reducible to a linear combination of the three complete elliptic 
integrals, if a transformation to standard form is effected. For small values of 
m, n this can be done without much difficulty, but for large values it is more 
practical to find for A mn asymptotic approximations. 

We give below the explicit form of A 00 . In the evaluation we use the appro- 
priate transformations prescribed by Erdelyi [15] . 
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r 




/2(u - v) 


K 


(/ 


(1 - v) (1 + u) 


2(u - v) 


77 2 /c 2 - 4d 2 




i k(Vzh 

/i - v) (i + u) \r (i - 


2(u - v) 


(18) 


v) (1 + u) 


1 




(1 - v) (1 + u) 


v /2(u - v) \ V 2(u - v) 

The ranges of cJ 2 for which (18) is valid are given in eq. (16). 

We proceed now to find asymptotic expressions for A and consider the 

mn 

three cases: (1) n < m > > 1; (2) m < n »1; (3) m * n » 1. 

Since A mn has different forms in different ranges of cJ 2 , we discuss below 

★ 

only the first integral in eq. (16) which can be written as follows: 

/ a — cJ 2 — 2 |c | cos cp\ 

.eos-iy cosnqpcosm cos 1 — — — 

1 / 1 O / IKI . O I 

I = 1 


2tt 2 /c 2 - 4d 2 


/ 

Jo 


N + 2 [d] cos qpy Hrn 


/(cos cp - v) (u - cos cp) 


(19) 


where 


-l<v<l<u. 


We see that I has a branch point singularity at cos cp = v in the denominator 

mn 

and it is desirable therefore to use an approximation procedure which removes 
this singularity. Such a procedure is provided here by the stationary phase 
method. 

Case 1; n < m > > 1 

Here we consider the integral J_ ^ , 


^nrn 


1 

/c 2 - 4d 2 



g(qp) e imf ('P) dqp 


The two other integrals in eq. (16) can be treated similarly. 


( 20 ) 
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in which 


g(<P) = 


cos ncp 

/(cos cp - v) (u - cos cp) 


f (cp) - COS 


a — c<j 2 — 2 | c j cos cp 
2( |b| + 2 |d| cos cp) 


( 21 ) 


Note that I = 1/2 tt 2 ft e (J mn } . The function f (cp) has a stationary point at 
cp = 0, and we obtain: 

a - co 2 - 2 | c | x (a - o?) |d| + |bc| 

f (o) = COS 1 ; — ; — — — ; f (o) = - 


2(|b| + 2 |d| ) 


/c 2 - 4d 2 /(I - v) (u - 1) (|b| +2|d|) 


( 22 ) 


It can also be easily verified that the function g(<p) has the property 


lim g(<P) - gCo) - Hm (23) 

CP-0 f'(<p) f"(cp) 

Moreover, the function (g(cp) - g(0))/f' (9) has no singularity at cos cp = v and 
therefore is well-behaved throughout the interval of integration. Then all the 
conditions necessary for the application of the method are satisfied, and we find 


^mn 



Ab\ 


2d 


(c 2 - 4d 2 ) 1/4 [(1 - v) (u - 1)] 1/4 


1 ~~ e imf (°) + 0 j • 


(24) 


The result in eq. (24) reflects the fact that sgn f"(0) = -1 and that the stationary 
point occurs at one end of the integration interval. 

Finally 


/|b| + 2 |df | co s m |cos’’ 


a - a? - 2 | c | ' 
2( ! b | + 2 |d|) 


+ sm m (cos 


-I 


a - co‘ 


' 2 - 2 id 


2(|b| + 2 | d | )j 


4/m7r 3 /(a - cJ 1 ) |d| + |bc| (c 2 - 4d 2 ) 1/4 [(1 - v) (u - 1)] 


1/4 


+ O 



(25) 


For this order of approximation the result does not explicitly depend on n . 
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Case 2: m < n » 1 


Here we first make a change of variable: 


cos 9 - 


- 0 ) 2-2 


c| cos cp 


2( | b | + 2 |d| cos <p) 

Then with this new variable, the integral I becomes 

mn 

/ a - a? - 2 |b| cos 6 

. cosmt'cosn Icos 

— f 

/b 2 - 4d 2 *'o 


-•cos 1 v 


I = 


2 C I c I + 2 |d cos 9) 


/(cos <9 - v') (u' - cos 9) 


in which 


(26) 


dO 

(27) 


v' ^ 


a - cJ 2 - 2 I cl 


u' ^ 


- a? + 2 |c| 


2(|b| +2|d|) 2( |b| —2 |d| ) 


(28) 


and 


- 1 < v' < 1 < u'. 

The result follows immediately since now n has the role of m in the pre- 
vious case also with respect to the form of the integrand. Then eq. (25) will be 
valid for this case too, when an appropriate identification of the parameters is 
carried out. 

Case 3; m % n » 1 

The author has been unable to find a meaningful asymptotic expression for 
this case by the stationary phase method. The method of steepest descent is not 
helpful either since here it becomes identical with the phase method. 
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APPENDIX D 


PARTIAL FRACTION DECOMPOSITION AND SUMMATION OF FINITE SUMS 

Let F (A.) be a polynomial of degree N, and f (k) one of degree M < N . Assume 
that F(/v) has j distinct roots A k with multiplicities p k , such that 


F(M = F N ~]""|~(\-\ k ) Pk . (!) 

k = l 


Then it is well-known that the following decomposition exists and is unique 


f (^) _ _[m § 
F(M F N 


M.N 


+ 


j 




LL 

k— 1 r = l 


c 


rk 


- KY 


( 2 ) 


where C rk are given by 


c rk = 




(P k - r)! d\ Pk 


f (\)(\-X k ) 
F(\) 


A. = K, 


(3) 


For simple roots, all p fe = 1 and the formula (2) becomes 


F(M F n 


N 





(4) 


157 



with 


q, = 


f(\) 

F'(\) 


(5) 


In the sequel we shall have use only for the simpler case (4) with the further 
restriction M < N . Several formulas are of particular interest: 


f'(\> _ y i 
n\) 4_. X - \ 


( 6 ) 


By repeated differentiation, eq. (6) yields 

N 

(- 1 )"- 1 d "- 1 y 

(n - 1 )! dX n - l \F(k) ) L -* 

The result in eq. (7) can be extended by analytic continuation to all real n , 
the l.h.s. of eq. (7) being transformed then into an appropriate fractional deriva- 
tive [43] of logF(-V). 

These results suggest the possibility of their validity also for functions other 
than polynomials, e. g. , functions possessing a discrete infinity of zeros which 
satisfy certain conditions of convergence. We shall not pursue this subject further 
here except for the following, well-known example: 

Let F(X.) = s in A.. Then formally 


" K) n 


n > 1. 


(7) 


sin k ~ 


+00 

| | - Tin) 

n=-oo 


( 8 ) 
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and 


F'(M 

F(\) 


+00 

COt k = 

n=-oo 



1 

X. - 7rn 


( 9 ) 


A rearrangement is performed to produce convergence with the final result 


cot A . = _ + 2 A . 

k 


L 

n~ 1 


k 2 - tt 2 n 2 


( 10 ) 


otherwise obtained from complex function theory. 

If the \ k 's are the characteristic roots of a matrix a n , then eq. (7) gives a 
formula for evaluating the trace of (a n -\I) -n , F(\) being the characteristic 
polynomial of A N 


F(\) = | A N - \ I 


( 11 ) 


The previous treatment can be extended in still another direction. Consider 
a sum of the type 



Mr 


( 12 ) 


To evaluate this sum we define a function S (1) (k) which coincides with 
when k-k : 


IN 

s (i > (k) =j2 

Ml 


k-k. k-k 


= 

F(M 


k-k 


(13) 
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The last result is obtained by substituting for the unrestricted sum from 
eq. (4), with M < N. S (1) (A.) can be rewritten as follows 


S<‘> (A.) = 


f(M -c r 


F(M 

\ - \ 

r 


FM 


(14) 


It is to be noted now that S (1) (A.) becomes an indeterminate expression of 
the form 0/0, when k -* k T . This is so because of the form (5) for the coefficients 
C f . Using the L'Hopital rule repeatedly we find 


S (i> _ f'(K) 
r F*(\ r ) 


F(A- r )F"(A. r ) 

2[F'(A r )] 2 


(15) 


A particularly simple relation is obtained for the case f(A) = F'(A) : 


N 





1 F"(X r ) 

2 F' (k r ) 


(16) 


The evaluation of sums of higher powers 



Mr 


(17) 


cannot be done simply by differentiating (15) or (16) with respect to k T , since 
F'(A r ), F"(A. r ) depend on k t not only through their argument but also through 
their coefficients. It is much simpler therefore to use again the same procedure 
as above. We illustrate this for the sums (17) with n = 2, S ( r 2) . 
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We have to evaluate the expression 


S< 2 > 


1 im 



C k 


(\-\ r ) 2 


1 im 


<--£ (in)- Cr \ 

dk \F (k)J (k _ K )2 

V T J 


1 im 


f(k)F(k ) _ f(\)F*(\) +c rm T 

L^-^r. 

[F(\)3 2 


(18) 


Again the limit process leads to the indeterminate form 0/0 and we make 
use of the L'Hopital rule several times to obtain: 


f-(\ r) f(^ r )[F"(\ r )] 2 

2F*(k r ) 2[F'(\ r )] 2 4[F'(\ f )] 3 


f(k T )F"'(k T ) 
+ 6[F'(^ r )] 2 


(19) 


For the case f (k) = F' (k) we obtain the simpler expression 


N 



Mr 


1 F*(\ r ) 

2 F'(\ r ) 


1 F w (\ f ) 
3 F*(\ r ) 
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For powers n > 4, the formulas become quite complicated. In Chapter IV 
the results in eqs. (16) and (20) are used to evaluate certain perturbation sums. 
Since one seldom considers perturbation approximations beyond the fourth order, 
the form of S^ n ) for arbitrary integral n is not needed. 

The results in eqs. (16) and (20) can be extended to the functions discussed 
in the paragraph preceding eq. (8). Thus if again F(\) = sin A. and in eq. (20) we 
take X. =0, then 

r 5 


oo 



or 


00 



a well-known result, here following immediately from 


F'(0) = 1; F"(0) = 0; and F'" (0) = - 1. 


In connection with applications below and in Chapter IV it is important to note 
that if the polynomials F(X) and f (A.) are known in c losed form and if the roots 
\ k are explicitly known also, then the formulas (4), (6), (7), (16) and (20), give 
closed form expressions for the corresponding finite sums. 

We proceed to apply some of these formulas to the exact evaluation of one- 
dimensional sums related to Green's functions of matrices appearing in the text. 
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The Green's matrix associated with a given matrix A is by definition the 


matrix G 


G (X) s (A - X I) -1 . ( 21 ) 

If X k and T denote respectively the eigenvalues and diagonalizing matrix of A , 
then the elements g r .(X) of G are given by 


IN 

<M = £ 

k = l 


T j - 1 
rk A kj 

- k 

k 


( 22 ) 


where we assume X not to coincide with any of the eigenvalues of A. 

It is readily seen that the sum in eq. (22) is of the type exhibited in eq. (4), 
with M < N and = - T rk T^ 1 . If a polynomial f r i (X) can be found such that 


T 


rk 



F'(M 


where F(X) is the characteristic polynomial of A, then 


(23) 


g rj M = 


f r,(M . 

F(X) 


(24) 


To illustrate the procedure we consider the Green's functions associated 
with a monatomic linear chain with fixed ends. 
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The dynamical matrix of this system is of the form shown in eq. (2) of 
Appendix E . Inserting in eq. (22) the \ and T r ^> there given, we find 




N 


IN 

- £ 
+ 1 4— < 

ksl 


COS 


7T(r - j)k 
N + 1 


cos 


77(r + j)k 

N + 1 


> ( 25 ) 


. 77k . 77k 

a — A. + 2b cos a - A. + 2b cos 

N + 1 N +T 

r, j = 1, . • • , N 


Eq. (25) shows that we have to evaluate the sum 


S(m, N) = _ 

N + 1 


L 

k = l 


COS 


77 m k 

N + l 


\ - a - 2b cos 


7ik 

N + l 


(26) 


Now eq. (7) of Appendix E gives for the characteristic polynomial of A 


yN+1 y-N-1 

= (- i) N b N+i r v 


U 


(27) 


with 


U = /(\ - a) 2 - 4b 2 ; V = X Z 3 + V. 

2b 


(28) 


It is not difficult to verify that if we take F(A.) and f (k) in eq. (4) as follows 


F(\) = U 2 J9 n (\) ; f (\) = v N+1_m + V 


-N-l+m 


(29) 
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where for the present 0 < m < N + 1 , then 


C. = = 1 cos 7Tmk ; k = 1 N (30) 

F'(X k ) N+l N + l 

and 

C* , = (*!)" * (31) 

F'(\*) 2(N + 1) 

where 

X* = a ± 2b . (32) 


Finally we can write 


c/ MN 1 v N+1 - m + v- N ~ 1+B1 

S(m, N) = _ 

U yN+i - V -N_1 


1 1 + (-l) m 

2(N + 1) |_X - a - 2b X - a + 2b_ 


(33) 


Since m' = N + 1 - m satisfies 0 < m' < N + 1, eq. (33) will be valid also for 
m', with the result: 


N 

S' (m, N) = 1 T 

N + l Z— i 


(-l) k cos 


7rm k 


N + l 


X - a - 2b cos 


77k 

nTi 


1 v m + v~ ra 
U v N+1 - V " N * 1 


2(N 


1 1 + (-l) N+1-m " 

+ 1) _A - a - 2b X - a + 2b_ 


These coefficients ore related to the polynomial in eq. (29). 


(34) 
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I 


The range of m in eq. (26) can be extended now to all integers by using eqs. (33) 
and (34). 

For m = 0 we obtain from eq. (33) 


— r 

N + 1 L — i 

k = l 


1 


X. - a - 2b cos 


77k 

N + 1 


i_ x - a _i v N+1 + y-N -1 

N + 1 u 2 U V N+1 -V""' 1 


(35) 


It is easily seen now that 


r S(r - j, N) + S(r + j, N) 
r + j < N + 1 

- g rj = *S (36) 

S(r - j, N) + S' ( r + j -N - 1, N) 

N + 1 < r + j < 2N 

Other sums with closed form expressions can be obtained from (33) and (34) 
by addition, subtraction or proceeding to appropriate limits with the parameters 
involved. For instance, in the sum (33) the left hand side can be regarded as the 
Riemann sum associated with the integral 


7 T 

cos m 6 , n 

— d U 

k - a - 2b cos 9 

and letting N “* 00 we obtain the value of this integral on the r.h.s. The value will 
depend on the position of V with respect to the unit circle. 
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We also remark that these results will be valid, if none of the denominators 
in the sums considered vanish, regardless of the nature of the parameters. If 
a, b were operators (e.g., matrices), the sole requirement for validity would be 
that these commute. 

In principle the preceding treatment applies to finite sums in any number of 
dimensions. However, in practice, the application of these results is limited by 
the fact that the explicit forms of the appropriate polynomials are not known. 

In most cases one summation can be performed exactly but the rest have to be 
treated according to the methods of Chapter in. 
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appendix e 


MATRICES ASSOCIATED WITH PERIODIC SYSTEMS 
Most of the material to be presented in this Appendix, while well-known, is 
rather scattered throughout the literature. Moreover the presentation there is 
frequently fragmentary, partial results being given as needed for specific appli- 
cations. Also most of these results are obtained by determinantal methods 
which do not yield easily explicit expressions for the eigenvectors. For this 
reason the treatment to follow is done mainly by matrix methods, as these allow 
us to take into account more conveniently the various symmetries involved. The 
matrices to be subsequently listed are those met in the text together with some 
variants and generalizations which, to the author's knowledge, have not been 
published elsewhere. Each one of the matrices presented is accompanied, when- 
ever possible, by its eigenvalues, eigenvectors and characteristic determinant. 

The matrices will be classified according to the dimensionality of the 
problem from which they arise. 


Part I. One Dimensional Matrices 


1. Continuants 

These matrices belong to particular subsets of the class of Jacobi matrices, 
the latter being defined by 


/ a i 


0 \ 


V N , 

C V N b , 

V n - 1 

x \ 

''O. 

\o V, / 


(i) 
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(i) Simple Continuant 
This matrix is defined by 


a b 0 

\ V 

X N 

b ' 'b 

x X 

\ 0 b v a 


( 2 ) 


Its elements can be written in the form 


(^ N )i j - b8 i-l,j + a S ij + b8 i,j"l : ’ N 


The characteristic polynomial is the determinant 


9 n e>- l A »- x 1 


It can be easily verified that J9 n (k) satisfies the recurrence relation 




(3) 


(4) 


(5) 


The solution can be found by putting J9 n = x N and solving the quadratic equation 
in x that results. The solution is 


1 


K (M = 

/(a - k) 2 -4 b 2 


a-A. + i/(^a-\) 2 -4b 2 \ _( a - k - /(a - A.) 2 -4 b 2 


2 / \ 2 
This can be written in a more symmetrical form as follows: 


9 „W = (-l> N b w * 1 - V " tl ~ Y - 


■N- 1 


where 


U = v/(\- a) 2 -4 b 2 ; V = 


U 


_ \-a +U 


N +1 


(6) 


(?) 


2b 


(8) 
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Various representations can be given to $ N (\) when certain relations are 
imposed on a, A. and b. For instance, if one can write 


a - A = 2b cos cp ; 0 £ cp < tt 

then 

K = ("1) N bN — - ( - + 1} ^ (9) 

N sin cp 

The polynomial sin ^ + 12 ( cos 1 is known as the Chebishev polynomial 

sin (cos - 1 x) 

of the second kind. 

The eigenvalues of A N can be immediately found from $ N (A.) = 0 by using 

eq. (7): 

X, = a + 2b cos — — - ; k = 1 , . . • , N (10) 

K N + 1 

The eigenvectors x k = (x lk - • \x Nk ) can be found by putting x jk = A k <f j 
in the eigenvalue equation 

<*N ‘ I) • = 0 < U > 

and solving the resulting quadratic equation in E, . 


The normalized eigenvectors thus obtained are the columns of the matrix T 
which brings A N to diagonal form: 


t r 2 • 77 j k 

)k Vn + i n + i 

T is a symmetric orthogonal matrix, T - 1 = T , such that 


( 12 ) 


T- 1 A n T = A n (13) 

with 

(\) lk = \ v < 14 > 

As T does not depend on the elements of A N , the simple continuant matrices 
commute, since all of them can be brought to diagonal form by the same similarity 
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transformation. In spite of this property, these matrices do not form a group 
with respect to multiplication since their powers are not simple continuants. 

(ii) Asymmetric Simple Continuant 
Let A n be here given by 


A 


N 





0 



\ 



(15) 


where the b. ,c. satisfy the relation 

b. c. = b 2 for i = 1, • • *, N - 1 (16) 

In particular for b. = c . (complex conjugate), a n will be a hermitian matrix. 

The eigenvectors of this matrix are the columns of the matrix ST, with S a 
diagonal matrix the elements of which are 


S jj = ( c i c 2 ' ' * c j - i bj b. + j ••'b N + 1 ) 1/2 

and T is the matrix given in eq. (12). 

The eigenvalues are 


(17) 


L = a + 2b cos — — ; k = 1, • • N (18) 

N + 1 

The characteristic polynomial (or determinant) is as in eq. (7). 


(iii) Generalized Continuant 

If one inquires what is the matrix with eigenvalues of the form 

N 

X k =a o + 2 22 cos ¥TT + ( - 1)ia N + i ;k = 1 ’ '" )N < 19 > 

j = 1 
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and with the diagonalizing matrix T of eq. (12), one obtains the following inter- 
esting result: 


/ a ° v---vA 




V 


Symmetric 


'S / 




I / 

I r 


?N 3 N + 1 


I / 

a N 


V 


✓ 

'n + i 


Symmetric 


(20) 


This matrix can be decomposed as follows: 


N+ 1 

A„= £ a.»«> 

j = 0 


where 


A (°) = I N ; A(»> 


/ 0 N 1 0 ' 

x A 

\o''l ''of 


. A (N + X) =. 


0 


1\ 


4 


K 0 / 


(21) 


A (j) = 
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Not all of these matrices are independent, in fact 


j = 0, - - * 


N . 


N even 


A<N+l-i) = _ A (N + D A (i) J 


i = 0, * * •, — — - ; N odd 
2 


(23) 


The number of independent matrices is (N/2) + 2 for N even and (N + l)/2 + 2 for N 
odd. 


An important particular case is obtained when we put a . = 0 for j = 3, . . . , N + 1 


/ 




\ 


0 



with eigenvalues 


, . 77 k - 2 77k 

K. = a. + 2 a, cos + 2 a cos _ — - 

k 0 i N + 1 2 N + l 


k = 1, • • • , N 


(24) 


(25) 


The determinant of A N - A I with A N given in eq. (24) can be easily evaluated 


by writing 

N 

K-mI =TT 

k = 1 

Using the identity 


cos 


2 77k 

N + l 


= 2 cos 2 


77k 

N + l 


-1 
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and the form (25) of ^ k , we can write 


where 



+ 2bj cos 


77 k \ 

nti; 



+ 2 b„ 


cos 


77k \ 

ntt; 




-4c (a - 2 c - \) 


2 i / c — 


b + ^b 2 - 4 c (a - 2 c - \) 


2 v/c~ 


bj = vc ; b 2 = /c” 


Then the characteristic determinant is 


|a n -M| =^( x > (MJ9<, 2 > (M 

in which the JQ^ 1 ) (X) are given by eqs. (6) or (7) with a - X replaced by a i and 
b by b. . 

(iv) Alternating Continuant 
Let A n have the form 



where the last term on the main diagonal is v for N even and u for N odd. 
Accordingly, we treat these two cases separately. 
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(1) N = 2 n 

One can write the following decomposition for A 2n : 



We apply the similarity transformation T of eq. (12) on A 2n , with the result 


T " 1 A 2n T = 


A', 


2n 



R 


2n 


where 


and 


i = 


u + v 
2 


+ 2b cos 


7rk \ 

2n + 1 j 





(28) 


(29) 


(30) 


The matrix given by eq. (28) belongs to the class of matrices {Z} defined by 



It is not difficult to prove that a similarity transformation P brings Z to the 


form: 




(32) 


where P is the permutation matrix 



8 . 


j , 2k “ 1’ k = 


8 . 


j , 2(2n-k +1)’ 


1 » • • • * n 
k = n + 1, 


. , 2n 


(33) 
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and 


Z< k > =| ) ; k = 1, • . • , n 

y 2n-k+l X 2n-k+l 


(34) 


The eigenvalues of Z 2n are then 

. x, 

^(Z) 


X k + X 2n~k+1 




+ y k y 2n~k+l 


(35) 


k = 1, • - • , n 


Using the above results we find for the eigenvalues of A 


2n 


*■ 


U - v\2 


2b cos 


7rk 


2n + 1 


(36) 


k — 1 , • n 


and these appear in the order A.* , A.“ , A.* , A." , . . . , A” • The eigenvectors of A 2n 
are the columns of the matrix TPS, with T given by eq. (12), P by eq. (33), and 
S is the matrix 

r \ 

S = S< v 2 > (37) 

V V 'S <">/ 

in which S (k) is a 2 x 2-matrix diagonalizing the matrix A (k) , 


/u + V 7Tk 

' — + 2b cos 


u - v 


A< k > = 


2n + 1 


u - v 


u + v 7r(2n + 1 - k) 

— I — + 2b cos ) 

2 2n + 1 / 


(38) 


The characteristic polynomial of A 2 n can be shown to have the form i9 2n (A.) 


if a - A. is replaced by /(u-\)(v-A) in eqs. (6) or (7). 
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(2) N = 2n + 1 


Making a decomposition similar to eq. (27) and following the same procedure 
as above we find that the eigenvectors of A 2n+1 are the columns of the matrix 
TPS, where T is given by eq. (12), P is defined as 

f 8 j,2k-i ; k = 1 " +1 

p kj = J (39) 

[ 8 j,2(2n + 2- k ) ; k = n + 2, .... 2n + 1 


S is the matrix 



where S (k) is the 2x2 matrix diagonalizing the matrix A (k) , 


/u + v 7rk 

— _ + 2b cos 


u - v 


2n + 2 


A< k > 


u - v 


u + v 0 , n( 2n + 2 - k) 

2 2n + 2 / 


(40) 


(41) 


The eigenvalues are 


* -H 1 * 1 ' ( 


2b cos 


7rk 

2n + 2 


k = 1, . . . , n 


(42) 


\ 


2n + 1 


= U 


and their order is , . . . , K + , X. , , , . 

1 1 n n 5 2n + l 

The characteristic polynomial of a can be found by using the recur- 

2n +1 

rence relation 


U 2n+1 -A.l| + b 2 | A 2n _j - XI | =(u-\) J9 2n (M 


(43) 
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where J9 2n (\) is the characteristic polynomial of A 2n described in the previous 
subsection. The solution is 

, , A 

|a 2 „ + i -\I| = (- 1)" (u - K) b 2 " £ 1)k — (44) 

k=o b 2k 

Performing the summation in eq. (44) by utilizing eq. (7), we obtain 

|4 *.« - KI 1 = /B V. <M (45) 


where J9 + (\) is given again by eq. (7), with a - A. replaced by V(u -k) (v-k). 


(v) Alternating Continuant Generalized 


Let A 2n+1 be the matrix 


■*2n+l 


/U 

b 



(46) 


To find the eigenvalues and eigenvectors of this matrix we make use of a 
decomposition of the eigenvalue equations, first employed by Born [25] , which 
is of some independent interest. Let the equations mentioned be 




where x is a 2n + 1-dimensional column vector. 
We rewrite the eq’s (47) as follows: 


b 

N 

C \ 


(u - X.) I n • X 


o _ 


N 


0 


\ 


x N 
x N 


. X 


v c N b 


(47) 


(48) 
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The n-dimensional column vectors x° , x e are defined by 



Inserting the value of x° from eq. (48) into eq. (49), we obtain 



Using eq. (50), we finally obtain: 



with 

a = b 2 + c 2 - (u - A.) (v - k) ; /3 = be (54) 

The matrix appearing in eq. (53) is a simple continuant and therefore equat- 
ing to zero its known eigenvalues we get n quadratic equations for X : 

b 2 + c 2 - (u - k) (v - k) + 2bc cos — =0; k = 1, . . • . n (55) 

n + 1 
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We have found in this manner 2n eigenvalues of A 2n+1 ; 






+ (b - c) 2 + 


2 /be cos 


2n + 2 


(56) 


k = 1, . . . , n 


For b = c these eigenvalues reduce indeed to those given in eq. (42). Also we 
have found for the components of x® : 


x e , = sin — - ; j, k = 1, . . . , n 


Jk 


n + 1 


and 


«£*■ 1 


V k -u 


/X . \ 


\ 

C \ 
\ 

\ 


\ 


X \ 

i 0 \ \ / 

\ \ bl 


x£; k = 1, . 


. n 


(57) 


(58) 


The remaining eigenvalue is easily seen to be equal to u , 


^2n + l = U 


(59) 


with the eigenvector 


X 2 k( 2n + 1 =0; k = 1 * n 


2n + l 


2k+l,2n+l 




k — 0, • • • j n 


(60) 


Note that the eigenvectors we found have not been normalized, but the normaliza- 
tion (if needed) presents no difficulties. From the mode of solution it is seen 
that the characteristic polynomial of A 2n+1 is given by: 

|A 2n+1 -kl\ = (u -M J9 n (M < 61 > 

in which $ n (\) is as in eq. (6), with a - A. replaced by b 2 + c 2 - (u-\)(v-A.). 
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For the matrix of even order & , 

2n 



no analytic results are available, since the eigenvalues cannot be explicitly found 
for b f c . 

(vi) Variants of the Continuant 

( 1 ) 



Using methods similar to those employed previously for the simple contin- 
uent, we obtain the eigenvalues 


K - a + 2b cos — ; k = 0, . . . , N - 1 (64) 

k N 

and the normalized eigenvectors as columns of the orthogonal matrix T 
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Note that all matrices A N of this type commute as they are brought to diagonal 
form by the same similarity transformation. 

The characteristic polynomial is 


\\ - M| = (a - X + 2b) ViM (66) 

where J9 n _j (X.) is given by eq. (6). 



The eigenvalues are 

X. = a + 2b cos Zjl ; j = 1 N (68) 


and the diagonalizing orthogonal matrix T is 




T PT . 77 (k - 1/2) 

T -J — sin — 

kj T N N 

j i N 
(- l) k 


k = 1, . . . , N 


> 


(69) 


T - 

kN 


/N 


J 


The characteristic polynomial is given by 


1^, -Ml =(a -X-2b) < 7 °) 

in which (X) is again as in eq. (6). 
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( 3 ) 


j a + b b 


b a ^ 

\ \ \ 

\ \ \ 

\ \ 

\ a 

0 \ 


’ b a - b . 


The eigenvalues are 


k. = a + 2b cos 77(3 ~ 1/2) ; j = 1 N 

J N 


The eigenvectors are columns of the orthogonal matrix T 


T / 2 _ 77(k - 1/2) (j - 1/2) 

T ki =/ n C ° S N 


k,j = 1, . . . , N 


The characteristic polynomial is 


\\-kl \ =J9 N (\) -b 2 J9 N _ 2 


a ± b b , 


b a 


\ 0 


\ \ ' V 

\ \ x 
x \ \ 

\ a b 

0 \ 

\ 

b a 
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Eigenvalues 


\ + 


= a + 2b cos 


77(2k - 1) 

2N + 1 


-\ 


k = 1, • • . , N y 


(76) 


K 


= a + 2b cos 


2irk 

2N + 1 


J 


The diagonalizing matrices T + , T“are orthogonal 




r _ 2 77(2k - 1) (2 j - 1) 

kj “/2f rn 2(2N + 1) 


k, j = 1 N 


\ (77) 


T l 

kj 


2 77C2k-l)j 

sin 

2N + 1 


/ 2N + 1 


and 


J 


I A* - \I| =(a-Ub) J9 n _i(M - b 2 i9 N _ 2 CM 


(78) 


(5) 


2n 


a + C 

b 


a c 0 
\ \ 

c \ \ 

C \ \ \ 

\ \ \ 
'c a b 


(79’) 


b a + c 


The elements on the first upper and lower diagonals are altematingly b and c . 
This matrix is associated with a monatomic linear chain with free ends and 
alternating spacing constants. 
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We give below the derivation of the eigenvectors and eigenvalues since it 
illustrates a frequently used device in problems with similar symmetries. 

The eigenvector equation 


(A 2n - \I) • x = 0 

when written in full, leads to two prototype equations 

bx 2 j - 1 + (a - ^ x 2j + CX 2 j + 1 = 
cx 2j + (a - K) x 2j + 1 + bx 2j+2 = o 

j = 1, . . . , n - 1 



(80') 


(81') 


and two boundary conditions 

(a - k + c) Xj + bx 2 = 0 
bx 2n _i + (a - \ + c) x 2n = 0 




> 




We put 

x 2j _i = A y 2j_1 ; x 2j =By 2 i 
and insert in eq's (81') to obtain 


(82') 


(83') 


(b + cy 2 ) A + (a-A.)yB = 0 

(a-\)yA + (c+ by 2 ) B = 0 
A non-trivial solution for A,B requires y to satisfy the equation 



(84’) 


(b + cy 2 ) (c + by 2 ) - (a - A.) y 2 = 0 


(85') 
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or 


y± 


= ^ ( a ~ M 2 - b 2 - c 2 + -J £(a - A .) 2 - b 2 - c 2 ] 2 - 4 b 2 c 2 ^ ( 86 ') 


Let us put now 


(a - A .) 2 - b 2 - c 2 = 2 bc cos 26 


(87') 


Then 


y± = e ±i5 


( 88 ') 


Inserting this in eq. (84') we find 


hp"^ i 

B De + ce . 

1 a - A. 1 


_ be ‘ 3 + «-■' A 
* a - A. £ 


> 




Introducing now 


(89') 


x -A e^i'D^+A e _i < 2 i -1 > ff ^ 

x 2j-l - n l e + 2 e 1 


i2)9 


(90') 


x 2j - Bj e i 2 i* +B 2 e“ i: 

into eq. (82') and inaking use of eq's (87') and (89'), we finally obtain a system 
of two equations for A x , A 2 . The requirement of a non-trivial solution for these 
yields an equation for 6 : 


sin 2 nd [b 2 - (a - A. + c) 2 ] = 0 ( 91 ') 

Then, by inspection of the original matrix, we find that the only admissible roots 
of sin 2ni9 = 0 are 
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f) _ 7fk . k _ l 

^ k -_, k-i,. 


• , n - 1 


(92') 


with the remaining two roots furnished by b 2 - (a-\+c) = 0. 


The eigenvalues can be finally given as 


= a ± i J b 2 + 2bc cos — + c 2 ; k = 1, • • • , n - 1 




r ( 93 ’) 


k 1 = a + c ± b 


J 


For each k there exists one linear relation between Aj and A 2 . Eliminating 
one of them and introducing its expression in eq's (90'), we obtain the components 
of the respective eigenvector. The matrix T 2n which brings A to diagonal form 
should have the normalized eigenvectors for its columns. 

( 6 ) 




/2 b 


0 


v2b a b 
\ \ 

\ \ 

K V N 

\ \ \ 

\ N b 
° \ \ 

'b N a 


The eigenvalues are 


k. = a + 2b cos 77 ^ — — ; j = 1, . . . , N 

J N 


(80) 


The eigenvectors are columns of the matrix ST, where 


T kj = cos - n(i - 1/2) . kfj= , N 


(81) 


and the matrix S is given by 
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/2 8 .. 

1 1 


S =< 

i j A 


s ij ; j =2 N 


(82) 


Note that while the rows of T are pairwise orthogonal, this not so for the columns 
(eigenvectors), and hence T is not an orthogonal matrix. 


The characteristic polynomial is 


|a n -M| = yv) -b 2 J 9 n (\) 


(83) 


where again J9 (\) is of the form (6). It is of interest to note that if one puts 
a - \= x, then it can be shown that |a n - \I | is proportional to the Chebishev 
polynomial of the first kind, defined by 

3 n (x) = cos N (cos -1 x) (84) 


( 7 ) 



(85) 


This matrix is associated with a diatomic linear chain with free ends. We give 
below the derivation of its eigenvalues and eigenvectors since again this exhibits 
an approach useful for problems of this sort. 


The equations to be solved are 


(A 2n - \I) -x = 0 


( 86 ) 
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where x is a 2n -dimensional vector. We divide now each equation in (86) by b and 
define: 




U - k ' 

b ’ 



We make also the substitution 


x = S-y 


(87) 


( 88 ) 


where the diagonal matrix S is given, in an obvious notation, by 

S = diag(v / 7 7 , /£, . . • , /tj , J (89) 

By inspection it is clear that neither u nor v are eigenvalues and hence £ , V / 0. 
We multiply now the matrix A 2n - A.I on the left by the diagonal matrix S' defined 
by 

S' = diag ( — , • • • , > ~ ~ ^ 

\ /£ /17 kf /t] / 

This operation does not change the eigenvalues because the null space of ( A 2n - 
\D S is identical with that of S' ( A 2n - A.I) S (the null space of S' contains only 
the 0-vector). After all these transformations, eq. (86) becomes: 


in which 



a = / ^77; 



(91) 


(92) 
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Eqs. (91) can be solved by putting 


y r = y r ; a = - 2 cos 0 (93) 

in the rth equation of (91). We find then 

y r = Ae ir(? + Be" ir5 (94) 

where the constants A , B will be determined from the first and last eq's in (91). 
The process of elimination finally yields the eigenvalues equation 

(ae' ifi - 1) (ae ie - 1) ( e i2n 9 - e~ i2n0 ) = 0 (95) 

By referring to the original matrix it can be shown that the only admissible 
solutions of 


e i2n 6 _ e “i2n 0 _ q 


(96) 


are 


9 - ' k - 1 

e k -_, k - 1 . ■ 


. , n - 1 


(97) 


The two remaining values of Q are solutions of the equation 

a? - 2a cos (9 + 1=0 (98) 


Eq. (97) yields the eigenvalues 


( 2bc ° s ^j ; k=i "- 1 <99) 

Eq. (98) when expressed in terms of g , 77 becomes 


+ € + T) = 0 


( 100 ) 
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with the solutions (the so-called surface modes) 


v± u + V 
\=— 



(101) 


The eigenvectors corresponding to the eigenvalues (99) are S* *y^ k) where 

/ 0 \ / \ 
i 

2n 

77 k (r - 1) 


v OO = a* 
y± - 


sm 


r sin 


sin 


2n 


sm 


77k(2n - 1) 


sin 


77 k 

2n 

77 k r 


2n 


7 


sm 


0k(2n - 1) 
2n 


( 102 ) 


and S* is the matrix S of eq. (89) with k replaced by , while those corresponding 
to the eigenvalues are S* -y£ 0 } where 


and 



(103) 



k = 0, 1, . . . , n - 1 


(104) 


The eigenvectors presented above are not normalized. The characteristic poly- 
nomial is 


K- m l = 




- b (/7T^ + /^-TX ) j9 2n-X^ ) +b2 j9 2n-2^> 


(105) 
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where all the polynomials J9 j in eq. (105) have the form shown in eq. (6), with 
a - k replaced by /(u -A.) (v - k) . 

We note that such analytic results do not exist for the matrix of odd order 


A 2n+1 ~ 



(106) 


though its characteristic determinant can be written down without difficulty. 

2. Associated Jacobi Matrices 
i) Let a n denote the matrix 



We exhibit below a similarity transformation which reduces a n to a diagonal block 
matrix composed of simple Jacobi matrices. We treat separately the cases N = 

2n , and N = 2n + 1 . 
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We define a permutation 


7T 


2n 


/ 1 2 3 4 5 • . . 2n 

\ 1 n + 1 2 n+2 3...2n 


( 108 ) 


When we apply this permutation on the columns of the unit matrix I 2n , we obtain 
a permutation matrix P 2n with the property 



(109) 


Note that P 2n 
transpose of P 2n . 


is orthogonal so that P“* = P 2n , where P 2n 
The elements of P 2n are given by 


denotes the 


(P 


2n^i j 


s. 


i.2j-l 


s. 


i,2(j-m) 


j = 1, . . . , n 


j = n + 1, • ■ • , 2n 


(HO) 



(2) N = 2n + 1 


On applying the permutation 


77, 


2n+l 



(111) 


on the columns of the unit matrix I 2 we obtain the permutation matrix 
which brings A 2n+1 to block-diagonal form: 


r s. 


i,2j-l 


j =1, . . . , n + 1 


^*2n+Pij “ 


S. 


i , 2 ( j ~n“ X ) ’ j - n + 2, • . • , 


2n + 1 


( 112 ) 


and 


P _1 A 

* 2n+l “2n+l 



(113) 


These reductions enable us to obtain eigenvalues and eigenvectors for a n in 
eq. (107) if the resulting block matrices belong to any one of the classes of 
matrices previously described. 
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(ii) If A n denotes the matrix 



it is still possible to effect a similar reduction. This is done by assuming N = i p + r 

(i.e., N = r mod(p)), where 0 < r < p - 1. The appropriate permutation here is 

Tf 0 defined by 
ip + r 


pk -. (p - 1) i i r + k; k = 1, • . . , -t 


n = 1, 


r =$ k = 0, . 


I 


pk+n-(p-l-n)'£+k+l 


> (US) 


n = r + 1 , • • •, p - 1 ^ k = 0, . •• , i-lj 

The permutation matrix P , obtained by applying rr^ p + r on the columns of the unit 
matrix, reduces the matrix a n to a block-diagonal form in which the first r blocks 
are (i + 1) x (l + l) Jacobi matrices and the remaining p - r are ^ -matrices 
of similar type. 
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3. Circulant Matrices 


i) Simple Asymmetric 


Let A„ be the matrix 

N 



(116) 


The rows of this matrix are cyclic permutations of the first row and hence the 
notation 


A N “ < s 0 S 1 S 2 * ■ * S N-l^cyc 

is self-explanatory. The elements of a n can be written also as 


(117) 


r-k 


r > k 




k, r = 0, . . . , N - 1 


(118) 


S N-(k-r) k > r 


The eigenvalues and eigenvectors can be found in several ways but the most elegant 
procedure is based on the observation that A N can be written as a linear combina- 
tion of powers of a single matrix, 


n-i 



s. 

i 


(119) 
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in which Q is of the form 


r 1 

f 0 10 

\ 

\ 

. \ 

\ 

0 0 \ 
\ 1 0 0 


( 120 ) 


Q is a permutation matrix arising from the application of the cyclic permutation tt c 


/ 1 2 3 • • • N - 1 N 

\ 2 3 4 . . . N 1 


( 121 ) 


on the columns of the unit matrix. The powers of v c form a (commutative) sub- 
group of the permutation group, and this property is shared also by the matrix 
representation Q. Since (tt c f* =3 (the identity permutation) we obtain 


qn = i (122) 

This shows immediately that the eigenvalues of Q are e i27Tk/N , k = 0, 1 

N - 1. A simple calculation then yields also the eigenvectors. Returning to a n 

we can write for the eigenvalues 


N-l 

\ s. e i2wk i' N ; k=0, ...,N-1 

j=0 


(123) 
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The normalized eigenvectors are the columns of the unitary matrix U which 
brings a n to diagonal form 

i 

U r/ p = — e N ; r, i = 0, .... N - 1 (124) 

/N 

It is clear, from eq. (119) or the form of U, that any two circulant matrices 
commute. 

The characteristic determinant cannot be written down in closed form for 
the general case. vVe give below its form for two particular cases: 


(1) Sj 0i j 3, . . . , N 1 



(126) 
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I 


Various generalizations of the asymmetric circulant are possible [26] . We 
mention here one of the simplest: instead of applying the permutation tt c of 
eq. (121) on the columns of the unit matrix, one applies it on the columns of an 
arbitrary diagonal matrix A = diag ( a t a 2 . . . a N ). The resulting generalized 
circulant a' n is of the form: 



Sj (0') J 


where 



The eigenvalues of A' n are then 



j=o 


and its jth eivenvector is given by the column 


(127) 


(128) 


(129) 
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1 



in which 

. 2rrj 

x j = V" S e N ; j = 0, . . . , N - 1 


(ii) Alternating Asymmetric 
We consider the matrix 



( 130 ) 


(131) 


(132) 


( 133 ) 
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where J is the same as the second matrix on the Lh.s. of eq. (27). The first 
matrix in eq. (133) is a simple circulant and we apply on a o „ the similarity 
transformation U of eq. (124). Then 


A 2n = U' 1 A 2n U = M + J ' 


( 1 ^ 


where 



(13 


and 



(12 


We further apply on A' 2n a similarity transformation P, 


8 . 


2 k ; k = 0, . . . . n - 1 


P Wi = < 


(1 


8 j . 2k-2n+i : k - n, ... , 2n 1 


such that 


P-i A' 2n P = diag(V A r . . • , A n . 1 ) 


(1 



where 



u - v 


\ 2 


u - v 


2 




r+n 


\ 

/ 


(139) 


It is of some interest to note that the permutation P of eq. (137) is a special case 
of the permutations (115) which reduce displaced Jacobi matrices. The effect of 
P on a matrix A of the type 



where the A 's are n x n diagonal matrices, is as follows: 


P' 1 A 2n P = diag(A r A 2 , ... , A n ) 

where 



The eigenvalues and eigenvectors of A r 
the eigenvalues of A are 



can be easily found. We summarize: 



r = 0, . . . , n - 1 


(140) 
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the order of appearance being K * A.“ ... ^-1 ^n-i * The eigenvectors are the 
columns of the matrix UPS, with U and P as described above, while S is given 

by 

S = diag(S 0 , Sj, • • • , S n .j) (141) 

in which S r is the 2 x 2-matrix diagonalizing the matrix A r of eq. (139). 

Analytic results, such as given above, have not yet been found for alternating 
circulants of odd order. 


(iii) Doubly Alternating Asymmetric 
Let A be the matrix 



In the sequel we shall employ the notation 
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L 



Let us define now the permutation matrix P 


Then 


k,2j ; j = 0 n - 1 


kj 


< 

w S k.2(j-n) + l ; J - n > ‘ 


2n - 1 


(143) 


/K B\ 

P' lA 2 nP = 

\ C n D n 

where A , B , C and D are simple circulant matrices given by 

n n n n 


(144) 


A n = ( S 0 s 2 s 4 • • • S 2n- 2 )cyc. 

B n = ( S 1 S 3 S 5 • • ' S 2n-l>cyc. 

C n = <^1 tr 3 °S • • • °2n-i>cyc. f 

D n = (^O ^2 °4 ' ' ' a 2n- 2 )cyc. J 

We apply now a similarity transformation U on eq. (144), 


(145) 


0 \ 



(146) 


in which U n is the unitary matrix of eq. (124). The resulting matrix is 
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(147) 


U ' 1 P' 1 a 2n P u = 



where a a is the diagonal form of A, A B that of B , etc. A further permutation Q , 
of the form exhibited in eq. (137), brings the matrix in eq. (147) to the form 


Q’ 1 U 1 P' 1 A 2N PUQ = diagCA,,, Aj A n . 1 ) 

where 



Finally the eigenvalues of A 2n are given by: 


(148) 


(149) 


± _n + \ D /K -\ dn2 


K- 




(150) 


The eigenvectors are the columns of the matrix PUQS, where 


S = diag(S 0 , Sj, . . . , S n _ t ) 


(151) 


and the ^ are 2 x 2-matrices diagonalizing the A k of eq. (149). 


(iv) Simple Symmetric 


(1) N = 2n 


A 2n ( S 0 S 1 ’ ’ • S n S n - 1 ' ' ' S l)cyc. 


(152) 
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This is a special case of the simple asymmetric matrix (116), its eigenvalues 


are 


K = s o + 2 



77k j 

s. cos 

J n 


+ S n 


(153) 


k = 0, - . . , 2n - 1 

while the eigenvectors remain the same. 


(2) N = 2n + 1 


^ 2 n+l ( S 0 S 1 ’ ' S n -1 * ' ’ S l ^ cyc . 


(154) 


^ = s ° + 2 Z] Sj 

1^1 


2vk j 

. cos ; k = 0, . . • , 2n 


(155) 


2n + 1 


Again the eigenvectors are the same as those for (166). 

For the simple case s. = 0, j = 2 there is no need to distinguish 

between the parities of N , the eigenvalues of 


A n = (* 0 s, 0 


0 


eye . 


(156) 


being given by 


K 


k 


+ 2Sj cos 


277k 

~N~ 


k = 0 N - 1 


(157) 
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The same is true for the slightly more general ease 


A N = ( S 0 S 1 S 2 0 ' • 0 S 2 S l)cyc. 

for which 


( 158 ) 


k. = s 0 + 2s. cos + 2s, cos — — 

k 0 1 N 2 N 

k = 0, . . . , N - 1. 159) 

The characteristic polynomials for (156) and (158) can be easily obtained b using 
eqs. (125) and (126). 


(v) Doubly Alternating Symmetric 


Let 


A 


2n 



S 


1 * 


n - 1 




a 

n 


a 

n - 


1 ’ ■ 



(160) 


be a special case of eq. (142). The eigenvalues and eigenvectors are easily ob- 
tained from those of (142). Below we exhibit only one particular case of (142) 
which appears more frequently: 


A 


2n 



(161) 
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The eigenvalues are 



s o +CT o 


+ 



+ sj + 2s t s 2 cos 


27rk 

n 


+ s 


2 

2 


k = 0, . . . , n - 1 


(162) 


The characteristic polynomial of A 2n in (161) is given by 

U 2n - Ml =® n M -S* S 2 J9 n 2 (\> -2( Sl s 2 ) n (163) 

where J9 n (\), J9 n 2 (\) are as in eq. (6), with a - k replaced by ( u - k) ( v - k) - 
(s 2 x - s 2 ), and b replaced by - s x s 2 . 

(vi) Skew-Circulant 
Let 



(164) 


Making use of methods similar to the ones described above we find the 
eigenvalues 


A. 


k 



j=0 


iT7(2k + l)j 

e N ; k = 0, . . 


N - 1 
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The eigenvectors are the columns of the unitary matrix U which brings a n to 
diagonal form 


U 


kj 


, ff(2j+l)k 

= _L e 1 N ; k, j= 0, . . . , N - 1 . 

i/'N 


(166) 


Here also some generalizations are possible but we shall not be further con- 
cerned with these, except for the following case of an alternating skew circulant 






U + V 

+ 

2 


2n ~ \ 7 T ( 2k+ 1 ) j 

j=l 


(169) 


The eigenvectors are columns of the matrix UPS, with U given by eq. (166), P 
by eq. (137) and S by 
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S = diag(S 0 , Sj s n . x ) 

where S k is a 2 x 2 -matrix diagonalizing the matrix A k 



( 170 ) 


( 171 ) 


We conclude this discussion of one dimensional matrices with the remark 
that no restrictions, as to relative values, were imposed on any of the parameters 
involved. Hence, all the results presented are universally valid. We shall see 
below that this is true also for higher dimensions. 
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Part n. Two Dimensional Matrices 

Introduction 

It is remarked in the text that whatever the dimensionality of the problem 
from which a matrix originates, the matrix is always an ordinary one, i.e., its 
elements are complex or real numbers and only two indices are needed to 
specify all of them. On the other hand it happens in many cases that the matrices 
considered may be handled more conveniently if some partition into blocks can 
be effected. This is particularly so when a large number of elements vanish. 

In this fashion the concept of a generalized matrix arises. Several problems 
confront us once a partition is obtained: for instance, the individual blocks may 
not commute and just as important, properties of the original matrix may not 
carry over to the partitioned matrix, e.g., symmetry in the ordinary elements 
does not necessarily imply symmetry in the matric elements. 

We exhibit below results which parallel those in Part I. It is clear that 
whenever a one-dimensional matrix can be brought to diagonal form by a 
similarity transformation which is independent of the matrix elements, the 
corresponding result for a generalized matrix can be obtained automatically, 
regardless of the commutativity of the appropriate block matrices. For other 
cases though, commutativity is essential. 


1. Continuant Matrices 


i) Simple Continuant 

Let A be the NjN 2 x N X N 2 -matrix 



(172) 
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in which F, G are N x x N x -matrices . A reduction of A to block-diagonal form 
can be effected immediately by using the similarity transformation T N2 ® I N , 
with T n ^ given by eq. (12), I N the unit matrix and the Kronecker (or direct) 


product ® being defined as follows 


(\ ® 'n.) = ( T ».) * V k. i =!.•••. N 2 


(173) 


The reduced matrix a' will have therefore the form: 


< 174 > 


Note that this reduction is independent of the commutativity of F with G . 
The matrices defined in eq. (174) may be termed the generalized eigenvalues of 
A . It is clear that the actual eigenvalues of A will be those of the N 2 Nj x N ( - 
matrices of eq. (174). Using the matrices of Part I, a whole body of results 
can be obtained for all those cases where F and G can be simultaneously diag- 
onalized. Here we discuss only the case of F and G being simple continuant: 


a n b s 

B, n n X x 0 
' \ n 

V \ N 

o N x N n N B 
n b n A 


<\\ 

D \ \ 0 

\ \ N 
\ \ \ 
o \\'d 

D N C 


(175) 


Then the transformation I ®T N , with T N as in eq. (12), brings A' to 

W 2 W 1 "l 

diagonal form. The eigenvalues of A are 


\ = A + 2B cos — + 2C cos— + 4D cos — cos ^ • - 

141 N x + 1 N 2 + 1 Nj + 1 N 2 + 1 


k = 1, . . . , Nj ; j = 1, . . . , N 2 


(175) 


The eigenvectors are the columns of the matrix ( T N ® I N ) ( I N ® T N ) 


T n ®T n . The elements of this matrix (in vector index notation) are: 
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(\*\) 


pq 


V(N 1+ 1)(N 2 + 1) N 1 +1 


q : 77 p 2 q 2 

sin sin 

N 2 + 1 


(176) 


where p = ( p t , p 2 ); q = (qj , q 2 ) and p l9 qj = 1, • • • ,N t ; p 2 , q 2 - 1, • • • , N 2 . 

The results above are still valid when A , B , C , D become arbitrary matrices. 


(ii) Alternating Continuant 
Let A be the matrix 



The matrices F t , F 2 , G are Nj x Nj -matrices, and when N 2 is odd the last 
F matrix in (177) is F t while for even N 2 it is F 2 . We have to treat these two 
cases separately. 

(1) N 2 = 2 n 2 

We follow the procedure of the one-dimensional case and apply the similarity 
transformations T 2 P 2 = (T 2n ®I N ) (P 2 n 2 ®I Nl )> wtiere T 2n2 is given by eq. (12) 
and P 2n2 by eq. (33). The result is 

P" 1 T' 1 AT, P 2 = diag - ^ ) (178) 

where 


H k = 


T, + F, 


- 2 + 2 G cos 77 k 


2n 2 +1 


F ! F 2 


F l- F 2 


F + F w(2n + 1 -k) 

- + 2G cos / (179) 

2 2n 2 + 1 / 
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If the three matrices Fj , F 2 and G commute, then the same similarity trans- 
formation S N (whenever this exists) will bring them to diagonal form, and we 
write 


S' 1 P~ x T" 1 A T 2 P 2 S = diag (A A 2 , 



where 


S = I 


2n- 



(180) 


(181) 


and 




(182) 


The Nj x Nj matrices A Fl , A*" 2 , A c are the diagonal forms of Fj, F 2 and G 
respectively. Each A k can now be brought to a block-diagonal form, the blocks 
being 2x2 matrices , 




>/ x + 
j 


77k 


i 2 cos „ , 

2 1 2n +1, 


The similarity transformation Q 2N ^ needed, is as defined in eq. (137) (there 
it is denoted by P ) . It is a simple matter now to write down the eigenvalues of 
A kj and its eigenvectors. Returning to the original matrix A , we see that its 
eigenvectors are the columns of the matrix T 2 P 2 S OR where T 2 , P 2 ,S have been 
already defined above, and 
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(184) 


Q = 



R = diag (R n ,R 


® q 2Ni 

12’ * ' * ’ ^n 2 Nj) 


(185) 


The R kj in eq. (185) are the 2 x 2-matrices which bring A kj to diagonal form. 

In certain applications (e.g., a diatomic lattice) we encounter the case where 
Fj,F 2 and G do not commute. If these matrices are entirely arbitrary no further 
reduction of A is possible. Below we treat one particular case which admits of 
a complete reduction even though commutativity does not exist: 



G = 


/C X A. 0 

\ \ 
o \ \ 

\ D 'C 


(187) 


The last A' s in F , F 2 will be A 2 ,A X respectively forN x even, andA x ,A 2 for 
N x odd. It is easily seen then that the matrix H k of eq. (179) can be written in the 


form 



(188) 


where H * are continuant matrices with elements 

k 




i - l) f 


Ai + A, 


+ 2 C cos 


7T k 


2n 2+ l 


S . 

r 3 


(189) 
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and 


J = diag 



- A 2 


2 


*1 




-a 2 


2 


) 


(190) 


Since the H* ’s are continuant, the transformation T N ^ of eq. (12) will bring 
them simultaneously to diagonal form. Also we have already shown how the J 
matrix transforms under T in eq. (30) . Hence on applying the transformation 
T x = I 2n ® T n on the matrix of eq. (178), we obtain 

T; 1 P' 1 T” 1 A T 2 p 2 T x = diag , H' 2 , • • * , H'„ 2 ) (191) 



The diagonal matrices contain the eigenvalues of . We employ now 
the permutation matrix P x = I nj ® P 2Ni » with P 2N ^ defined in eq. (33), to bring 
H'to a block-diagonal matrix, the blocks being the 2 x 2-matrices A' k . 



with 


a 1+ a 2 

^kj = — t ~ 


± 2 C cos 


7rk 

2 r» 2 + 1 


+ 2 


^ B ± 2 D cos 


77k \ 77 j 

2^TT) COS vT (194) 


Finally the eigenvalues of the original matrix A are given by 


\*. = 
kj 


^kj +^kj 




^1 ~^2 \ 2 


k = 1 , • • • , n 2 ; j = 1 , • • • , Nj 


(195) 
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The eigenvectors are the columns of the matrix T 2 P 2 Tj P x S , where 
T 9 , P, , T, , Pj have already been defined and Sis the matrix 

S = diag <S n ,S 12 , * ’ * , S n oNl ) ( 196 ) 

the S. . being the 2 x 2-matrices diagonalizing A/ . e( b (193). 

k j J 


(2) N 2 = 2n 2 + 1 

It is readily shown, in a manner identical to the case of N 2 even, that the 
eigenvalues of A for non-commuting F lt F 2 , G are 



The other Nj eigenvalues are those of the matrix F t , which for Nj even are 
given by eq. (36) and for Nj odd by eq. (42), when u, v are replaced by Aj, A 2 
respectively. 

Trivial changes in the transformation matrices of the even N 2 -case will 
quickly yield the eigenvectors for the present case. 

Similarly, for commuting , F 2 and G , the eigenvalues of a will be the 
roots of the matrices 



k = 1, • • • ,n 2 ; j = 1, • • • ,N 1 
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To these roots one has to add the eigenvalues of Fj in order to obtain the 
complete spectrum of A . 

We conclude this section with the following remark: the reductions performed 
in the non-commuting case will be still valid to a certain extent if the numbers 

A, B, C and D are replaced by matrices. Thus if these are arbitrary matrices 
eq. (193) remains valid, further reduction being dependent on the nature of A, 

B, C, D. We shall make use of this result in our treatment of higher dimensional 
matrices. Here we mention one case arising in a two-dimensional lattice problem: 
A, B, C and D are 2 x 2-diagonal matrices. Then eqs. (195), (197) are valid if 
these matrices are replaced by their respective diagonal elements. 

(iii) Variants of the Continuant 
(1) Let A be the matrix 



and F, G are N x x -matrices. 

We define now a similarity transformation T 2 = T N ® I N , where T N is the 
orthogonal transformation of eq. (65). Then 

T" 1 A T 2 = diag (201) 

with 

H. = F + 2 G cos — ; k = 0, • • • ,N, - 1 (202) 

N 2 * 

This reduction is valid for arbitrary F, G . Further reduction will depend on 
the nature of F, G. If in particular F andG commute, then the same similarity 
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transformation S will bring them to diagonal form and the eigenvalues of A 


will be 


k = 0, • • • ,N 2 - 1 


\ = \? + 2 A9 cos 

> J N 

j = 1, • * ‘ .Nj 

where and are the eigenvalues of F and G respectively. 
We make use of eq. (203) for the special case 


/ A 


\ 0 

\ 


+ B B 
B A 

F - I \ v \ 

r - 1 \ \ \ 

0 \ A B 
\ 

B A + B 


; G = 


c 

\ 

\ 

\ 

\ 

0 \ 


(203) 


(204) 


Then we can identify S with T from eq. (65), and the eigenvalues of A will 

N i N i 

be given by: 


k = 0, • • • ,N, -1 


7T 1 „ „ 77 k 

. -A + 2B cos — -*-+ 2 C cos — — 
N, N, 


(205) 


j =0, • • • ,Nj-l 

The eigenvectors are the columns of the matrix 

* - ® V Un, ® V = t k, ® Tk. 

Note that the result in eq. (205) remains valid even when A, B and C become 
arbitrary matrices. 


(2) Here we consider the alternating matrix 



(206) 
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Again F x , F 2 and G are Nj x N x -matrices. If we employ the results for 
the one-dimensional case (85), we can show — provided the matrices Fj , F 2 , G 
commute and G possesses an inverse — that the eigenvalues of A are the roots 
of the equations 

(F, - \I) (F 2 - M) -4G 2 cos 2 — = 0 (207) 

N 2 

k = 1, • • • ,n 2 - 1 

and 

(F t - A I) (F 2 - A I) + G EF 1 - M + F 2 - AI] = 0 (208) 


The eigenvectors will follow just as in the one-dimensional case. Since by 
assumption F ls F 2 and G commute, the eigenvalues of A will be the roots of 
the equation 


7Tk 


(A 1 - k) (\ 2 - k) - 4 (X?) 2 cos 2 _ = 0 
j i 1 n„ 


k = 1, * • • ,n 2 - 1 


j = 1, • • • ,Nj 


> (209) 


(\ f i _ K) (k * 2 - k) + \? r k Fi - \ + \ F2 - \i =o 

j j J L j j J 

The more interesting case of non-commuting F t , F 2 where 


J 


/ 


A I+ B \ 


B A 2 

\ \ 

\ \ 

\ \ 

0 \ \ 


\ 


\ 0 
\ 

\ 


B 


y b \ 2 + b / 


A 2 + B B 
B A, 


‘ \ 


; F„ 


\ A 

\ A 2 
\ \ 

0 \ \ 

\ \ 


' o 
\ 

s 


B 


B 'Aj + B / 


G = 



\ 


\ 0 


V 


(210) 


has up to the present resisted all efforts toward its diagonalization. 
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2. Circulant Matrices 


(i) Simple Asymmetric 

Here we consider the matrix 

* = (*. A. ••••*„,- (211) 

in which the A k are N x xN ( -matrices. 

The similarity transformation U 2 = U N ® I N , with U N given by eq. (124), 
reduces A to the form 


Uj 1 A U 2 = diag (A 0 , Aj, 



( 212 ) 


where 


n 2 -i 

A r = J2 A i ei ^T 1 : r = 0, •••’N 2 -l (213) 

j =o 2 

Further reductions will depend on the nature of A^ If for instance all A j 
commute then the transformation S = L. ® —where S N brings the A. 

IN 2 N j IN j ) 

simultaneously to diagonal form - will completely diagonalize A : 

S" 1 ^ 1 AU 2 S = diag(X 00 \ 01 ••• ^ N2 _ 1iNi ) (214) 

in which 


n 2 -i 




j =o 


X A i e i iZLli ; 
k N„ 


k = l, 


r = 0, • • * ,N 2 - 1 


(215) 


The eigenvectors are the columns of the matrix II ® S„ . 

M 2 Nj 

The special case in which the matrices A v are circulant is of particular 

N 2 

interest: the matrix S N is then the U N of eq. (124) and the eigenvalues are 

N 2 -1 Nj-l 


\ 


rk 




i 2 77 k ^ i 2 77 r i 

U N. N 


(216) 
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The result in eq. (216) remains valid when the A j i become arbitrary matrices 
since the diagonalizing transformation does not depend on the matric elements 
of A . 


(ii) Doubly Alternating Asymmetric 


The matrix to be considered here is 


A = 


A 0 A x 


®2n 2 - 1 ®1 * * * 


l 2n 2 -1 


(217) 


B 


2n 2 2 / eye. 


We have already used this notation in eq. (142). Here the A’s and B’s are 
NjXNj -matrices. We define the similarity transformation P 2 = P 2 „ 2 ® I„ 
with P 2n given by eq. (143) . Then 

'a(!> A< 2 > 


p^ap 2 = 

a< 3 > a< 4 > 

where A ( r ) are the generalized simple circulants 

A^ 1 ^ = (A 0 A 2 • * * A 2n ^_ 2 ) cyc 

A<2) = ( A 1 A 3 **• A 2n 2 -l)eyc. 

A (3) = (B X B 3 ••• B ) cyc . 


(218) 


> 


A (4) = (B 0 B 2 • • * B 2n 2 ) ( 


eye . 


(219) 


We apply now on eq. (218) the similarity transformation U 2 = I 2 ® U 2 
where U 2 is as in eq. (212). The result will be 


Up PI 1 A P 2 U' = 


/ K (l) 

K< 2 > 

\ K (3) 

k< 4 > 


( 220 ) 
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where K (r) is the block-diagonal form of A (r) , r = 1, 2, 3, 4. A further permuta- 
tion Pj = 0 2n ® I N » with 0 2n as in eq. (137), brings the matrix in eq. (220) to 
the form 


p-i 


u- 


P^AP, U' Pj = diag (H ( 


o’ ^1’ 


H _,) 
n 2 l 


( 221 ) 


The matrices H k are given by 



( 222 ) 


and K[ r ^, r = 1, • • • , 4; k = 0, • • • ,n 2 - 1 is the k th block of K (r) , similar in 
form to A k of eq. (213). 

To reduce the H k we need at this point to know the commutation properties 
of the Vs and B's. Assume first that all these matrices commute: then, if IL 

— w i 

denotes their common diagonalizing matrix, we obtain the matrix H k 




A ( k X) 

Ai 2 >\ 

A ( k 3) 

Ai 4 V 


(223) 


where A ( k r) is the diagonal form of K ( k r) . Now, the permutation P 2N defined in 
eq. (137), when applied on leads to the matrix 


P 2 Nl H k P 2 Nl -diag(J ki , J^, ••*,J kNi ) 


where 


'kr 


Aiv m, 2 A k = "2-i 

vtv k iv) ' -- 1 ", 


> 


(224) 


(225) 
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Finally the eigenvalues of the original matrix a are 


k r _ 


^kV - w* 


+ >S» X< k >> 


(226) 


k = 0, . • • , n 2 - 1; r = 1, . . Nj 


The eigenvectors are the columns of the product of the various transformations 
employed. 

We return now to the original matrix a and assume that not all of the A 's 
and B's commute. For simplicity we shall assume that Aj , Aj, • • •, A 2n ^. 1 
coincide with B,, B„ » • • • » B„ , and all are simple circulant matrices, while A 

1 2 2n ^ - 1 u 

and B 0 are alternating circulant of the following type: 

A)1 A)2 A>N,-1 


A 0 = 


A)1 A)2 \nj-2 j 


eye. 


B„ 


@ \l A) 2 ■ • ^Nj-l 


y^oNj-i a An An ■ ■ A)Nj 


eye . 


(227) 


Let us denote by Aq the generalized simple circulant 

\ +B 0 


*0 - ( ~ ' ^1’ ^2’ ’ " ’ ^2n 2 -l 


eye. 


(228) 


Then the matrix A can be written also as 


A- — B 0 _ . 

A = Ap + — diag (I Nj . - I Ni I Ni - - I N 1 > 2 n 2 ' 


(229) 


Then the transformation U 2 of eq. (212), with N 2 = 2n 2 , reduces A Q to block- 
diagonal form, and on applying it on a it yields 


U- 1 AU 2 =diag(A<°>, A<*> A (2na ' l> )+J 


(230) 


225 


in which 



A further reduction of the matrix in eq. (230) is provided by the permutation 

P 2 = P 2 n 2 ® J Nl * With P 2n 2 & Ven ^ eC b ( 137 )> 

P 2 J u 1 aU 2 p 2 = dia S(H 0 - H r .... H^.j) (233) 


where 

A ( k) B o 

2 

(234) 

\ ~ «0 (k+n 2 ) 

u 

2 / 
k = 0, . . . , n 2 - 1 

At this point we must assume n = 2n x since no analytic solution is known for 
the odd case. Hence, if this is so, we can write: 



Aq-Bq 

2 


diag (1, -1, . . . , 1, - 1), 


(235) 


On the other hand the A (k) ’s are simple circulant and can be brought simultaneously 
to diagonal form by the transformation U defined in eq. (124) . If we apply on 

2n i 

H, the transformation T ® U„ we obtain 

k l 2n 


226 


(I 2 ® l! 2, 1] )‘ l "kdj •«,„,) = diag(A< 0 k >, A ( j k \ A < 1 k *' V ) + 

0 i„,\ 


0.-/3 



(236) 


in which A ( Q r) contains the first n x eigenvalues of A ( r > and A (r) the n t suc- 
ceeding ones. 

The similarity permutation I ® P. , with 

n x 

I 1 0 0 0 ^ 

(237) 

\o 1 0 0 / 

brings the matrix of eq. (236) to the form 
a ni ® « k (I 2 ®U 2ni )(I ni ©p 4 ) = diag(H( k °), H^ 1 ) ) (238) 


0 0 10 
0 0 0 1 


where 


= 


r A< k) —I 

0 2 i 


0-/3 T Ck + n 2 ) 


/ A < k ) ^ I 

/ 1 2 "il 




I A. 

2 n i 1 


V^ 1 


Tk+n 2 ) 


(239) 


n, 0 


We apply now the permutation P , defined in eq. (137), on Hi 0) and Hi 1 ) to 

2n k K 


obtain 


K. »£ 0> P 2 „. = diag (If" , *<„»>,) 


"\ 


> 




(240) 
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where 


1 


1 

/ X (k) 
j +n l 

a - /3\ 

2 \ 

II 

o 


P* 

■ *-* 

II 

l a - fi 

\ 2 

(k+ n 2 ) 1 
' / 


In eq. (241) A< r) is b y definition the i t h eigenvalue of the matrix A (r) of 
eq. (231). Again the entire reduction up to this point remains valid if the A.j r >'s 
become matrices. It is an easy matter now to write down the eigenvalues of A, 
after solving the quadratic equations arising from K| 0) and K: 1 \ The eigenvectors 
will be the columns of the matrix U 2 P 2 Uj 0 Pj S , where U 2 is defined in eq. 

(230), P 2 in eq. (233), U, - (I 2lIj ® U 2n ) and U 2iij as in eq. (124), Q - (1^^ ®P 4 ), 
Pj ■ (I 2n ® P 2n ) and P 2n ^ as in eq. (137), and finally S = diag (S<°> , S<°>, 
eU) ), in which S (r) is a 2 x 2-matrix diagonalizing K. (4) , r = 0,1. 

n 2 “1 , n,-r kj XJ 

(iii) Symmetric Circulants 

As in the one-dimensional case, the results for the symmetric circulants 
follow automatically from those for the asymmetric ones and therefore will not 
be treated further here. 

Part III. Three Dimensional Matrices 

The results to be obtained in the following depend largely on those of Parts 
I and II. Since we have shown in Part II how partial diagonalizations can be ef- 
fected in several cases involving arbitrary matrix elements, here we shall pro- 
ceed directly to the matrices of interest. 

1. Continuants 

(i) Simple Continuant 

Let A be the N 2 N 3 matrix 
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0 


with 




(242) 


(243) 


(See page 230 for Equation 243a) 


All of the matrices above are generalized or ordinary simple continuant 


matrices, and we can write down the eigenvalues and eigenvectors immediately: 
= A + 2Bj cos j 6 + 2B 2 cos kcp + 2B 3 cos r\p 


+ 4Cj cos j 0 cos kcp + 4C 3 cos j 6 cos n p + 4C g cos kcp cos r\p 
+ 8Dj cos j 6 cos kcp cos n p 


> 

(244) 


n 7T 77 , TT 

d = - : <p=- ; ; »/,= 


Nj + 1 


n 2 + 1 


n 3 + 1 


j = 1 , . • . , NjJ k = 1, - . . , N 2 ; r = 1 , . . • , N 3 


J 


This result remains valid if A, B. , C. and Dj become arbitrary matrices. 
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CO 

TP 

<N 



II 




II 

S3 



m m 

/C J CQ 



II 
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The eigenvectors of A are the columns of the matrix T N ® T M ® T w , 

w 2 "3 N 1 

where T N are given by eq. (12). 

(ii) Alternating Continuant 

The matrix to be considered here is of the form 


where 


A = 




\ 

\ 

\ 


\ \ 


\ \ 


0 


T 




( 245 ) 


S. 



N 2 = 2n 2 + 1 
N 2 = 2n 2 


( 246 ) 


Also 



( 247 ) 


( 248 ) 
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and 


I \ \ 

B 1 A 2 \ 


\ 


\ 


\ B 1 

\ 

' \ , 

B I A '/ 


A’ = <^ 


( A x ; Nj = 2n x + 1 

A 2 ; N i = 2n i 


= F 2 = 


A 2 B t 

\ 

Bj Aj \ 0 

\ N \ 

\ \ X 

\ \ B. 

° \ \ 1 


\ 


\ \ / 

B i A"/ 


A 2 ; Nj = 2nj + 1 


A" =• 






Aj ; Nj = 2n x 


(249) 


(250) 


The matrices T and G are as in eq. (243) . 

For simplicity we shall choose N 2 = 2n 2 , N 3 = 2n 3 and Nj = 2n r All 
other possible choices lead to soluble cases and the procedure to be followed is 
similar to the one used below. 

First note that (Sj , S 2 ) and (Fj, F 2 ) are pairs of non-commuting matrices. 
We proceed as in the two-dimensional case, and apply successively the trans- 
formations T 2 P 2 , 


where 


P^T" 1 AT 2 P 2 =diag(Q 1 , Q 2 Q n > 


T 2 _ (T 2 n 2 ® V ® \ ’ P 2 “ < ' P 2 n 2 ® ® 1 


with T 2n , P 2n as in eq. (178), 



and 


Si + S 

A k = — - + 2T cos kcp; 


9 = 


Ti 

2n 2 + 1 


(251) 


(252) 


(253) 


(254 


232 



Recalling the definitions of the various matrices in eq. (254), it is easily 
seen that the ’s are generalized continuants the matrix elements of which are 

( Vpq = ( ~ 2 + 2Hcosk ?) s pq + (G + 2K cos kcp) (8 ptq+1 + s p q-1 ) (255) 

and 

(256) 



We apply now successively the transformations T = ( I ® T, ) ® I. and 

P 3 = ^2n ® P 2n ) ® 1 2n — With T 2n &S in eq - ( 12 ) End P 2n aS 111 ^ ( 33 ) _ ° n 

eq. (251) to obtain 


<T 2 P 2 T 3 “V* * T 2 P 2 T 3 P 3 = diag (Q u . Q, 


2 y 


* ‘ ’ a n 2 .2n 3 ^ 


in which 



and 


(257) 


(258) 


F + F 

= -i - + 2H cos iy + 2 [G + 2K cos Zy] cos mi p 


(259) 


'C — 1 , • . . , 2n , m — 1, . « .f 2n , i p — 


77 


2o 3 + 1 


Again we observe that the A^ m of eq. (258) are continuants with the elements 

fA, +A, 1 

= j 2 + 2B 2 COS + 2 [B 3 + 2C S COS ^ COS } 8 pq + 


+ {Bj + 2Cj cos <ty + 2 [C 3 + 2D X cos ty] cos mi//} (S 


p,q + l 


p.q 




p. q = 1 , • • • , 2n t 


(260) 
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If we apply on eq. (257) the successive transformations T = ( I ® I ) ® T 0 , 

X zn 2 3 ^ ri i 

P. = (I, ® I. ) ® P, , with T and P as in eq.'s (12) and (33) respectively, 

1 zn 2 Zn ^ *n j Z n ^ z n ^ 

we get 


< T 2 P 2 T 3 P 3 T , P .)'' A^PjTjPjT.P, = diag (a ln , n 112 


> (262) 


where 


krj A !- A 2 


A 1 “ A 2 


2n 2~ k +1 , an^r + l^rij-j +1 


(263) 


A + A 

A^ mn = — i- — — - + 2Bj cos ni9 + 2B 2 cos *tcp + 2B 3 cos nup 


+ 4Cj cos nd cos 'ttp + 4C 3 cos nd cos m <p + 4C S co s 'tcp cos m0 


8Dj cos ni9 cos •Cep cos mi p 


- 1 , • • • , 2n 2 , m = 1 , . • . 2n 3 ; n = 1 , . . • , 2^; 6 - — - ^ (264) 


It is easy now to find the eigenvalues of S2 kr j , provided that all capital 
letters in eq. (263) and (264) denote either scalars or diagonal matrices. Thus 





in which // kr is identical with A krj of eq. (264) and /uT kr has the same form as 
M j kr but each cosine enters witha(— ) sign. 

(iii) A Simple Variant 
Let 


A = 


/s+T T x 

T S \ 0 

\ \ ' 

\ \ ' 


\ 




\ \ 

S T 


\ 


T S+T 


I 


where 


(266) 



Just as for the case (i) above the eigenvalues can be written down im- 
mediately: 
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\ = A + 2Bj cos j 6 + 2B 2 cos kcp + 2B 3 cos r\p 

+ 4C X cos j 9 cos kcp + 4C 3 cos j 9 cos r\ p + 4C S cos kcp cos r\p 

+ 8D X cos j 9 cos kcp cos r\p 


n 7T 77 , 77 *A 

9 = tr ; ¥ = n : ; ' =0, 

‘l 2 3 


,Nj-l; k=0 Nj-1; r = 0, ...,N 3 -1 


(268) 

The eigenvectors are the columns of the matrix T N ® T N ^ ® T N ^ , with T N given 
by eq. (65). 


2. Circulants 

(i) Simple Asymmetric 

It is clear from the previous treatment that general results can be 
written down for arbitrary three-dimensional circulants, the elements and sub- 
elements of which are also arbitrary circulants. For the sake of simplicity we 
restrict the following discussion to the simpler and more frequently met matrix 


A = 




\ ^ 
\ 

0 \ 


\ 

\ \ 

\ \ 


\ 


\ 


\ 


\ 


\ 


T 



(269) 
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in which 


S = 




G 





H K 

\ \ 

L \ x 0 


\ \ 

\ \ 


\ 


\ 






' B 3 C 3 C 4 

\ \ 

c ;\ \° 

« = i ' \ \ 

0 \ \ c, 
c < x B =, 


\ 


\ 





I (270) 




We assume A to have ordinary symmetry* (in the elements A, B,C,D) and 
therefore the transposed block matrices in eq. (269) and eq. (270) are transposed 
as ordinary matrices also. This means, for instance, that 


T = 


H L 

K \ \ 
\ \ 


\ 


\ 


\ \ L 



\ 


k 



(271) 


*The asymmetry in the title of this section refers to the form of A as a block matrix. 
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with similar expressions for H, K and L if these in turn are generalized 
matrices. 

The eigenvalues of a are 


+ 

+ 


A jkr = A + 2Bj cos j 0 + B 2 e ikC P + B 2 e' ik( P + B 3 e _i ^ + § 3 e _i ^ + Cj e^i***” 
r e -i(j5+kcp) c 0 i (kcp-j 0) + q e -i(k<p-j0) + c e i(i0 + W') + C, 

2 2 3 ^ 

q e -iC jfl-ri/') + q e i(j0 " r ' / ' ) +C e i ( kcp+r ' / ' ) + C, e ' i < k(p+r ^ ) + 

4 4 5 5 

+ C e i(kC(>-r<//) + Q 0 - i ( kCp- r-A ) + Q e i ( j -9+kCp+nA) + e _i < > < 9+ k<P+r>/') + 

6 6 1 ^ 

, d e " i ^ 6 ’" kCp+r ' / ' ) +D e i (j |9 " k< l )+r ^> +D e i ("j^ +k( ^ +r ^'l + D e i(je_kCp " r ' / ' ) + 

2 2 3 3 

+ D 4 e^^+kV-WO + D 4 e -i(j<9+kCp-r-/>) (272) 


where 


2tt 2tt , 277 

(9 = — ; cp = -rj— - 'P = t: — 

N, V N 2 N 3 


■x 


(273) 


j=0 Nj - 1 ; k = 0 ,N 2 -1; r=0,...,N 3 -l 


The eigenvalues \ jkr in eq. (272) were written under the assumption that 
A, Bj » • • • » D 4 were matrices. If we assume these to be either diagonal 
or symmetric, we obtain the simplified expression 


A = A + 2Bj cos j 9 + 2B 2 cos kcp + 2B 3 cos r p + 

+ 2Cj cos (j 6 + kcp) + 2C 2 cos (j 6 - kcp) + 2C 3 cos (j 9 + xp) + 

+ 2C 4 cos (j 9 - xp) + 2C 5 cos (kcp + xp) + 2C g cos (kcp - xp) + 

+ 2Dj cos (j 9 + kcp + up) + 2D 2 cos (j 9 -kcp + up) + 2D 3 cos (j 9 -kcp - rp) 

(274) 

+ 2D 4 co s ( j 9 + kcp - r\p) 
with 9 , cp, p and j , k , r as in eq. (273). 

The eigenvectors will be the columns of the matrix U M ® IL ® U with 

N 2 N 3 1 

U ’s as in eq. (124). 


f 
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ii) Alternating Asymmetric 

The matrix considered here is of the form 



with Sj , S 2 regularly alternating along the main diagonal of A and 



( 275 ) 


(276) 


(277) 


Fj and Fj alternate regularly along the main diagonals of , S 2 and so do A } 
and A 2 in F x , F 2 . The matrices T and G are as defined in eq. (270). 


We remark that the matrices Sj , Fj do not commute with S 2 , F 2 , respectively. 
To diagonalize a we put it in the form 


A = A 0 + diag 
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in which ±1/2 ( Sj - S 2 ) alternate regularly and 




S, + s 2 


T, 0 . . • 0 T 


eye . 


(279) 


The transformation U 2 which reduces A Q to block-diagonal form is given by 


^2 “ («N 2 ® ® J Nj 


(280) 


with IL as in eq. (124). Then we can write 
N 2 

^AU 2 = diag (A<°), A (1 >, •••• + J (281) 

where 

A 0O =L(S i + S 2 ) + Te ik( P + Te' ikC P 
2 

k=0 2n 2 -l; cp = 7r/n 2 (282) 

and 

■ ■ C, ':) 


with 


J, = i„ 2 ®j< s . - V 

A further reduction is provided by the permutation 

P 2 = (P N2 ® I N ; ® I Nj 

with P given in eq. (137) (for n = n ). The result of this permutation 


(284) 


(285) 


is 


(U, 


p y 


A U 2 P 2 = diag (H 0 , Hj, 


, H -i ) 
n 2 1 


(286) 
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in which 


4<‘> i(S,-S 2 )' 


A = 


k = 0, . 




A (k+ "2 ) 


• . n 2 - 1 


(287) 


Now, the matrix 1/2 (S x - S 2 ) is of the form 


y (Sj - S 2 ) = diag (*1 


F, - F 2 F x - F 2 


F 1 F 2 


(288) 


where ±1/2 (F x - F 2 ) alternate regularly, while the are all generalized 

simple circulants. Therefore is we apply on 1^ the transformation V x , 


v x =(I 2 ®U N )®I 


(289) 


with Ujj as in eq. (124), we obtain 


V; 1 V x = diag (A< k >, A< k >, A< k+ "2\ A ( x k+n2> ) + J' 


(290) 


in which a^ contains the first n 3 (generalized) eigenvalues of A ( ^ ) and A<^> 
the n 3 succeeding ones. The matrix J' is given by 


/ . ,: \ 


J' 


j; 


j; 


(291) 


\j; 


with 


K = I„ 3 < F , - 


(292) 


The similarity transformation V 2 = I n ® P 4 , with P 4 as defined in eq. (237), 
brings the matrix in eq. (290) to the form 


(V x V 2 r x H, V X V 2 = diag (H<°> , H<»>) 


(293) 
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where 




( 0 ) _ 


'<> j; \ 


/ A< . 1> j; 


• H (1) = 

j; a“-»»/ ‘ Vi; /C’' 


We apply now the permutation 


P, = P, ® I N 

3 2n 3 Nj 


with P 2n as in eq. (137), on H£ 0) and H£ x) : 


P-' H<«>P 3 =dia 6 (KW, K<f> 


p ; 1 H £ ,)p 3 =<ii.g(K(j). k<;> 


in which 


K <°> = 

k r 


\ r ^ F i- F P 


- (F. - F,) A. + 

2 1 2 k+n 2> r+n 3 


• K< x > - 

’ "*kr “ 


A. . -(F. -F,) 

k , r +n ^ 2 1 2 




and where 


with 


k = 0, . . . , n 2 - 1 ; r = 0, . 


• - n 3 - 1 


A = 1 (F, + FJ + H e ip( P + H e“ ipC P + G e 1 ^ 

pq 2 1 2 

+ G e~ iq ^ + K e i ^ p( P +q ^^ + K e " i ^ pC P +q ^^ 

+ L + L e" i ^ p< ^ ql ^^ 


cp = 77 '/ n 2 ; \p - 77/ n 3 ; p = 0 , . . . , 2n 2 -1; q = 0 , . . . , 2n 3 

The matrices A are again circulant, while 

pq 0 > 




\ - A 2 


, . . . , 


(294) 


(295) 


(296) 


(297) 


(298) 


_ 1 


(299) 
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in which ± 1/2 (Aj - A 2 ) alternate regularly. Then if we apply on 
the transformation 

°3 = l 2 ® % 

we obtain 


U ' 1 K<°> U, = diag (A<°> A* 1 ) A<°> . A<V . ) + W 

3 kr 3 kr kr k+r^.r+n^ k+Oj.r+n^ 

where A (0) contains the n, first eigenvalues of A , and A (1) the n, 
pq 1 0 pq pq 1 

ones. Also 


W 


/ 0 w '\ 


W. 


w. 


\ w > 


* W 1 -V 


/ 


The similarity permutation V 3 = I n ® P 4 , with P^ as in eq. (237) 
the matrix of eq. (301) to the form 

(U3V3)-' K<;)U 3 V 3 =dia g (X(;), YW) 

in which 


A(°> 

k r 


W 1 \ / W *1 


X<°) = 

kr 


• y<°> - 

kr _ 


W 1 A k+n 2 .r+n 3/ 


W A (1) 

"l A k + n 2 ,r + n 3/ 

Applying the permutation P 2n defined in eq. (137), we finally obtain 


P : 1 X<°> P 0 = diag (Q( x ’°) , Q< x -°) 

2n j kr 2^ 6 v 0 kr lkr "j-lkr 


P;i Y(°>P,„ = dia B (Q<f°>. 


2nj kr 2rij 


Ok r lkr 


n j - 1 kr 


> and Kj 1 > 

k r 

(300) 

(301) 

succeeding 

(302) 

, brings 

(303) 

(304) 

(305) 
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in which 


i 

i A. 

J *r 

y (\ - V ' 

1 1 

' "j+njkr 


= 

i \ <*. - v 


> 

/ — > 
<< 

© 

> — ✓ 

II 



jkr 

A j+nj k +n 2 r+n 3 ^ 

j jkr 1 

A j k+n j r+n 3 


j = o, . . . , 

rij - 1 ; k = 0, . 

* • • j *^2 — 1 j 

•-1 

II 

© 

n 3 - 1 ( 2 


and the expressions A{ mn coincide with the ones given in eq. (272), if A there is 
replaced by 1/2 (A x + A 2 ). 

The same procedure gives for the matrices 


n (x,D = 

jkr 


A jkr+n 3 

- A 2>\ 

JO 

►-» 

V-' 

II 

' +nj k r+n 3 

j - v\ 



jkr 



U<*, -V 

^j+nj k+n 2 r j 

1 

l j(*. -V 

^jk+n 2 r j 


(307) 


This is the reduction required. It is simple now to find the eigenvalues if the 
matrices Aj , . . . , B 4 are diagonal or can be simultaneously diagonalized. 


♦ 


4 
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